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The  Entrainment  and  Homogenization  of  Tracers 
within  the  Cyclonic  Gulf  Stream  Recirculation  Gyre 

Robert  $.  Pickart 

The  various  distributions  of  tracer  associated  with  the  Northern  Recir¬ 
culation  Gyre  of  the  Gulf  Stream  (NRG)  are  studied  to  try  to  obtain  informa¬ 
tion  about  the  flow.  An  advective-diffusive  numerical  model  is  implemented 
to  aid  in  the  investigation.  The  model  is  composed  of  a  gyre  adjacent  to  a 
boundary  current  in  which  a  source  of  tracer  is  specified  at  the  upstream 
edge  of  the  current.  This  set  up  attempts  to  simulate  the  lateral  transfer 
of  properties  from  the  Deep  Western  Boundary  Current  (DWBC)  to  the  NRG  in 
the  region  where  the  two  flows  are  in  close  contact  west  of  the  Grand  Banks. 

The  results  of  the  model  are  analyzed  in  some  detail.  Tracer  is  en¬ 
trained  into  the  gyre  as  a  plume  which  extends  from  the  boundary  current  and 
spirals  across  streamlines  toward  the  gyre  center.  The  maintenance  of  the 
spiral  during  spin-up  and  its  relationship  to  the  occurrence  of  homogeniza¬ 
tion  at  steady  state  is  examined.  An  asymmetry  in  the  spiral  exists  due  to 
the  ellipticity  of  the  gyre,  which  also  effects  homogenization. 

The  anomalous  properties  that  are  fluxed  into  the  NRG  include  salt, 
oxygen,  and  freon.  These  particular  tracers  are  independent  from  each  other, 
the  former  two  because  they  are  characterized  by  different  vertical  profiles 
in  the  deep  layer.  This  results  in  a  decay  of  oxygen  but  not  salt,  due  to 
the  presence  of  vertical  mixing  as  discussed  by  Hogg  et  al .  (1986,  Deep-Sea 
Research,  33,  1139-1165).  Their  analysis  is  expanded  upon  here.  The  effect 
of  vertical  mixing  on  the  gyre/boundary  current  system  is  examined  within  the 
context  of  the  numerical  model.  Results  are  applied  to  recently  collected 
water  sample  data  from  the  region  which  leads  to  an  estimate  of  the  lateral 
and  vertical  eddy  diffusion  coefficients  and  an  estimate  of  the  amount  of 
oxygen  in  the  NRG  that  has  diffused  from  the  DWBC. 

The  accumulation  of  freon  within  the  NRG  is  considered  in  addition  to 
salt  and  oxygen.  Appreciable  levels  of  freon  have  been  present  in  the  ocean 
only  since  1950,  and  the  atmospheric  source  functions  have  been  increasing 
steadily  since  then.  A  simple  overflow  model  is  presented  of  the  manner  in 
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which  freon  may  be  stirred  In  the  Norwegian-Greenland  basin  prior  to  over¬ 
flowing  and  entering  the  DWBC.  Once  In  the  boundary  current  the  concentra¬ 
tions  are  diluted  by  way  of  mixing  with  surrounding  water.  Two  different 
schemes  are  considered  in  which  the  immediate  surrounding  water  accumulates 
a  substantial  amount  of  freon  as  time  progresses.  These  models  suggest  that 
the  freon-11 :freon- 12  ratio  may  not  be  a  conserved  quantity  for  the  water  in 
the  core  of  the  l)W8C.  It  is  found  that  the  level  of  freon  in  the  NRG  is 
barely  above  the  existing  background  level. 
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Introduction 


Over  the  past  several  years  a  clearer  picture  has  emerged  of  the  mean 
abyssal  circulation  of  the  western  North  Atlantic.  In  particular  the  exis¬ 
tence  of  a  deep  cyclonic  gyre  situated  between  the  New  England  seamounts  and 
Grand  Banks,  whose  eastward  flowing  portion  contributes  to  the  deep  Gulf 
Stream,  has  been  addressed  by  Hogg  In  a  series  of  works.  Tnis  is  in  contrast 
to  the  deep  anticyclonic  subtropical  gyre  described  by  Worthington  (1976). 

Worthington  postulated  using  a  historical  data  set  that  roughly  60 
Sverdrups  of  water  colder  than  4°C  recirculate  to  the  south  of  the  Gulf 
Stream.  In  order  to  produce  a  consistent  flow  pattern,  the  constraint  of 
geostrophy  was  relaxed  in  certain  regions  of  the  gyre  where  instead  the 
continuity  of  water  properties  was  used  as  a  guide  for  flow  lines.  Several 
years  later  Wunsch  and  Grant  (1982),  with  the  same  data  set,  produced  a  very 
different  deep  flow  pattern  using  Inverse  methods.  They  postulated  a  cyclon¬ 
ic  gyre  transporting  roughly  25  Sverdrups  (which  did  satisfy  geostrophy 
everywhere).  In  addition  to  this  evidence  based  on  hydrographic  data,  there 
Is  also  direct  evidence  for  cyclonic  recirculation. 

In  1983  Hogg  compiled  all  the  available  deep,  long-term  current  meter 
measurements  in  this  region  of  the  North  Atlantic.  From  this  data  he  infer¬ 
red  a  consistent  streamline  pattern  which  includes  a  cyclonic  gyre  that  he 
estimates  to  be  carrying  ~20  Sverdrups.  It  is  somewhat  different  from  the 
Wunsch  and  Grant  gyre  however;  it  is  smaller  In  extent  both  zonal ly  and  mer¬ 
idionally.  Also  in  contrast  to  the  Wunsch  and  Grant  pattern  is  the  presence 
of  a  small  anticyclonic  gyre  just  to  the  south  of  the  cyclonic  recirculation. 
Hogg  estimates  this  flow  to  be  10  Sverdrups. 

Recently  additional  current  meter  data  has  become  available  in  this 
region  of  the  Gulf  Stream.  Hogg  et  al.  (1986)  incorporated  these  measure¬ 
ments  into  the  existing  historical  data  set,  and  refined  the  earlier  picture 
presented  by  Hogg  (1983).  In  particular,  the  north-south  length  scale  of  the 
cyclonic  recirculation  appears  to  be  even  smaller  yet.  A  schematic  stream¬ 
line  of  the  gyre,  which  Hogg  et  al .  have  termed  the  Northern  Recirculation 
Gyre  (NRG),  is  shown  In  Figure  1.1.  One  of  the  experiments  that  yielded  the 
new  current  meter  data  was  the  Abyssal  Circulation  Experiment  (ABCE).  In 


1983-84  ABCE  was  carried  out  to  learn  more  about  the  NRG.  It  consisted  in 
part  of  a  moored  array  of  current  meters  centered  near  67°W,  40UN.  In  addi¬ 
tion  to  the  evidence  for  tighter  recirculation  it  was  found  that  the  westward 
return  flow  of  the  gyre  is  nearly  depth  independent. 

At  this  point  it  seems  evident  that  there  is  a  substantial  cyclonic  re¬ 
circulation  of  the  deep  Gulf  Stream.  It  Is  not  evident,  however,  what  drives 
this  gyre.  There  has  been  a  considerable  amount  of  modelling  work  that  has 
addressed  in  one  form  or  another  the  presence  of  deep  flow  in  this  region.  A 
feature  of  eddy  resolving  Gulf  Stream  numerical  simulations  is  the  existence 
of  one  or  more  regions  of  closed  circulation  in  the  deep  layer.  Some  models 
exhibit  only  anticyclonic  deep  flow  while  others  contain  both  anticyclonic 
and  cyclonic  deep  gyres  (Holland  and  Lin,  1975).  Harrison  (1982)  proposes 
that  the  cyclonic  recirculation  in  these  models  may  be  inherently  related  to 
quasigeostrophic,  adiabatic  considerations. 

One  such  numerical  experiment  which  contains  deep  cyclonic  flow  is  that 
of  Holland  (1978).  The  model  is  a  two  layer  quasigeostrophic  flow  driven  by 
symmetric  wind  stress.  Holland  and  Rhines  (1980)  analyzed  this  model  in  some 
detail,  and  showed  that  the  deep  gyres  (both  anticyclonic  and  cyclonic)  are 
driven  by  eddy  thickness  fluxes  (or  equivalently,  heat  fluxes)  in  the  surface 
layer.  Hogg  (personal  communication)  has  applied  this  idea  to  the  NRG  using 
a  limited  heat  flux  data  set.  This  first  attempt  has  suggested  that  this  may 
not  be  the  driving  mechanism  In  the  ocean. 

The  mean,  lower  layer  potential  vorticity  field  (Q)  in  the  Holland 
(1978)  simulation  shows  an  area  of  uniform  Q  In  the  region  of  the  two  most 
intense  counter-rotating  gyres.  Rhines  and  Young  (1982a)  have  shown  that  in 
the  presence  of  weak  eddies,  homogenization  of  Q  will  tend  to  occur  within 
closed  streamlines.  Hogg  and  Stommel  (1985)  used  as  a  premise  the  condition 
of  uniform  potential  vorticity  and  derived  deep  cyclonic  recirculation  to  the 
north  of  the  Gulf  Stream  In  an  analytical  framework.  Their  model  relies  on 
the  presence  of  variable  bottom  topography  and  the  thermocline  topography 
associated  with  the  Gulf  Stream  (both  of  which  were  modelled  realistically) . 
In  the  model  they  show  that  the  southward  extent  of  the  gyre  flow  is  deter¬ 
mined  by  the  position  of  the  surface  Gulf  Stream,  which  lends  support  to  the 
idea  of  tight  recirculation  suggested  by  the  newer  data. 
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One  of  the  inconsistencies  between  the  numerical  models  and  the  data  is 
that  in  the  ocean,  directly  beneath  the  surface  Gulf  Stream,  is  the  westward 
return  flow  of  the  NRG.  In  the  models  the  separated  Gulf  Stream  extends  from 
top  to  bottom.  Richardson  (1985)  constructed  from  a  combination  of  current 
meter,  surface  drifter,  and  SOFAR  float  data  an  average  velocity  section  at 
55 °W  which  seems  to  resolve  this  inconsistency.  He  shows  a  Gulf  Stream  that 
is  top  to  bottom  but  in  the  vertical  is  sloped  to  the  south.  On  either  side 
of  the  Stream  is  westward  flow,  consistent  with  the  double  gyre  scheme  of 
Hogg  (1983).  The  other  inconsistency  between  the  model  ocean  and  real  ocean, 
that  regarding  the  driving  mechanism,  remains  to  be  sorted  out.  Hogg  and 
Stomrnel  (1985)  have  revealed  some  elements  that  seem  important  in  regard  to 
cyclonic  recirculation,  but  this  has  no  direct  bearing  on  what  forces  the 
flow. 

The  homogenization  that  accompanies  flow  within  closed  streamlines 
reveals  itself  in  another  context  as  well,  that  of  passive  tracers.  The 
uniformity  of  Q  is  crucially  tied  to  the  structure  and  dynamics  of  the 
circulation  field,  which  also  makes  it  difficult  to  address.  In  terms  of  a 
passive  tracer  however  homogenization  is  more  easily  studied,  which  suggests 
that  the  case  of  a  passive  tracer  be  carefully  examined.  It  is  the  hope 
that  the  ideas  and  insights  developed  in  these  simpler  surroundings  can  then 
be  applied  to  the  more  complicated  case  of  a  dynamically  active  quantity. 

Musgrave  ( 1 985)  did  a  numerical  study  of  the  homogenization  of  passive 
tracers  in  the  thermocline  of  a  subtropical  gyre.  The  process  he  modelled 
was  that  of  tracer  being  subducted  into  the  gyre  from  outcropping  lines  in 
the  northern  regions.  The  abyssal  gyres  are  not  subject  to  this  type  of 
ventilation.  However,  as  depicted  in  Figure  1.1,  a  portion  of  the  Northern 
Recirculation  Gyre  passes  very  closely  to  the  North  Atlantic  Deep  Western 
Boundary  Current  (DWBC)  which  flows  along  the  continental  slope.  The  water 
in  the  boundary  current,  having  recently  come  from  high  latitudes,  has  very 
distinctive  characteristics,  and  so  the  DWBC  represents  a  source  of  tracer 


into  the  deep  layer. 

The  other  component  of  ABCE  was  a  hydrographic  cruise  covering  a  siz¬ 
able  portion  of  the  NRG,  in  which  water  sample  data  was  also  collected.  Hogg 
et  al .  (1986),  using  data  from  several  previous  cruises  as  well,  mapped  out 


the  corresponding  lateral  distributions  of  various  tracers.  The  conspicuous 
feature  in  all  of  the  distributions  was  a  region  of  very  weak  gradients, 
believed  to  be  in  response  to  the  stirring  action  of  the  cyclonic  recircu¬ 
lation.  Hogg  et  al.  presented  a  simple  model  of  the  manner  in  which  tracer 
diffuses  from  the  boundary  current  and  subsequently  becomes  homogenized  with¬ 
in  the  gyre.  In  this  work  we  examine  this  process  in  some  detail.  We  anal¬ 
yze  more  thoroughly  various  aspects  of  the  model  used  by  Hogg  et  al . ,  and 
apply  the  results  to  the  ABCE  data  set. 

The  DWBC  is  characterized  by  a  distinct  signal  in  salt,  oxygen,  silica, 
freon,  and  various  other  tracers.  As  described  in  Hogg  et  al .  (1986)  it  is 
believed  that  the  eastward  flow  of  the  NRG  pulls  a  plume  of  these  tracers 
away  from  the  boundary.  This  type  of  process  was  mentioned  by  McCartney 
et  al.  (1980)  in  describing  a  meridional  section  of  silica  at  55°W.  Here  we 
consider  only  the  distributions  of  salt,  oxygen,  and  freon.  By  studying  the 
interaction  between  diffusion  from  the  current  and  advection  from  the  nearby 
recirculation  we  can  obtain  information  about  the  entrainment  and  homogeniza¬ 
tion  that  occurs  in  a  gyre,  and  about  specific  characteristics  of  the  NRG. 

To  aid  in  this  study  a  simple  two-dimensional  advective-diffusive  model 
was  implemented  which  was  designed  to  represent  the  UWBC/NRG  system.  The 
velocity  field  in  the  model  is  specified,  and  a  source  of  tracer  is  intro¬ 
duced.  In  chapter  one  we  discuss  some  of  the  results  of  the  numerical  model. 
We  split  the  analysis  into  two  parts,  the  first  part  focussing  on  the  process 
oy  which  tracer  penetrates  the  streamlines  of  the  gyre  during  spin-up.  A 
single  numerical  experiment  is  analyzed  to  understand  the  details  of  how  this 
occurs.  The  entrainment  is  characterized  by  a  plume  of  tracer  which  extends 
from  the  boundary  current  and  wraps  into  the  gyre,  spiralling  across  stream¬ 
lines  towards  the  gyre  center.  We  first  examine  what  factors  cause  the  plume 
to  cross  streamlines  in  the  simpler  context  of  a  rectilinear  shear  flow. 
Results  obtained  for  this  idealized  flow  pattern  are  then  applied  to  the 
full-blown  case  of  the  gyre  in  the  numerical  model. 

The  second  part  of  the  numerical  model  analysis  focusses  on  the  occur¬ 
rence  of  homogenization  within  the  gyre  at  steady  state.  We  discuss  how  the 
velocity  structure  of  the  gyre  and  the  character  of  the  spiralling  plume  are 
tied  into  this  process.  Results  from  several  different  numerical  runs  are 


compared  as  well  to  illustrate  the  effect  that  varying  the  diffusivity  has 
on  homogenization. 

Chapter  one  deals  with  the  model  gyre  circulation  alone,  and  results 
obtained  apply  to  any  closed  circulation  of  this  form,  not  just  the  NRG.  In 
chapter  two  we  take  a  broader  perspective  and  include  the  fact  that  the 
tracer  filling  the  interior  diffused  from  the  boundary  current.  We  use  the 
numerical  model  to  help  analyze  tracer  distributions  from  ABCE  in  order  to 
learn  more  about  the  NRG.  To  make  the  comparison  between  model  and  data  more 
revealing  a  simplified  box  model  representation  of  the  numerical  model  is 
solved  analytically.  We  first  apply  the  box  model  to  the  numerical  model  in 
order  to  help  interpret  some  of  the  results.  Tne  information  obtained  by 
this  comparison  is  then  used  to  interpret  the  data.  Among  other  things  we 
learn  about  the  size  of  the  NRG,  how  much  tracer  it  entrains,  and  what  values 
of  diffusivity  (lateral  and  vertical)  are  associated  with  the  flow.  We  also 
learn  what  factors  dictate  how  quickly  the  gyre  is  filled  and  relate  this  to 
the  analytical  results  of  Rhines  and  Young  (1983). 

In  their  work  regarding  the  tracer  distributions  associated  with  the 
NRG,  Hogg  et  al .  (1986)  discussed  a  subtle  difference  between  those  distribu¬ 
tions  of  salt  and  oxygen.  They  explained  the  difference  as  a  result  of  ver¬ 
tical  mixing.  In  chapter  two  this  distinction  is  considered  in  greater 
detail.  We  are  able  to  distinguish  between  a  salt- type  tracer  and  oxygen- 
type  tracer  in  the  numerical  model  and  box  model.  (Chapter  one  deals 
exclusively  with  a  salt- type  tracer,  so  the  results  apply  directly  to  the 
ABCE  salinity  data.)  In  chapter  two  we  compare  the  results  of  two  numerical 
runs  that  include  vertical  mixing  with  two  of  the  runs  of  chapter  one  (with¬ 
out  vertical  mixing),  and  show  that  some  of  the  information  about  the  NRG  can 
only  come  from  an  oxygen- type  tracer.  A  further  run  without  vertical  mixing 
is  presented  to  illustrate  the  effect  of  western  intensification  of  the  gyre. 

Freon  is  also  an  oxygen- type  tracer,  but  it  is  unique  in  that  the 
amount  of  freon  in  the  world  oceans  is  increasing  very  rapidly.  While  it  is 
believed  that  the  oxygen  and  salt  distributions  in  the  region  of  the  NRG  are 
close  to  steady  state,  the  distribution  of  freon  is  constantly  changing.  In 
order  to  understand  the  evolution  it  is  necessary  to  know  how  the  DWBC  source 
strength  changes  in  time.  It  is  evident  then  that  to  study  freon,  a  regional 
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model  is  not  sufficient;  we  must  consider  the  evolution  of  freon  in  high 
latitude  source  waters,  and  throughout  the  DWBC. 

Chapter  three  addresses  this  previous  history  of  the  freon  and  now  it 
affects  the  distribution  found  in  the  NRG.  The  treatment  is  in  three  parts: 
the  overflow  process  in  the  Norwegian-Greenland  Sea  which  leads  to  the  forma¬ 
tion  of  the  DWBC,  the  advection  and  mixing  that  occurs  in  the  boundary  cur¬ 
rent,  and  finally  the  diffusion  into  the  NRG  (within  the  regional  domain  of 
the  previous  chapters).  The  treatment  in  chapter  three  is  analytical. 

Simple  models  are  used  to  determine  freon  concentrations  in  the  overflow 
water  and  DWBC,  and  this  information  is  used  to  drive  the  box  model  of  the 
NRG  region.  (It  is  not  feasible  to  use  the  numerical  model  of  the  previous 
chapters.)  Freon  has  been  in  the  ocean  since  roughly  1950.  Its  atmospheric 
source  function  is  well  defined  and  in  recent  years  measurement  techniques 
have  improved,  both  of  which  make  freon  an  appealing  tracer  to  study.  Cor¬ 
respondingly  freon  studies  have  become  more  common,  and  experiments  have  now 
been  carried  out  in  various  regions  including  the  Arctic,  Antarctic,  and 
Mediterranean.  Studies  in  the  Atlantic  include  TTO  (Transient  Tracers  in 
the  Ocean)  and  A BCE. 

Using  the  atmospheric  freon  concentrations  and  seawater  solubilities, 
estimates  have  been  made  of  the  rate  in  which  water  is  being  transported 
from  northern  latitudes  along  the  western  boundary  of  the  Atlantic  (Smethie 
and  Trumbore,  1984;  Weiss  et  al . ,  1985).  When  this  type  of  calculation  is 
applied  to  the  DWBC,  it  predicts  a  core  speed  of  ~1  cm/sec.  These  estimates 
do  not  take  into  account  water  formation  processes  and  also  rely  on  the 
assumption  that  the  dilution  which  occurs  is  with  freon-free  water.  In  chap¬ 
ter  three  we  address  both  of  these  points.  We  find  that  the  overflow  process 
has  significant  bearing  on  the  calculation,  and  using  two  separate  boundary 
current  models  arrive  at  larger  core  speeds  for  the  DWBC.  Estimates  of  the 
diffusivities  that  come  out  of  these  models  are  similar  to  the  independent 
estimate  from  the  NRG  calculation  of  chapter  two.  We  find  also  that  freon 
is  only  now  beginning  to  accumulate  within  the  Northern  Recirculation  Gyre. 


CHAPTER  ONE:  GYRE  PROCESSES 


Preliminaries 

The  northern  recirculation  gyre  (NRG)  is  believed  to  be  a  weakly  depth 
dependent  flow  transporting  roughly  20  Sverdrups,  driven  by  eddies  from  the 
surface  Gulf  Stream  (Hogg  et  al.,  1986).  During  a  portion  of  its  circuit  the 
water  passes  closely  to  the  DWBC  (Figure  1.1)  at  which  point  it  is  ventilated 
by  lateral  diffusive  transfer  of  various  water  properties  from  the  current 
(ventilated  in  the  sense  that  the  gyre  is  replenished  by  younger  boundary 
current  water).  A  simple  two-dimensional  numerical  model  was  constructed  to 
represent  this  process.  The  streamlines  of  the  model  are  shown  in  Figure  1.2 
and  consist  of  a  boundary  current  situated  alongside  a  gyre.  Flow  speeds  are 
representative  of  the  current  meter  data  in  Hogg  (1983),  and  the  size  of  the 
gyre  is  roughly  that  suggested  by  Hogg.  At  the  northern  edge  of  the  model 
boundary  current  a  steady  Gaussian  source  of  tracer  is  specified,  and  at  the 
southern  edge  tracer  is  allowed  to  advect  out  of  the  domain.  Everywhere  else 
along  the  boundary  there  is  no  flow  (velocities  there  were  set  identically 
equal  to  zero)  and  open  boundary  conditions  enable  tracer  to  diffuse  out  of 
the  region. 

The  evolution  of  tracer  in  the  interior  is  governed  by  a  finite- 
difference  approximation  of  the  two-dimensional  advective-diffusive  equation. 


~  +  u  '  ve  =  v  •  kV9 
at 


(1.1) 


where 


e(x,y)  =  tracer  concentration. 


u(x,y)  = 


eddy  diffusivity  (constant), 
velocity  vector. 


and  !  -  I  lx  +  i  ly  • 

The  circulation  is  steady,  and  initially  the  domain  is  tracer  free.  We  set 
the  diffusivity  K  =  10  cm  /sec.  Details  of  the  set  up  and  numerics  of  the 
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model  are  discussed  in  Appendix  A.  The  scheme  that  was  used  is  based  on  that 
of  Smolarkiewicz  (1983),  with  the  addition  of  a  centered- differenced  diffu¬ 
sion  term,  and  so  includes  a  corrective  step  to  minimize  implicit  diffusion. 
As  the  simulation  is  allowed  to  progress,  tracer  advects  downstream  and 
spreads  laterally.  Some  of  it,  having  diffused  into  the  edge  of  the  gyre,  is 
pulled  eastward  forming  a  plume  which  wraps  around  the  gyre  (Figure  1.3). 
Tracer  slowly  fills  the  gyre  in  this  manner  until  at  steady  state  a  homogene¬ 
ous  pool  forms  within  the  gyre. 

The  focus  of  this  chapter  is  an  analysis  of  the  processes  involved  in 
the  subsequent  entrainment,  i.e.  once  tracer  has  entered  the  edge  of  the 
gyre.  The  ideas  that  are  discussed  apply  then  not  only  to  the  NRG  but  to 
any  closed  circulation  of  this  form  near  an  external  source.  The  discussion 
first  centers  on  the  initial  penetration  of  tracer  into  the  gyre.  This  pro¬ 
cess  is  isolated  in  the  context  of  a  simpler  flow  field.  Then  the  occurrence 
of  homogenization  at  steady  state  is  discussed. 


Initial  Penetration 


Closer  inspection  of  Figure  1.3  shows  that  the  plume  of  tracer,  as  it 
winds  around  the  gyre,  migrates  across  streamlines  toward  the  gyre  center. 

The  reason  for  this  spiral  is  that  the  portion  of  the  plume  which  spreads 
inward  enters  a  region  of  stronger  velocity  and  advects  around  more  quickly. 
Note  also  that  the  spiral  is  asymmetric  in  that  where  the  flow  is  zonal  the 
spiral  is  not  as  pronounced  as  in  the  meridional  flow.  The  spiral  character¬ 
izes  the  entrainment  of  tracer  into  the  gyre  and  it  is  of  interest  to  con¬ 
sider  it  in  some  detail.  To  understand  why  the  asymmetry  exists  it  must  be 
understood  what  factors  govern  the  spiral.  To  do  this,  a  problem  involving 
diffusion  in  a  simple  shear  flow  is  considered. 

The  effect  that  velocity  shear  has  on  the  spreading  of  a  passive  tracer 
has  been  studied  considerably,  in  particular  the  process  of  shear  dispersion 
whereby  cross-stream  shear  enhances  the  spreading  of  tracer  along  streamlines 
(Rhines,  1983).  Here  a  different  aspect  in  which  cross-stream  shear  influ¬ 
ences  the  diffusion  of  tracer  is  addressed. 

For  a  given  distribution  of  tracer  consider  the  parameter  which  is  the 
ratio  of  the  alongstream  gradient  to  the  cross-stream  gradient: 
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Figure  1.3:  Plume  of  tracer  spiralling  into  gyre  (instantaneous  distribution 
during  spin-up).  The  dark  lines  are  the  bounding  streamlines  of  the  gyre  and 
boundary  current. 


The  value  of  is  one  measure  of  the  extent  to  which  shear  dispersion 
occurs.  For  the  same  shear  and  diffusivity,  a  large  «L  means  prevalent 
shear  dispersion  whereas  a  small  means  this  effect  is  negligible.  Shear 
dispersion  acts  on  a  distribution  of  tracer  to  reduce  its  (Rhines,  1983). 

Here  we  are  interested  in  the  effect  that  cross-stream  shear  has  in  the 
spreading  of  tracer  across  streamlines,  when  the  distribution  of  tracer  is 
characterized  oy  a  small  a^.  The  analysis  applies  to  situations  in  which 
there  is  a  localized  source  of  tracer.  Such  a  distribution  in  a  linear  shear 
flow  is  analogous  to  tne  plume  of  tracer  penetrating  the  edge  of  the  gyre. 

A)  Linear  Shear  Flow 

The  equation  governing  the  horizontal  evolution  of  tracer  is  the  two- 
dimensional  advective-diffusive  equation, 

et  *  “x  *  v«y  *  "V  *  K0yy  •  a-2» 

where  x  =  zonal  distance, 

y  *  meridional  distance, 
e  a  concentration  of  tracer, 
u  a  zonal  velocity, 
v  a  meridional  velocity, 
and  <  a  diffusivity. 

Solutions  were  obtained  numerically  using  the  finite  grid  approximation  dis¬ 
cussed  in  Appendix  A.  In  the  region  of  inflow  a  Gaussian  concentration  of 
tracer  is  assigned,  and  where  there  Is  outflow  tracer  advects  out  of  the 
domain.  Where  the  cross  boundary  flow  is  insignificant  the  diffusive  open 
boundary  conditions  are  applied. 

The  set  up  of  the  problem  is  depicted  in  Figure  1.4.  At  t  =  0  a 

step  function  source  is  turned  on,  and  tracer  progresses  downstream  while 

spreading  laterally.  The  center  of  mass  of  the  tongue  proceeds  to  migrate 
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Figure  1.4:  Graphic  depicting  a  plume  of  tracer  In  a  linear  show  flow. 
Length  scales  of  the  plume  are  as  shown;  tracer  Is  Introduced  at  u  =  Uf 


across  streamlines  as  with  the  gyre  flow.  It  Is  relevant  to  define  four 
length  scales  for  this  problem: 


L  =  The  cross-stream  penetration  of  tracer,  defined  as  the  maximum 
meridional  extent  of  a  chosen  concentration  of  tracer. 

L  =  The  alongstream  length  of  the  tongue,  defined  as  the  zonal  dls- 

d 

tance  to  where  the  meridional  penetration  Is  greatest. 

L  *  The  displacement  of  the  center  of  mass  of  the  tongue  across 
cm 

streamlines  at  the  point  where  the  meridional  penetration  is 

greatest.  This  measures  migration  of  the  tongue.  (Note  that 

migration  requires  cross-stream  shear.) 

L  =  L  -  L  .  This  measures  spreading  of  the  tongue, 
s  p  cm 


Tne  quantities  L  and  L  are  the  respective  x  and  y  length  scales  of 
a  p 

the  tracer  distribution,  Lcm  and  L$  are  the  first  and  second  y-moments 
(Figure  1.4). 

For  a  northward  diffusing  particle  of  tracer,  consider  the  balance 
between  advectlon  and  diffusion  where  u  *  UQ  ♦  <*y  (v  *  0), 

(U+ay)e3<e+Ke,  ( 

'  o  ^ '  x  xx  yy 


where  UQ  =  (constant)  reference  velocity, 

and  a  =  cross-stream  shear. 


We  estimate  the  order  of  magnitude  of  each  term  In  (1.3)  using  the  x 
and  y  length  scales,  and  define  the  following  non-dimensional  parameters: 

L 

6  =  -j-^-  Is  tne  aspect  ratio, 
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o  D  La  .  ..  _  n  alongstream  diffusive  time  scale 

- '  ■  1s  the  alongstream  Pec  let  number  =  - *  advectlve  time  scale - 


2  „  ,  .  cross-stream  diffusive  time  scale 

P  =  P  «  Is  the  cross-stream  Peclet  number - adve'cfl ve ~ffme~s'cTre - 


In  terms  of  these  parameters  (1.3)  becomes 


2  _ 


62  ♦  1 


(1.4) 


(1)  Large  alongstream  Peclet  number:  relationships  between  length  scales. 

Consider  first  the  limit  of  small  diffusivlty  where  Pa  >>  1  but  Pc 
remains  0(1).  In  this  limit  the  aspect  ratio  will  necessarily  be  small, 

6  <<  1,  and  the  dominant  balance  in  (1.4)  Is 

P  62  ~  1  .  (1.5) 
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Note  that  the  alongstream  Peclet  number  is  composed  of  two  parts  which  can  be 
thought  of  as  two  separate  alongstream  Peclet  numbers,  one  for  the  shear  part 
of  the  flow  and  one  for  the  uniform  part.  We  define  the  parameter  S  as  the 
ratio  of  these  two  Peclet  numbers,  which  Is  a  measure  of  the  shear  which  the 
tracer  experiences. 


Using  this,  (1.5)  can  be  rewritten. 
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Limits:  S  <<  1 

This  condition  causes  the  shear  to  be  negligible,  and  Lp  obeys  the 

rule. 


*L  1/2 
(tt-1) 


(1.6) 


S  >>  1 

In  this  case  the  shear  is  so  strong  that  the  reference  velocity  is 
negligible.  Here  L  obeys  the  rule. 


Both  the  shear  and  the  reference  velocity  are  important. 


1, 


L 


P 


Three  different  examples  of  distributions  in  which  P  >>  1,  P  -  1 

d  C 

appear  in  Figure  1.5.  Each  one  represents  a  snapshot  as  the  tongue  evolves. 
In  the  first,  corresponding  to  S  <<  1, 


S  =  . 1,  and 


Lp  Ls  Lcm 


In  the  second,  corresponding  to  S  >>  1, 

S  .  2.3.  and  Lp  -  Lc„  »  L$  . 

In  the  third. 


S  =  .6,  and  L  >  L  -  L  . 

p  cm  s 

(A  complete  listing  of  parameters  appears  in  Table  1.1.) 

As  seen  in  Figure  1.5,  as  the  plume  in  the  first  example  progresses 
downstream  Lp  Is  consistent  with  the  S  <<  1  law,  and  in  the  second 
example  it  is  consistent  with  the  S  >>  1  law.  (In  the  former,  where  the 
shear  is  negligible,  an  analytic  solution  is  obtainable  which  agrees  with 
tne  numerical  result  that  Lp  traces  out  a  parabola.)  In  the  third  example 
Lp  corresponds  to  neitner  of  these  laws.  Initially  the  slope  is  close  to 
\U  as  the  shear  is  not  yet  felt,  but  approaches  that  of  1/3  as  penetration 
Increases. 
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TABLE  1.1:  Parameters  associated  with  the  different  examples  In  the  linear 
shear  flow  problem. 


S  P.  a  U  (cm/sec)  L_  (km  x  100) 
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When  the  shear  is  negligible,  spreading  of  the  tongue  accounts  for  most 
of  the  penetration  of  tracer  across  streamlines.  On  the  other  hand,  in  a 
strongly  sheared  flow  the  penetration  is  mostly  due  to  migration  of  the 
tongue.  In  the  third  example,  which  is  between  these  extremes,  spreading  and 
migration  are  both  substantial,  but  with  increasing  penetration  the  situation 
approaches  that  of  the  shear  extreme,  and  correspondingly  Lcm  becomes  more 
closely  correlated  with  Lp.  Note  that  this  example  does  not  correspond 
exactly  to  the  S  ~  1  case  discussed  above.  This  means  that  when  u  and  au 
are  of  equal  magnitude,  spreading  and  migration  do  not  contribute  equally  to 
tne  penetration  but  rather  migration  is  somewhat  more  prevalent. 

These  results  can  also  be  obtained  analytically  through  an  analysis  of 
a  slightly  different  problem,  that  of  a  point  discharge  of  dye  in  a  linear 
shear  flow.  Smith  (1982)  solved  this  case,  and  while  the  discharge  is  not  a 
continuous  source  but  rather  an  initial  spot  that  evolves,  the  same  informa¬ 
tion  regarding  penetration  can  be  derived.  In  terms  of  the  present  notation 
the  flow  field  considered  by  Smith  is 

u  =  UQ(  t)  +  a(  t)y  , 

and  the  level  of  discharge  is  y  =  0.  In  our  case  both  the  reference 
velocity  and  cross-stream  shear  are  constants,  UQ( t )  =  UQ  and  a(t)  =  a. 

The  solution  is  a  Gaussian  in  the  cross-stream  direction  of  the  form 

e(x,y,t)  =  c(x,t)  exp 

where  v(t)  =  variance, 

y0(x,t)  =  position  of  the  origin, 

and  c(x,t)  =  amplitude. 

In  each  section  across-stream  the  Gaussian  is  centered  progressively  further 
to  the  north  looking  downstream,  with  a  different  amplitude.  The  variance 
is  independent  of  alongstream  direction  but  varies  in  time.  A  snapshot  of  a 
spot  of  tracer  progressing  downstream  appears  in  Figure  1.7. 

The  analogy  to  the  continuous  source  case  is  that  the  leading  edge  of 
tne  dye  spot  evolves  the  same  as  the  leading  edge  of  the  plume  discussed 
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Figure  1.7:  Instantaneous  distribution  of  tracer  in  the  linear  shear  flow  of 
Figure  1.5b.  Initially  the  distribution  was  a  delta  function  at  the  origin. 


previously.  The  penetration  here  is  defined  exactly  as  before:  the  maximum 
cross-stream  extent  of  a  given  concentration  isoline.  Note  that  for  the 
example  shown  this  does  not  occur  at  the  extreme  downstream  edge  of  the  dye 
spot,  but  rather  just  upstream  of  this  (Figure  1.8).  As  before  we  define  the 
alongstream  position  that  corresponds  to  the  maximum  cross-stream  penetration 


as  x  =  L 


In  the  limit  of  small  shear 


L  corresponds  to  the  center  of 

a 


the  dye  spot,  whereas  in  the  opposite  limit  La  corresponds  to  its  down- 

a 


stream  edge. 

The  penetration  is  equal  to  the  sum  of  the  distance  to  the  origin  of 

v  (yQ  and  /v 


the  Gaussian  and  the  square  root  of  its  variance,  Lp  =  yQ  + 


are  analogous  to  Lcm  and  L$ 


S  is  defined  as  before. 


respectively,  Figure  1.8).  The  shear  parameter 
aL- 

In  Figure  1.9  the  relationship  between 


L  and  L  is  plotted  as  time  progresses  for  a  strongly  sheared  flow, 

pa 


A1  so 


shown  is  S  versus  L  .  Consistent  with  what  was  found  earlier,  when  the  pen¬ 


etration  is  so  small  that  the  shear  is  negligible  with  respect  to  the  refer¬ 


ence  velocity  (S  <<  1)  Lp  obeys  a  square  root  law.  At  longer  times 


when  tiie  dye  spot  has  diffused  far  enough  across  stream  that  the  opposite  is 


true  (S  >>  1),  Lp  obeys  a  cube  root  law.  It  is  evident  that  in  the  weak 
snear  limit  Lp  ~  / v ,  whereas  in  the  strong  shear  limit  Lp  ~  yQ. 


(ii)  Small  alongstream  Peclet  number:  enhancement  of  spreading. 

In  the  first  set  of  examples  it  is  seen  that  for  S  >>  1  migration  of 
the  plume  (i.e.  movement  of  its  center  of  mass)  is  more  prevalent  than 
spreading,  and  for  S  <<  1  the  opposite  is  true.  In  each  of  these  cases 
Pa  >>  1.  With  a  smaller  P  the  system  becomes  less  sensitive  to  the  vel- 

a  a 

ocity  and,  more  importantly,  changes  in  the  velocity.  Tnus  we  might  expect 
that  a  reduced  P  will  diminish  the  importance  of  migration  versus  spread- 

d 

ing  in  contriDuting  to  the  penetration,  as  is  the  case  with  a  reduced  S. 

The  distinction  between  S  <<  1  versus  P  <<  1  should  remain  clear  however: 

d 

in  the  first  instance  the  cross-stream  change  in  velocity  is  not  important 
because  it  is  small,  in  the  second  instance  it  is  not  important  because  the 
system  does  not  recognize  it. 
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The  relationship  between  La  and  LE 


Figure  1.9:  (a)  The  relationship  between  La  and  Lp  for  the  distribution 

in  the  point  discharge  problem.  The  early  stage  of  evolution  is  character¬ 
ized  by  a  slope  of  one-half  and  later  stages  by  a  slope  of  one-third. 
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(b)  The 
stream. 
In  both 


value  of  the  shear  parameter  as  the  distribution  progresses  down- 

The  enlarged  circle  and  square  represent  the  same  stage  of  evolutl 
t a;  and  ( b) . 


A  second  set  of  examples  appears  in  Figure  1.10.  In  1.10a  P.  "  10 

d 

whereas  in  1.10b  P  ~  .1  (S  is  comparable  in  each,  see  Table  1.1).  Indeed 

d 

with  a  large  alongstream  Peclet  number  LC(J)  >  L$,  and  with  a  small  along- 

stream  Peclet  number  Le  >  L„.  Mote  in  1.10b  that  L„  ~  Li(  so  that  for 

this  example  the  balance  of  terms  in  (1.5)  is  no  longer  applicable,  i.e.  the 

aspect  ratio  is  now  0(1)  and  the  alongstream  diffusive  flux  term  must  be 

retained.  Here  both  the  alongstream  and  cross-stream  Peclet  numbers  are 

small,  whereas  previously  P  >>  1  and  P  ~  1.  The  dominant  balance  in 

d  c 

(1.4)  for  this  example  is  thus 


(ill)  Discussion 

It  is  seen  that  variation  In  the  alongstream  Peclet  number  P  alters 

d 

the  Importance  of  the  alongstream  diffusive  term  versus  the  advective  term 
in  balancing  the  cross-stream  diffusion.  Variation  In  the  shear  parameter  S 
on  the  other  hand,  serves  to  enhance  or  diminish  advection  by  a  constant  vel¬ 
ocity  field  versus  a  sheared  velocity  field.  This  means  that  two  criteria 
must  be  satisfied  in  order  to  obtain  migration  of  the  plume.  First,  P, 

d 

must  be  large  enough  so  that  the  system  is  sensitive  to  the  velocity  field. 
This  condition  is  necessary  but  not  sufficient.  In  addition,  S  must  be 
large  enough  so  that  the  cross-stream  shear  is  significant. 

The  parameter  S  depends  on  the  meridional  length  scale  Lp  which  is 
only  known  after  the  tongue  has  evolved.  As  with  all  problems  In  scaling,  a 
discussion  of  relevant  balances  requires  some  knowledge  of  the  solution.  It 
would  be  desirable  however  if,  given  the  values  of  a,  UQ  and  <,  we  could 
say  whether  or  not  the  tongue  will  migrate  (provided  the  diffusivity  is  small 
enough  that  the  advective  flux  is  important).  What  we  can  do  is  cast  the 
answer  in  terms  of  L  . 

d 

The  idea  is  to  substitute  for  i  in  the  expression  for  S.  When 

S  <<  1  we  use  (1.6)  with  the  definition  of  S  to  obtain 
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Figure  1.10:  (a)  Instantaneous  distribution  of  tracer  In  which  Pa  >>  1 
Pc  “  1.  The  value  of  S  Is  such  that  migration  Is  more  prevalent  than 
spreading. 


Provided  UQ  Is  not  Identically  zero  (it  can  be  arbitrarily  small)  we  have 
seen  that  for  every  problem  the  shear  is  at  first  negligible,  but  eventually 
dominant.  We  thus  set  S  =  1  in  (1.8)  and  determine  the  value  of  L_, 

a 

L.  =  L.  ,  when  the  transition  occurs,  i.e  when  (1.8)  breaks  down.  For 

a  a-j- 

those  cases  in  which  the  tongue  migrates,  the  transition  occurs  "instantly" 

2 

(in  the  example  of  Figure  1.5b  L  ~  10  km)  and  for  those  in  which  the 

aT 

tongue  spreads,  the  transition  “never"  occurs  (in  the  example  of  Figure  1.5a 

L  -  105  km). 
aT 

Besides  discussing  trends  in  L  and  its  relationship  with  L  and 
3  p  r  cm 

L$,  it  is  important  to  say  a  few  words  about  the  extent  of  the  penetration. 

As  an  alternative  to  Lp  as  a  measure  of  cross-stream  penetration,  consider 

the  integral  of  tracer  in  the  region  y  >.  0,  x  <_  Lx,  i.e.  the  total  amount 

of  tracer  that  has  penetrated  northward.  Here  lies  a  further  distinction 

between  the  S  »  1  case  and  S  «  1  case.  For  two  plumes  --  one  in  a 

strongly  sheared  flow  and  one  in  a  uniform  flow  --  in  which  L  and  L  are 

■ 

the  same,  there  is  significantly  less  tracer  north  of  y  =  0  in  the  shear 
flow.  This  is  because  the  effect  of  northward  shear  on  a  distribution  of 
tracer  is  to  increase  Its  northward  gradient,  and  (Figure  1.11)  this  causes 
a  southward  flux  of  tracer  across  part  of  the  y  =  0  line.  So  whereas 

tracer  has  penetrated  just  as  far  across  stream  in  the  shear  flow,  there  is 

less  of  it. 

Although  the  P,  <<  1  case  resembles  that  for  Pa  >>  1,  S  «  1  in 

d  a 

that  spreading  of  the  plume  dominates  migration,  these  instances  represent 
opposite  extremes  in  penetration.  For  a  given  flow  field  (UQ  and  a)  and  a 
given  La,  the  value  of  lp  depends  on  the  value  of  tc .  A  sufficiently  small 
<  means  that  Lp  is  not  large  enough  for  the  plume  to  notice  the  shear  -- 
large  P  spreading  occurs.  With  increased  <  (and  L  )  the  large  S 

a  P 

regime  is  approached  and  migration  becomes  important.  Small  Pfl  spreading, 
on  the  other  hand,  occurs  with  large  enough  <,  and  this  represents  the 
upper  extreme  of  penetration. 

8)  Application  to  Gyre  Flow 

We  return  now  to  the  gyre  problem.  The  process  in  which  the  plume  of 
tracer  gets  caught  in  the  edge  of  the  gyre  and  diffuses  into  a  region  of 
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Figure  1.11: 


Value  of  ae  along  the  line  y  =  u- 

(a)  Uniform  flow.  The  meridional  flux  of  tracer 

northward. 


is  everywhere 


stronger  flow  resembles  the  shear  flow  example  analyzed  above,  and  some  of 
the  ideas  that  were  developed  are  now  applied  to  this  problem.  Keep  in  mind 
however  that  there  are  differences  between  the  two  examples.  For  instance 
here  the  cross-stream  shear  a  varies  both  alongstream  and  across  stream, 
there  is  alongstream  shear  as  well,  the  flow  is  curved  rather  than  rectilin¬ 
ear,  and  the  input  of  tracer  into  the  gyre  is  not  a  step  function  in  time. 

Figure  1.12  shows  a  time  history  of  tracer  entering  the  gyre.  The 
advancement  of  the  plume  in  each  one  year  segment  can  be  thought  of  as  dif¬ 
ferent  example  of  the  shear  flow  problem  examined  above,  with  the  following 
definitions: 

L  =  the  distance  that  the  leading  edge  of  the  plume  travels 

a 

alongstream  in  a  year. 

L  =  movement  of  the  center  of  mass  of  the  leading  edge  of  the  plume 
across-stream  in  a  year. 

a  =  cross-stream  shear  at  the  midpoint  of  La  (the  alongstream 
shear  is  negligible). 

UQ  =  velocity  at  =  0,  L  =  0. 

These  quantities  are  analogous  to  those  similarly  named  in  the  previous  shear 
flow  problem.  The  source  at  the  northern  edge  of  the  DWBC  is  applied  as  a 
step  function  at  t  =  0,  but  by  the  time  tracer  reaches  the  gyre  it  is  no 
longer  characterized  by  a  front,  i.e.  the  "source"  for  each  of  the  above  ex¬ 
amples  grows  in  amplitude  and  width.  This  means  we  are  unable  to  define  the 
analog  to  Lp,  which  in  turn  means  we  are  unable  to  measure  directly  the 
values  of  Pa  and  S.  We  can  however  estimate  the  size  of  Pa  by  noting  that 

a  a 

6,  the  aspect  ratio,  is  much  less  than  one  for  each  single  year  segment,  and 
this  necessarily  implies  that  P,  >>  1.  So  the  first  condition  for  migration 

a 

is  satisfied  everywhere  around  the  gyre. 

Because  a  spiral  does  occur  it  is  natural  to  assume  that  the  second 

condition  for  migration,  S  >>  1,  is  satisfied  as  well.  Recall  that  in  the 

migration  limit  L  ~  L  .  Since  we  can  measure  we  are  able  then  to 

3  p  cm  cm 

check  this  assertion.  For  each  single  year  segment  we  can  substitute  the 
values  of  a,  La,  and  <  into  (1.7),  where  Lp  is  replaced  by  Lcm  (the 


Figure  1.12:  Time  sequence  of  tracer  diffusing  from  the  boundary  current  and 
becoming  entrained  into  the  gyre.  Tne  dark  lines  are  the  bounding  stream¬ 
lines  of  the  two  components  of  flow. 
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proportionality  constant  for  (1.7)  was  determined  numerically).  This  pre¬ 
dicted  value  of  Lcm  can  in  turn  be  compared  to  the  measured  value.  The 
results  of  this  comparison  are  displayed  in  Figure  1.13.  It  is  seen  that 
there  is  good  agreement  between  the  predicted  and  measured  values  where  the 
flow  is  meridional.  However,  while  the  predicted  curve  does  display  less 
penetration  in  the  zonal  flow,  the  degree  to  which  it  is  less  is  not  nearly 
as  pronounced  as  that  which  actually  occurs,  suggesting  that  something  else 
is  happening  while  the  plume  travels  zonal ly. 

Figure  1.14  graphs  the  values  of  UQ  and  a  that  the  plume  experiences 
as  it  travels  around  the  gyre.  Also  shown  is  the  extent  of  the  corresponding 
spiral.  The  fact  that  both  a  and  UQ  are  larger  in  the  zonal  flow,  together 
with  the  small  extent  of  migration  there  suggests  that  this  corresponds  to 
the  advective  limit.  In  particular,  in  these  regions  not  only  is  Pa  >>  1, 
but  Pc  >>  1  as  well  so  that  isolines  of  tracer  nearly  coincide  with  stream¬ 
lines.  It  is  more  accurate  then  to  think  of  the  plume  as  mirroring  stream¬ 
lines  when  it  travels  in  the  stronger  zonal  flow,  while  spiraling  across 
streamlines  in  the  manner  of  the  shear  flow  example  when  it  travels  in  the 
weaker  meridional  flow. 

C)  Discussion 

The  preceding  analysis  focussed  on  the  entrainment  of  tracer  into  a 
gyre  which  initially  was  tracer-free.  This  process  is  characterized  by  a 
plume  of  tracer  spiralling  asymmetrically  inward  across  streamlines  as  a  re¬ 
sult  of  the  cross-stream  snear.  In  particular,  in  the  zonal  flow  the  spiral 
is  minimal  as  the  strong  flow  causes  the  plume  to  follow  streamlines.  In  the 
meridional  flow  the  spiral  is  of  considerable  extent  conforming  to  the  ideas 
developed  in  a  simpler  shear  flow  analysis. 

Previous  work  has  been  done  on  the  mixing  of  tracer  within  a  subtrop¬ 
ical  gyre.  Musgrave  (1985)  analyzed  steady  state  solutions  in  which  the 
northern  boundary  is  maintained  at  a  uniform  positive  concentration  while  the 
southern  boundary  is  kept  uniformly  negative.  He  discusses  the  presence  of  a 
spiral  that  extends  from  the  boundary  to  the  stagnation  point  of  the  flow  in 
the  center  of  the  gyre.  The  spiral  arises  because  of  the  choice  of  boundary 
conditions  (the  cross-stream  shear  of  the  gyre  is  of  the  wrong  sense  to  cause 
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the  type  of  spiral  discussed  here,  l.e.  a  is  everywhere  <  Oj.  As  tracer 
enters  from  the  northern  boundary  it  travels  anticyclonical ly  and  spreads 
into  the  interior.  Upon  encountering  the  negative  plume  that  extends  from 
the  south,  the  region  of  positive  concentration  shifts  away  from  the  bound¬ 
ary,  hence  the  spiral. 

For  the  type  of  gyre  considered  here  the  Inward  extent  of  the  spiral 
depends  on  the  flow  parameters  and  diffusivity.  Far  enough  into  the  gyre  the 
velocity  reaches  a  maximum  before  decreasing  to  zero  at  the  gyre  center.  The 
maximum  distance  over  which  the  spiral  can  extend  is  to  the  region  where  the 
cross-stream  shear  vanishes.  As  steaciy  state  is  approached  the  tracer  within 
the  gyre  becomes  homogenized.  This  process  is  discussed  in  the  next  section. 

Homogenization 

A  feature  of  closed  circulation  is  the  tendency  for  properties  to 
become  homogenized  in  steady  state  in  the  interior  of  the  flow.  The  homogen¬ 
ization  of  potential  vorticity  has  been  studied  by  Rhines  and  Young  (1982a) 
and  is  an  ingredient  in  their  theory  of  wind-driven  ocean  circulation  (Rhines 
and  Young,  1982b).  The  occurrence  of  such  homogenization  is  a  regular  fea¬ 
ture  of  numerical  flow  simulations  (e.g.  Holland  and  Rhines,  1980),  and  evi¬ 
dence  for  this  Is  found  in  data  from  the  North  Atlantic  (McDowell  et  al . , 
1982).  In  the  context  of  passive  tracers,  Musgrave  (1985)  conducted  a  numer¬ 
ical  study  of  homogenization  In  a  subtropical  gyre.  The  situation  he  consid¬ 
ered  however  is  quite  different  than  that  presently  being  addressed.  Rhines 
and  Young  (1983)  examined  the  time  history  of  the  process  by  which  gradients 
are  expelled  from  a  closed  circulation,  iesults  of  which  are  applicable  here. 
In  the  ocean.  Mi i 1 er  (1982)  discussed  the  homogeneity  of  the  salinity  field 
In  a  portion  of  the  subtropical  North  Atlantic.  The  tracer  data  presented 
by  Hogg  et  al.  (1986)  for  the  region  of  the  northern  recirculation  gyre  also 
exhibits  homogenization,  a  result  which  this  study  addresses. 

It  is  useful  to  consider  homogenization  in  the  context  of  a  passive 
tracer,  as  results  may  provide  insight  into  the  more  complicated  case  of  a 
dynamically  active  quantity  such  as  potential  vorticity.  Here  the  occurrence 
of  homogenization  is  related  to  the  penetration  process  that  was  analyzed 
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above,  before  proceeding  with  the  present  problem  though  we  first  review  the 
argument  for  homogenization  following  the  presentation  in  Rhlnes  and  Young 
(1982a). 

A)  Review  of  Homogenization 

Consider  the  steady  state  balance  of  advection  and  diffusion  in  a  gyre, 
governed  by  the  steady  form  of  (1.1), 

u  ‘  ve  =  v  *  <ve  . 

Integrating  over  the  area  Dounded  by  a  streamline  and  applying  the  divergence 
theorem, 

!  ue  *  n  ds  =  '  n  ds  ,  (1.9; 

js  -  -  Js  ~ 

where  S  Is  the  boundary  streamline  and  n  Is  the  unit  normal  to  the 
streamline.  Note  that  the  lefthand  side  of  (1.9)  Is  Identically  equal  to 
zero  because  u  and  n  are  perpendicular. 

In  the  limit  of  strong  advection  the  isolines  of  tracer  nearly  coincide 
with  streamlines,  l.e.  e  *  e(  ) .  This  gives 

*  e(  *  )  v  ij/  , 

*  - 

and  since  the  integral  is  around  a  streamline, 
e(  \i/ )  icV  r ds  =  0  . 

*  is  - 

The  quantity  inside  the  integral  is  positive  definite,  which  further  implies 
that 

e(  =  0  ,  i.e.  e  =  constant  , 

Homogenization  Is  thus  obtained  in  a  strongly  advective  system. 


B!  Spatially  Decaying  Gyre 

The  flow  field  that  Musgrave  (1985)  considered  in  his  numerical  study 
of  homogenization  consisted  of  a  Stommel-type  gyre  in  which  the  strongest 
flow  occurs  at  the  edge  of  the  gyre.  He  defined  a  Peclet  number,  P  =  UL/<  , 
using  the  length  scale  of  the  basin  and  the  characteri Stic  velocity  of  the 
gyre,  and  discussed  the  extent  of  the  homogeneous  pool  versus  P  as  well  as 
the  meridional  flux  of  tracer  versus  P.  In  terms  of  (1.1)  the  Peclet  number 
determines  to  what  extent  advection  balances  diffusion,  and  L  should  be 
defined  in  terms  of  the  tracer  distribution.  It  is  unclear  how  to  discuss 
results  in  terms  of  a  Peclet  number  so  defined,  especially  in  regards  to 
homogenization  when  locally  the  length  scale  becomes  infinite.  What  we  do 
here  is  define  P  in  terms  of  the  plume  of  tracer  which  penetrates  the  gyre, 
as  was  done  in  the  previous  section. 

The  gyre  presently  being  considered  has  its  maximum  velocity  relatively 
close  to  the  center,  decaying  from  this  point  to  the  edge.  (A  velocity  sec¬ 
tion  through  the  gyre  is  shown  in  Figure  1.15.)  For  simplicity,  for  the  time 
being  we  consider  a  symmetric  (i.e.  circular)  gyre.  Think  of  the  gyre  as 
being  divided  into  two  regions:  the  outskirts  of  the  gyre  were  the  flow  is 
weak,  and  the  inner  part  where  tne  flow  is  more  intense  (close  to  the  gyre 
center  the  flow  once  again  becomes  weak.) 

As  the  plume  of  tracer  enters  the  outer,  weaker  part  of  the  flow  it 
spirals  across  streamlines  (provided  the  shear  is  strong  enough)  in  the 
manner  discussed  in  the  previous  section.  This  region  is  characterized  by 
P  >>  1,  P  -  1.  Eventually  the  plume  reaches  strong  enough  flow  (we  clar- 
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ify  below  what  is  meant  by  strong  enough)  that  it  is  nearly  pulled  right 
around  a  streamline.  At  this  point  the  spiral  has  "collapsed"  to  a  stream¬ 
line,  and  here  P,  >>  1,  P_  >>  1.  Recall  that  these  latter  conditions  imply 

a  C 

that  we  are  in  the  advective  limit,  which  is  the  necessary  condition  for 
homogenization.  Specifically  then,  the  outer  region  of  the  gyre  is  where  the 
spiral  occurs,  and  the  inner  region,  delimited  by  the  collapsed  spiral,  is 
where  homogenization  occurs. 

For  the  asymmetric  gyre,  we  saw  earlier  that  where  the  plume  first 
enters  the  gyre  heading  east  the  fast  flow  keeps  it  nearly  tracked  to  a 
streamline.  However,  after  the  plume  has  turned  the  corner  to  the  north. 


Figure  1.15:  North-south  velocity  section  through  the  center  of  the  gyre 
Figure  1.2. 
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the  flow  along  that  same  streamline  weakens  and  the  plume  proceeds  to  spiral 
significantly  inward.  In  this  case  the  division  between  the  two  regions  of 
the  gyre  is  not  so  clear  cut.  However,  homogenization  will  not  occur  until 
the  plume  tracks  a  streamline  around  an  entire  circuit.  So  even  though 
portions  of  the  spiral  may  collapse,  it  is  only  where  the  cumulative  spiral 
collapses  that  marks  the  transition  between  the  two  regions. 

We  now  examine  more  clearly  the  condition  that  the  flow  be  strong 
enough  to  keep  the  plume  from  diffusing  appreciably  across  stream  in  the  time 
it  takes  to  recirculate.  The  advective  limit  corresponds  to  Pa  >>  1  and 

a 

Pc  >>  1.  The  more  stringent  of  these  is  Pc  >>  1,  or  in  terms  of  scales, 


In  considering  a  circuit  around  the  gyre  the  relevant  L,  is  L,  =  Lc  = 

a  _  a  S 

perimeter  of  the  streamline,  and  the  relevant  U  is  U  =  v$  =  average  vel¬ 
ocity  around  the  streamline.  We  define  the  homogenization  function  H  as 
the  ratio  of  these  two  quantities,  which  gives 

L  2 

Pc  =  H  (  ip  )  , 

v  ( 

where  H(  <i>  )  =  ^  ^  ■yp  The  function  H,  which  is  the  inverse  of  the 

circulation  time,  can  be  thought  of  as  a  measure  of  the  tendency  for  homogen¬ 
ization  to  occur  based  only  on  flow  characteristics.  A  larger  H  means  a 
greater  likelihood  for  homogenization. 

Consider  again  the  symmetric  gyre,  whose  streamf unction  is  given  by 

rf 

M  *-2 

*  =  *0U  -  e  )  , 

where  =  amplitude,  L  =  e-folding  scale  of  the  gyre.  In  this  case 
v$  =  v$.  From  (1.10), 


(1.10) 
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The  perimeter  L  =  2»r,  which  gives 


H(  *  )  =  -L  (  t  -  »  ) 
i\T  0 


This  says  that  the  greatest  tendency  for  homogenization  occurs  at  the  center 
of  the  gyre,  decreasing  linearly  with  increasing  \p  .  So  the  part  of  the  gyre 
near  the  center,  where  the  flow  becomes  weak  again,  is  included  in  the  advec- 
tive  region  because  the  flow  is  strong  enough  in  the  sense  that  the  circula¬ 
tion  time  is  small  (the  perimeter  of  a  streamline  is  small).  Contrast  this 
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to  a  solid  body  rotation  gyre,  t  =  ^r  ,  where  the  circulation  time  is  con- 

stant  for  each  streamline.  Here  v$  =  2  <pQr  and  H(  ^  ,  and  the  ten¬ 

dency  for  homogenization  is  the  same  everywhere. 

In  the  asymmetric  gyre  the  velocity  along  a  streamline  varies  around 
the  gyre.  The  streamfunction  is 


4>  =  <l»  (1  -  e 
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where  L  and  L  are  the  x  and  y  e- folding  length  scales.  Note  that 
x  y 


u  *  ds 


rsJS 


where  r  is  the  circulation  around  the  streamline  S.  From  (1.11), 
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which  gives 


H  (  4> )  =  ip  Q  -  <M  , 


(1.12) 


L  2  +  L  2 

where  Hq  =  (-^ — ^ ~ )  •  Figure  1.16  graphs  HQ  as  a  function  of  the 
Lx  Ly 

skewness  of  a  gyre.  It  shows  that  for  a  given  value  of  ip ,  homogenization 
is  more  likely  to  occur  in  a  more  skewed  gyre.  Note  that  this  is  true  even 
though  the  velocity  at  the  two  widest  sections  of  the  gyre  approaches  zero  as 
the  gyre's  skewness  increases.  This  is  because  the  portion  of  the  streamline 
in  these  sections  gets  vanishingly  small  as  well. 


C)  Limits  of  Diffusivity 

We  examined  homogenization  within  the  asymmetric  gyre  of  Figure  1.2  for 

various  values  of  the  diffusivity  <.  The  smallest  value  considered  was  <  ~ 
6  2 

10  cm  /sec,  and  the  homogeneous  pool  of  tracer  that  formed  in  steady  state 
is  shown  in  Figure  1.17.  This  is  the  final  state  of  the  same  example  that 
was  previously  analyzed  in  terms  of  the  asymmetric  spiral.  Figure  1.17b 
shows  the  path  of  the  spiral  and  how  it  indeed  closes  in  on  the  region  that 
eventually  becomes  homogenized. 

6  o 

When  the  diffusivity  is  increased  to  <  ~  5  x  10  cm  /sec  for  the  same 
gyre,  this  in  effect  causes  the  flow  to  appear  weaker  to  the  incoming  plume 
of  tracer.  In  particular,  the  zonal  flow  where  the  plume  first  enters  the 
gyre  is  not  fast  enough  anymore  to  be  in  the  advective  limit,  so  a  pronounced 
spiral  occurs  there  as  well  as  in  the  meridional  flow  (the  asymmetry  no  long¬ 
er  exists,  Figure  1.17dj.  Consistent  with  (1.12; ,  the  plume  now  has  to  pene¬ 
trate  further  into  the  gyre  before  it  encounters  flow  strong  enough  to  induce 
homogenization.  Correspondingly ,  the  steady  state  pool  is  reduced  in  extent 
(Figure  1.17c). 

Upon  increasing  <  even  more  (K  ~  lQ^cm'Vsec)  an  abrupt  transition 
occurs  in  the  manner  in  which  tracer  fills  the  gyre.  What  happens  is  that 
the  meridional  flow,  which  is  weaker  than  the  zonal  flow,  is  essentially 
turned  off.  That  is  to  say  the  diffusive  flux  there  is  now  of  the  order  of 
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Figure  1.16:  The  value  of  the  homogenization  function  versus  the  skewness  of 
a  gyre.  The  skewness  n  Is  defined  as  the  ratio  of  the  major  axis  of  the 
ellipse  to  the  minor  axis.  The  area  of  the  gyre  is  the  same  in  each  case  and 
set  equal  to  it. 
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Figure  1.17:  (a)  Steady  state  distribution  of  tracer  for  <  ~  106cm2/sec. 

The  shaded  region  corresponds  to  that  area  of  the  gyre  in  which  the  gradient 
of  tracer  is  _<  .1  concentration  units/km  x  100.  This  is  taken  as  the  cri¬ 
terion  for  homogenization. 
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(b)  The  region  of  homogenization  in  (a)  shown  in  relation  to 
the  Inward  spiral  of  the  plume  of  tracer  during  spin-up.  The  solid  lines 
are  the  bounding  streamlines. 
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(c)  Steady  state  distribution  for  <  ~ 
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(d) 


Homogenized  region  and  spiral  for  <  ~  5  x  10  cm  sec 


the  advective  flux.  Thus  the  plume  enters  the  gyre,  gets  pulled  eastward, 
turns  the  corner  and  stagnates.  By  the  time  tracer  diffuses  northward  and 
gets  caught  in  the  zonal  flow  and  subsequently  advected  westward,  the  west¬ 
ward  diffusing  tracer  from  the  stagnation  point  has  penetrated  the  center  of 
the  gyre.  So  whereas  in  the  previous  two  cases  tracer  was  advected  complete¬ 
ly  around  the  gyre  and  filled  the  center  in  a  bowl-like  fashion,  here  it  is 
advected  to  the  east  and  proceeds  to  fill  the  gyre  from  east  to  west  (Figure 
1.19).  A  small  amount  of  homogenization  does  occur  at  the  center  of  the  gyre 
(Figure  1.18). 

The  final  case  considered  can  be  thought  of  as  the  diffusive  limit 
7  2 

(  <  ~  5  x  10  cm  /sec).  Here  the  presence  of  the  zonal  flow  is  hardly  felt 
as  well  and  the  manner  in  which  the  gyre  is  ventilated  undergoes  yet  another 
change.  As  shown  in  Figure  1.19  tracer  basically  diffuses  from  west  to  east 
across  the  gyre,  with  an  undulation  corresponding  to  the  eastward  and  west¬ 
ward  flows. 

D)  Conclusions 

Homogenization  is  the  steacty  state  manifestation  of  tracer  penetrating 
a  closed  circulation,  provided  the  system  is  strongly  advective.  For  the 
spatially  decaying  gyre  considered  here,  the  occurrence  cf  homogenization  is 
closely  tied  to  the  characteristics  of  the  spiral  of  the  incoming  plume  of 
tracer  that  forms  during  spin  up.  In  particular,  where  the  spiral  collapses 
to  a  streamline  marks  the  outer  extent  of  the  homogeneous  pool  that  eventual¬ 
ly  develops.  As  the  diffusivity  is  increased  the  size  of  this  pool  shrinks. 
This  is  consistent  with  the  idea  that  homogenization  occurs  more  readily 
nearer  the  center  of  the  gyre  for  this  type  of  flow,  based  on  the  shorter 
circulation  times  there. 

Rhines  and  Young  (1983)  have  investigated  how  long  it  takes  homogeniza¬ 
tion  to  occur  in  a  closed  circulation.  They  showed  that  the  process  occurs 
in  two  stages,  the  slower  of  which  is  the  diffusive  time  of  the  gyre.  The 
problem  they  solved  is  simpler  than  the  one  being  analyzed  here.  Their 
initial  state  consisted  of  a  given  distribution  of  tracer  in  the  gyre,  which 
evolved  within  an  insulating  boundary.  In  the  present  problem  in  order  to 
understand  the  time  ni story  of  homogenization,  the  nature  of  the  sources  and 
sinks  of  tracer  must  be  considered.  This  is  done  in  chapter  two. 
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:  Instantaneous  distribution  of  tracer  during  spin-up  illustrat 
in  which  tracer  fills  the  gyre.  Shown  above  is  a  zonal  section 

center  of  the  gyre. 

(a)  <  ~  lQ^cm^/sec. 
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CHAPTER  TWO:  THE  NORTHERN  RECIRCULATION  GTRE 


Prel iminaries 

Up  until  now  we  have  focussed  on  the  gyre  circulation  of  the  numerical 
model,  isolating  the  processes  involved  in  the  penetration  of  tracer  into  the 
gyre.  Having  understood  the  entrainment  and  homogenization  that  occur  we  now 
consider  the  coupled  gyre/boundary  current  system;  i.e.  we  now  examine  how 
the  DWBC  supplies  tracer  to  the  NRG. 

As  mentioned  earlier,  the  DWBC  is  a  source  of  different  tracers  to  the 
abyssal  North  Atlantic.  Having  been  in  recent  contact  with  the  atmosphere 
its  waters  are  rich  in  oxygen,  tritium,  and  freons.  Its  waters  are  also  rel¬ 
atively  cold  and  fresh  and  marked  by  a  distinct  signal  in  silica.  Recently 
Hogg  et  al.  (1986)  presented  results  from  the  OCEANUS  134  hydrographic  survey 
of  the  region  where  the  NRG  arid  DWBC  are  in  close  contact,  highlighting  sev¬ 
eral  of  these  tracers.  (Some  data  from  previous  cruises  was  also  included.) 
They  discussed  the  data  in  relation  to  this  associated  flow  pattern,  remarx- 
i ng  on  the  area  of  nearly  uniform  tracer  concentration  in  the  region  of  the 
NRG.  It  was  also  shown  that  there  is  no  point  in  considering  each  of  the 
tracers  individually  as  they  do  not  all  give  independent  information.  In 
particular,  the  distribution  of  salinity  was  independent  from  those  of  the 
other  tracers,  the  rest  of  which  show  nearly  identical  features.  This  dif¬ 
ference,  the  presence  of  a  slight  minimum  near  the  location  of  the  NRG  for 
all  of  the  tracers  except  salinity,  was  explained  in  terms  of  vertical 
processes. 

Here  we  consider  three  different  tracers:  salinity,  oxygen,  and 
freons.  Using  the  numerical  model  in  association  with  a  simple  analytical 
model,  we  explore  what  these  tracer  distributions  can  tell  us  about  the  NRG 
and  its  relationship  to  the  DWBC  and  about  various  other  properties  of  the 
system.  We  examine  more  closely  the  distinction  between  salt  and  oxygen 
suggested  in  Hogg  et  al .  (1986).  Freons  are  considered  as  well  because  they 
are  in  a  transient  state  and  the  atmospheric  forcing  functions  are  known. 

An  investigation  of  time  dependent  input  is  the  focus  of  chapter  three. 
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We  start  with  an  analysis  of  numerical  model  results  only.  In  order  to 
understand  these  results  more  clearly  a  simpler  representation  of  the  model 
is  solved  analytically.  Some  of  the  ideas  that  are  derived  from  this  analy¬ 
sis  are  then  applied  to  the  same  data  set  discussed  in  Hogg  et  al .  (1986). 


Model  Study 

It  is  assumed  that  a  turbulent  transfer  of  properties  from  the  DWBC  to 
the  NRG  occurs  along  density  surfaces,  the  entire  gyre  being  ventilated  in 
this  fashion.  It  is  often  the  case  that  along-i sopycnal  mixing  in  the  ocean 
is  presumed  to  dominate  cross-i sopycnal  mixing.  However,  it  is  the  presence 
of  cross-i sopycnal  mixing  that  Hogg  et  al .  suggest  may  be  the  reason  for  the 
difference  in  the  distributions  of  salinity  and  oxygen. 

Consider  the  vertical  profiles  of  Figure  2.1  for  a  station  from  the 
OCEANUS  134  data  set.  At  the  density  level  of  the  DWBC  core  (average  depth 
£  3600  m)  the  salinity  distribution  decreases  monotonically  with  depth  (dense 
water  of  Antarctic  influence  keeps  near  bottom  levels  fresh).  However,  the 
oxygen  profile  has  a  relative  maximum  at  this  level.  Thus  in  light  of  the 
associated  gradients,  for  oxygen  there  is  a  cross-i sopycnal  flux  out  of  the 
deep  layer  into  the  water  above  and  the  water  below,  whereas  for  salt  there 
is  a  flux  into  the  layer  from  above  and  out  of  the  layer  into  the  bottom 
water.  This  suggests  that  the  effect  of  cross-i sopycnal  mixing  in  the  deep 
layer  may  be  more  pronounced  for  oxygen  than  for  salinity. 

Consider  the  three-dimensional  form  of  the  advective-diffusive 
governing  equation  (1.1), 
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2.1:  (a)  CTO  profile  of  salinity  versus  depth  below  2000  m  for 

0C 134  station  40  (in  the  gyre). 
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We  are  interested  in  the  evolution  of  tracer  in  the  deep  layer  which  we 
assume  is  oounded  above  and  below  by  density  surfaces  which  are  nearly  flat 
(thus  along-isopycnal  is  synonymous  with  horizontal,  cross-isopycnal  is 
synonymous  with  vertical). 

In  a  finite-difference  sense  we  represent  the  vertical  structure  with 
three  grid  points,  but  only  the  center  one  is  active,  corresponding  to  the 
concentration  in  the  center  of  the  layer.  The  upper  and  lower  points  are 
fixed  boundary  conditions  (representing  reservoir  values  above  and  below  the 
deep  layer.  Figure  2.2)  .  Applying  this  approximation  to  (2.1)  (w  =  0), 
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where  eM(x,y)  =  value  of  tracer  at  the  center  of  the  layer, 
eu  =  value  of  upper  reservoir  (constant), 

eL  =  value  of  lower  reservoir  (constant), 

H  =  layer  thickness  (vertical  grid  spacing  =  -*) , 


or  in  terms  of  the  anomaly  61  =  (eM 
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This  equation  is  quasi-three-dimensional  in  that  it  contains  a  param¬ 
eterization  of  a  vertical  process.  Note  that  this  vertical  flux  term  has 
tne  form  of  a  radioactive  decay  term  (although  for  a  property  such  as  silica, 
which  is  characterized  by  a  relative  minimum  in  the  deep  layer,  it  is  a 
growth  term).  As  discussed  in  Hogg  et  al.  (1986),  for  the  scales  involved, 
this  type  of  decay  in  oxygen  concentration  in  the  abyssal  ocean  overwhelms 
any  consumption  that  may  be  occurring. 

Equation  (2.2)  is  that  which  was  analyzed  by  Hogg  et  al.  in  differenti¬ 
ating  between  salt  and  oxygen.  That  analysis  is  expanded  upon  here.  For  the 

9  +  e, 

salt  case  it  is  assumed  that  - j — ~  ^  eM  (e'  £  0)  so  that  the  vertical 
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Figure  2.2:  Schematic  showing  the  vertical  resolution  of  the  deep  layer  in 
the  numerical  model.  Three  grid  points  are  used  to  represent  a  continuous 
profile  as  such. 
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flux  term  is  not  large  enough  to  affect  the  distribution.  Setting  it  ident¬ 
ically  equal  to  zero  reduces  (2.2)  to  the  equation  which  was  studied  in  chap¬ 
ter  one,  (1.1).  The  previous  results  concerning  entrainment  and  homogeniza¬ 
tion  can  be  thought  of  then  as  applying  to  the  salt  case.  We  consider  this 

e  +  eL 

case  further.  For  oxygen  it  is  assumed  that  - 2 - <  and  the  vertical 

flux  term  is  retained.  Using  the  same  two-dimensional  finite-difference 
scheme  that  was  applied  to  (1.1),  equation  (2.2)  is  then  analyzed  numerically 
and  compared  to  the  salt  case.  (Hereafter  the  prime  is  dropped  and  e 
refers  to  the  deep  anomaly.) 

A)  Without  Vertical  Flux 

We  turn  our  attention  once  again  to  the  same  numerical  simulation  that 
was  analyzed  in  the  previous  chapter  (the  small  diffusivity  limit  only),  but 
here  we  take  a  broader  perspective  and  consider  the  fact  that  the  tracer 
which  collects  in  the  gyre  originally  came  from  the  boundary  current.  As 
tracer  first  advects  downstream  from  the  northern  source  and  spreads  later¬ 
ally,  the  eastward  flux  of  tracer  into  the  interior  is  inhibited  by  the  west¬ 
ward  flow  of  the  gyre  and  accentuated  by  the  eastward  flow  further  to  the 
south.  All  along  the  region  of  contact  tracer  diffuses  from  the  boundary 
current  directly  into  the  gyre.  However,  because  of  the  strong  gyre  flow, 
tracer  does  not  penetrate  substantially  into  the  gyre  before  the  flow  turns 
off  shore.  For  this  reason  the  plume  of  tracer  which  extends  from  the  bound¬ 
ary  does  not  coincide  with  the  region  of  maximum  eastward  flow  of  the  gyre, 
rather  the  plume  is  well  south  of  this  region  (Figure  2.3).  This  is  an 
example  then  of  when  a  tongue  of  tracer  does  not  coincide  with  the  core  of 
the  current  (rather  it  defines  the  current's  edge). 

Each  successive  plume  originates  from  the  boundary  a  bit  further  to 
the  south,  as  it  is  made  up  of  tracer  that  diffused  further  from  within  the 
current  and  so  was  subject  to  stronger  southward  advection.  Throughout  the 
simulation  the  level  of  tracer  across  the  gyre  is  nearly  flat.  (There  is  a 
slight  minimum  in  the  center  of  the  gyre,  Figure  2.4.)  The  entire  system 
eventually  reaches  a  steady  state  as  input  at  the  northern  edge  of  the  cur¬ 
rent  is  balanced  by  advective  output  at  its  southern  edge  and  diffusive  out¬ 
put  along  the  remaining  (quiescent)  part  of  the  boundary.  The  gyre  itself 
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Figure  2.4:  Zonal  section  of  tracer  through  the  boundary  current  and  cente 
of  the  gyre  at  four  different  times  during  spin-up. 
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at  steady  state  is  marked  by  perfect  homogenization,  there  is  no  input 
anywhere  around  the  gyre  nor  is  there  output. 

Two  different  cases  were  considered  in  some  detail  k  ~  lO^cm^/sec  and 

c  p 

<  "  5  x  10  cm  /sec  (hereafter  referred  to  as  and<^).  Figure  2.6 
shows  the  accumulation  of  tracer  at  the  center  of  the  gyre  for  the  two  cases. 
(The  interior  portion  of  the  gyre  is  the  last  place  in  the  domain  to  reach 
steady  state.)  There  are  two  things  to  note.  First,  it  took  over  three 
times  longer  for  steady  state  to  be  reached  in  the  smaller  <  example  (an  e- 
folding  time  of  ~80  years  versus  ~25  years).  Second,  the  value  of  the  homo¬ 
genized  pool  at  steady  state  is  the  same  for  both  cases.  (All  simulations 
were  halted  when  the  yearly  accumulation  rate  fell  below  2  percent  of  its 
earlier  maximum  rate.)  In  <5  tracer  readily  diffuses  off  the  boundary  and 
spin  up  occurs  relatively  quickly,  whereas  in  tracer  diffuses  from  the 
current  slowly  but  spin  up  occurs  much  later  and  tracer  accumulates  in  the 
gyre  for  a  longer  time.  It  is  not  obvious  why  these  effects  exactly  balance 
each  other  to  produce  the  same  level. 

In  their  work  on  homogenization  of  passive  tracers  in  gyres  Rhines  and 
Young  ( 198 J )  showed  that  the  time  scale  for  homogenization  to  occur  is  the 
diffusive  time  scale  of  the  gyre.  It  is  interesting  then  to  compare  esti¬ 
mates  based  on  this  to  the  e-folding  times  observed  above.  Since  the  e- 
folding  times  differ  in  the  two  examples  by  only  a  factor  of  J,  it  is  not 
sufficient  to  merely  consider  order  of  magnitude  estimates.  We  therefore 
make  use  of  results  from  the  following  example  of  tracer  diffusing  into  a 
closed  region. 

Consider  a  circular  domain  of  area  AQ  within  which  tracer  is  free  to 
diffuse.  The  edge  of  the  domain  is  maintained  at  a  constant  value  e  =  eQ 
and  initially  there  is  no  tracer  in  the  interior.  The  steady  state  distribu¬ 
tion  will  be  e  =  eQ  everywhere,  and  we  are  interested  in  how  long  it  takes 
to  reach  this  state.  The  governing  equation  is  the  radially  symmetric  diffu¬ 
sion  equation, 
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where  A  =  area  corresponding  to  a  given  radial  distance. 
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(The  equation  here  is  written  in  its  area  form  rather  than  in  terms  of  radial 
distance. ) 

This  problem  was  solved  by  Young  (personal  communication) ,  the  solution 

being 

_  t_ 

e(A,t)  =  eQ  Me  G>J0( 2/737^  , 

where  eQ  =  homogenized  value, 

7,  =  spin  up  time, 

and  JQ  =  zeroth  order  Bessel  function. 

The  expression  for  the  e-folding  time  is, 

( .  05) Aq 


In  our  case  the  gyre  is  not  circular  but  elliptical,  so  Aq  =  irL^, 

Lz  and  l_M  being  the  zonal  and  meridional  scales  of  the  domain.  Also, 
because  the  distance  between  streamlines  varies  around  the  gyre  (which  influ¬ 
ences  the  cross-stream  gradients),  <  is  replaced  by  an  effective  diffusiv- 
ity  V  The  details  of  why  this  is  so  are  given  in  Rhines  and  Young  (1983). 
In  that  work  they  show  that  the  elliptical  analog  to  (2.3)  is 

®t  =  *Dr  (d  Dr*  *  (2'4) 

where  A  =  area  within  a  streamline  defined  by  <P  =  constant, 
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r  =  jj  v  dA  =  circulation  around  a  streamline. 

A 

The  example  that  they  solve  is  an  elliptical  solid  body  flow  defined  by 
the  streamf unction 
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\l>  being  the  amplitude  of  the  flow,  L  and  L  the  x  and  y  length  scales. 


For  the  gyre  that  we  are  considering,  (1.11),  it  turns  out  that  the  expres¬ 
sion  for  D  is  exactly  the  same  as  for  the  solid  body  case. 
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With  this,  (2.4)  becomes 
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(2.5) 


Rhines  and  Young  call  the  quantity  in  brackets  the  effective  diffusivity 


<  ,  as  (2.5)  is  just  (2.3)  with  <  replaced  by  <  . 


The  expression  then  for  the  diffusive  time  of  the  elliptical  gyre  of 
the  model  i s 
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(The  term  inside  the  brackets  represents  the  deviation  from  the  simple  scale 
analysis  estimate.) 


The  scales  L  and  LM  define  the  size  of  the  domain  in  question,  so 
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the  pair  of  them  take  on  different  values  for  the  two  examples  being  consid¬ 


ered,  as  the  homogenized  region  is  smaller  for  larger  k.  In  homogeni¬ 

zation  occurs  to  roughly  the  ^  =  16  streamline,  while  in  only  to  the 
=  9  streamline  (Figure  1.17).  Substituting  the  appropriate  values  into 
(2.6)  gives  a  diffusive  time  of  5  years  for  and  .5  years  for  Com¬ 

paring  these  values  to  the  actual  times  we  see  that  in  the  first  case  spin  up 
is  15  times  longer  than  the  diffusive  time  of  the  gyre,  and  in  the  second 
case  it  is  50  times  longer. 

In  light  of  the  character  of  the  spin  up  process  that  occurs  in  the 
model,  this  result  is  not  surprising.  It  is  evident,  since  the  level  through¬ 
out  the  gyre  rises  uniformly,  that  the  rate  of  spin  up  is  being  controlled  by 


v.v.y /  .-.v:  ;  .  , 


' ' ' 


,  .  .  ■  .  i 


a  ^ 


the  diffusion  of  tracer  from  the  boundary  current  (i.e.  it  takes  little  time 
for  tracer  to  diffuse  throughout  the  gyre  once  it  reaches  the  edge).  In 
order  to  quantify  this  idea,  as  well  as  understand  what  factors  dictate  the 
level  of  the  homogeneous  pool,  a  simple  analytical  model  is  examined  below. 

Diffusive  Transfer  Model 

It  is  instructive  to  consider  the  domain  of  the  numerical  model  as 
being  composed  of  several  subregions.  The  motivation  behind  this  is  to  con¬ 
sider  regions  which  can  only  communicate  with  one  another  diffusively,  so  as 
to  create  a  simplified,  coarsely  resolved  system  with  one  less  degree  of 
freedom  (i.e.  no  advection).  The  idea  is  that  in  the  strongly  advective 
limit  of  the  numerical  model  it  is  the  slow  diffusive  processes  which  regu¬ 
late  the  spin  up  and  dictate  the  net  transfer  of  properties  into  the  inter¬ 
ior.  In  conjunction  with  the  simplified  domain  we  consider  a  simplified  spin 
up  process.  We  assume  that  the  plume  of  tracer  which  penetrates  the  gyre 
immediately  conforms  to  a  streamline,  i.e.  that  the  advective  limit  applies 
everywhere  in  the  gyre  (which  implies  homogenization,  as  discussed  in  chapter 
one) . 

Trie  four  regions  that  are  considered  are  depicted  in  Figure  2.6.  The 
first  region  corresponds  to  the  boundary  current.  Advection  from  the  north¬ 
ern  source  fills  this  region  quickly  and  continually  acts  to  maintain  the 
amount  of  tracer  so  contained  within  it.  The  second  region  is  the  outer 
strip  of  the  gyre  (subsequently  referred  to  as  the  edge  of  the  gyre).  Tracer 
diffuses  into  this  region  from  the  boundary  current,  and  quickly  gets  pulled 
around  the  circuit  forming  a  ridge.  From  here  the  tracer  proceeds  to  spread 
laterally,  diffusing  inward  to  the  center  portion  of  the  gyre  (the  third 
region)  and  diffusing  outward  to  the  vast  area  surrounding  the  gyre  (the 
fourth  region). 

We  call  this  model  the  diffusive  transfer  model.  Each  of  the  subre¬ 
gions  of  the  domain  is  represented  by  a  single  value  of  concentration,  and 
the  diffusive  transfer  between  them  is  character! zed  by  a  set  of  time  scales 
which  represents  the  amount  of  time  it  takes  tracer  to  diffuse  from  one 
region  to  another.  The  exchange  occurs  such  that  the  strength  of  the  flux 
is  proportional  to  the  difference  in  concentration.  (Note  that  the  regions 
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Figure  2.6:  Domain  of  the  diffusive  transfer  model  consisting  of  four 
regions.  Tracer  passes  into  adjacent  regions  by  way  of  diffusion  only 


expressed  as  such  do  not  have  "sizes";  this  information  is  contained  within 
the  time  scales.)  Because  of  the  strong  advective  input  into  the  boundary 
current  region,  its  value  is  fixed  throughout  the  spin  up.  Also,  we  take  tne 
fourth  region  to  be  an  infinitely  large  background  reservoir  whose  value  can¬ 
not  be  altered  and  so  is  set  equal  to  zero.  As  such  then  there  are  only  two 
active  regions  in  the  model,  those  which  comprise  the  gyre.  The  governing 
equations  for  the  gyre  edge  and  gyre  center  are 
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where  £>e  =  concentration  of  gyre  edge, 

=  concentration  of  gyre  center, 

*  =  concentration  of  boundary  current  region, 

tb  =  boundary  time  scale  (diffusion  between  boundary  current  and 
gyre  edge) , 

tg  =  gyre  time  scale  (diffusion  between  gyre  edge  and  center) , 

=  reservoir  time  scale  (diffusion  between  gyre  edge  and 
background  reservoir). 

The  initial  conditions  are  that  the  entire  gyre  is  free  of  tracer,  i.e. 

=  <Je  =  d  dt  t  = 

The  system  (2.7)  is  a  pair  of  ordinary  first  order  linear  differential 
equations  which  is  solved  straightforwardly  using  the  operator  method  (Zill, 
1979).  It  is  convenient  to  express  the  results  in  terms  of  the  nondimension- 
al  parameters  representing  the  relative  sizes  of  the  three  diffusive  time 
seal es. 
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It  is  not  obvious  a  priori  what  the  magnitudes  of  these  parameters  are. 
However,  the  results  of  the  numerical  model  can  be  used  to  constrain  them. 
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The  solutions  for  ri  and  t  are  of  the  form 

e  c 

t_  t_ 

"  T1  "  T2 

«S(t)  =  d0  +  <5xe  1  *  <s2e  c  (2.8) 


where  6Q  is  the  steady  state  value  and  the  other  two  terms  are  transients. 
Of  the  transients,  the  slower  mode  (i^)  dictates  the  spin  up  of  the  system. 
We  define  the  spin  up  time  then  as  t  =  In  terms  of  the  nondimensional 

parameters. 


T  =  S^g.A^)  TG 
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aGB  .  aGR  1 


1  +  2  *  2  '  I  ^GB^  +  AGR^  +  2&GBAGR  +  4 


(2.9) 


The  function  S  measures  how  long  the  spin  up  is  compared  with  the  diffusive 
time  of  the  gyre.  Recall  that  the  time  scale  over  which  homogenization 
occurs  in  the  numerical  model  is  many  times  greater  than  the  diffusive  time 
of  the  gyre  (t  >>  x^)  due  to  the  fact  that  spin  up  is  controlled  by  the 
diffusion  of  tracer  from  the  boundary.  This  (by  the  definition  of  agb) 
implies  that  agb  <<  1.  Consider  now  the  effect  that  varying  agr  has  on 
the  value  of  S  (keeping  A^g  <<  1). 


l-  agr  0 

This  condition  corresponds  to  the  case  when  there  is  no  background 
reservoir  (as  if  there  were  an  insulating  boundary  around  the  gyre).  When 
AqB  >  0  (2.9)  implies  that  S  »  However,  with  no  reservoir  it  is  evident 
that  in  steady  state  the  gyre  will  equilibrate  to  the  value  of  the  boundary 
current  region,  which  is  not  the  case  (Figure  1.17). 


2.  aqr  1  • 

This  says  that  as  tracer  progresses  around  the  «dge  of  the  gyre  it  just 
as  readily  diffuses  inward  as  outward.  In  this  case  as  agb  >  0  the 
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function  S  >  2.6,  i.e.  spin  up  can  at  most  be  roughly  two  and  a  half  times 
longer  than  the  diffusive  time  of  the  gyre.  However,  in  and 

in  <5  t  ~  50  so  this  case  is  unacceptable  as  well. 

The  reason  for  the  difference  in  the  limiting  value  of  S  for  the  two 
cases  can  be  understood  as  follows.  When  there  is  no  background  a  large 
amount  of  tracer  must  accumulate  in  the  gyre  (#c  >  $)  so  it  takes  a  long 

time.  However,  when  the  reservoir  is  characterized  by  a^r  ~  1,  the  corres¬ 

ponding  condition  that  a^b  >  0  implies  that  t  R  is  small,  i.e.  there  is 
strong  diffusion  into  the  reservoi r.  This  together  with  the  fact  that  the  res¬ 
ervoir  concentration  is  maintained  at  zero  (by  definition)  means  that  this 
region  acts  as  a  strong  sink.  Thus  only  a  small  amount  of  tracer  accumulates 

in  the  gyre  (tfc  <<  *)  and  so  spin  up  does  not  take  as  long.  Note  that  in 

the  extreme  limit  no  tracer  accumulates  in  the  gyre,  which  also  means  this  case 
is  unacceptable. 

3.  a^r  >>  1  >>  tr) . 

This  condition  leads  to  the  same  problems  as  in  case  2,  as  diffusion 
into  the  reservoir  is  now  even  stronger. 

We  see  then  that  the  condition  imposed  by  the  spin  up  implies  that  both 
A^g  <<  1  and  a^r  <<  1.  However  it  tells  us  nothing  about  the  relative  mag¬ 
nitudes  of  these  two  parameters.  For  this  information  we  examine  another 
result  of  the  numerical  model,  that  regarding  the  level  of  tracer  in  the 
gyre.  In  (2.8j  the  expression  for  for  the  center  of  the  gyre  is 
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The  function  P  represents  the  fraction  of  the  boundary  input  value  that  the 
gyre  equilibrates  to.  Note  that  P  depends  only  on  abr,  which  measures 
the  relative  importance  of  a^b  and  a^r.  We  saw  earlier  that  the  level  of 
the  homogeneous  pool  in  the  numerical  model  was  independent  of  diffusivity. 

In  terms  of  the  diffusive  transfer  model  this  suggests  that  fbn  not  depend 


must  take  on 
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on  any  of  the  diffusive  time  scales.  For  this  to  be  so 
the  same  value  in  each  example.  To  determine  this  value  we  can  match  (2.10) 
to  the  level  of  the  homogenized  pool. 

In  order  to  do  this  we  have  to  relate  the  boundary  value  <t>  to  an 
appropriate  quantity  in  the  numerical  model.  Note  that  the  set  up  of  the 
transfer  model  is  such  that  $  is  independent  of  diffusivity  as  well.  Since 
the  total  amount  of  tracer  in  the  boundary  current  region  varies  with  the 
size  of  K  (Figure  1.17)  we  take  the  peak  value  of  the  input  at  the  north¬ 
ern  boundary  as  the  value  for  i1  .  (The  transfer  is  based  on  the  presumption 
that  the  system  is  strongly  advective.  In  the  extreme  advective  limit  the 
entire  boundary  current  would  equilibrate  to  the  input  distribution,  in  which 
case  it  is  evident  that  the  peak  value  of  the  Gaussian  would  be  the  correct 
choice. ) 

The  homogenized  level  of  the  gyre  is  approximately  equal  to  one  half 
the  peak  value  of  the  input  (Figure  1.17)  which  implies  then  that  ~  1. 
This  condition  in  turn  says  that,  in  each  case,  the  time  scale  in  which  the 
edge  of  the  gyre  equilibrates  to  the  boundary  value  is  also  the  time  scale 
in  which  the  edge  of  the  gyre  decays  into  the  background  reservoir.  It  is 
not  obvious  why  this  is  so,  but  the  following  simple  argument  suggests  that 
thi s  i s  the  case. 

Consider  the  circular  domain  in  Figure  2.7,  which  is  meant  to  represent 
the  edge  of  the  gyre.  The  outer  perimeter  is  divided  into  two  parts:  the 
section  in  contact  with  the  boundary  and  that  in  contact  with  the  background 
reservoi r.  It  is  assumed  that  the  advection  serves  only  to  stir  up  the 
tracer  such  that  the  concentration  e  varies  only  in  the  radial  direction. 
The  governing  equation  is  the  radially  symmetric  diffusion  equation  (2.3) 
(written  in  radial  form) 
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Integrating  over  the  area  of  the  edge,  Ag, 
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where  the  bar  denotes  the  average  over  Ag.  The  first  term  on  the  right 
hand  side  is  the  contribution  associated  with  the  boundary  and  the  second  is 
that  associated  with  the  reservoir  (x  is  the  azimuthal  coordinate).  From 
(2.12)  it  follows  that 
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where  LD  and  L„  are  the  arc  lengths  of  the  boundary  and  reservoi r  sections 
respecti vely . 

The  expressions  inside  the  brackets  represent  the  flux  of  tracer  into 
and  out  of  the  edge  of  the  gyre  in  the  two  regions.  In  line  with  the  trans¬ 
fer  model  we  approximate  the  gradients  in  a  finite-difference  sense.  In 
particular,  where  the  gyre  is  adjacent  to  the  boundary 
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and  where  it  borders  the  reservoi  r 
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where  AOg,AOR  =  typical  change  in  concentration  from  boundary  to  edge, 

edge  to  reservoir, 

4Ly,aLR  =  boundary,  reservoir  length  scales  associated  with 
AOg.AOp  (Figure  2.7) . 

With  the  approximations,  (2.13)  becomes 
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from  which  it  is  evident  that 
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As  seen  in  Figure  1.17,  the  gradient  off  the  boundary  is  much  steeper 
than  the  gradient  from  the  gyre  to  the  reservoi r.  However  this  is  compen¬ 
sated  for  in  (2.16)  by  the  fact  that  most  of  the  gyre  borders  the  reservoir. 
Using  the  steady  state  distributions  to  estimate  ALg  and  aLr  (and  the 
streamline  pattern  to  compute  Lg  and  LR)  this  results  in  an  estimate 
for  AgR  of  .8  in  and  1.2  in 

Consideration  of  the  functions  S  and  P  thus  results  in  the  following 
ordering  of  the  diffusive  time  scales. 
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To  simplify  the  transfer  model  we  set  AgR  equal  to  the  constant  value  of 
one  (xg  =  tr )  which  in  effect  leaves  us  with  a  single  non-dimensional 
parameter 


(2.18) 


The  expressions  (2.9)  and  (2.10)  simplify  to 


S(  A)  = 
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(2. 19b) 


Tne  parameter  a  expresses  the  relative  strengths  of  the  two  funda¬ 
mental  time  scales  of  the  system.  We  can  now  quantitatively  explain  how 
diffusion  from  the  boundary  controls  the  spin  up.  As  depicted  in  Figure  2.8 
when  a  >  0,  t  >  x^  and  the  system  spins  up  on  the  boundary  time  scale 


a 

$ 

'wmmm. 


c>r 

'  V. 


83 


(once  tracer  diffuses  from  the  boundary  to  the  edge  of  the  gyre  it  quickly 
spreads  into  the  center).  This  is  the  case  in  the  numerical  model.  On  the 
other  hand,  when  a  >  00 ,  x  >  which  corresponds  to  the  Rhines  and  Young 
example  in  which  spin  up  occurs  on  the  gyre  time  scale.  In  this  case  the 
edge  of  the  gyre  becomes  equilibrated  quickly,  and  tracer  then  proceeds  to 
diffuse  into  the  center.  By  matching  (2.19a)  to  the  numerical  simulations 
using  the  gyre  diffusive  times  estimated  earlier  for  the  and  exam¬ 
ples,  we  determine  Aj  and  Ag  such  that  the  value  of  S(a)  matches  the 
associated  t.  (The  magnitude  of  a  is  determined  in  this  fashion  rather 
than  from  the  definition  (2.18)  because  it  is  uncertain  what  value  of  Ag 
to  use  in  (2.15a).)  With  this  information  then  we  can  use  the  diffusive 
transfer  model  to  predict  the  amount  of  tracer  at  the  center  of  the  gyre  in 
the  numerical  model  as  a  function  of  time. 

It  was  discussed  aoove  that  both  i  and  6  are  sums  of  a  steady  term 

(«  C 

and  two  transient  modes.  It  is  the  case  that  in  the  small  a  limit  which 
applies  here  the  faster  transient  decays  so  quickly  relative  to  the  slower 
one  that  it  can  be  ignored.  Consider  the  function  R1(a),  defined  as  the 
ratio  of  the  time  scales  of  the  two  transient  modes  in  (2.8), 
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This  measures  the  relative  importance  of  the  two  modes.  The  graph  of  R1 

2 

versus  a  is  plotted  in  Figure  2.9  and  for  a  >  0,  R1  >  »  ( ri  £  £) .  Note 
also  that  when  a  >  00  the  same  mode  still  dominates  ( R1  %  2a).  The  smallest 
amount  of  discrepancy  between  the  decay  time  of  the  two  modes  occurs  when 
A  ~  1  and  ^  ~  5  T^. 

This  means  that  in  the  a  <<  1  limit  the  time  dependent  part  of  i>c 
(and  <ie)  is  approximately  a  single  exponential,  and  suggests  that  we  re¬ 
draw  Figure  2.5  and  normalize  the  time  axis  by  the  appropriate  spin  up  times 
of  the  two  examples.  This  is  done  in  Figure  2.10.  Also  shown  in  the  figure 
are  the  corresponding  6C  curves  from  the  diffusive  transfer  model,  normal¬ 
ized  in  a  similar  fashion.  Note  that  the  and  Kg  examples  collapse  to  a 
single  case.  They  agree  reasonably  well  with  the  predictions.  (Keep  in  mind 


Figure  2.9:  The  relative  decay  of  the  two  transient  modes  of  the  trans 
model  solutions.  The  faster  mode  can  be  Ignored  for  both  large  and  sma 
values  of  &. 
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Figure  2.10:  The  accumulation  curves  of  Figure  2.5  only  the  time  axis  has 
been  normalized  for  each  curve  by  the  associated  spin-up  time.  Tne  solid 
lines  are  the  predictions  for  the  gyre  center  of  the  transfer  model. 


that  the  transfer  model  is  a  crude  simplification,  meant  only  to  reveal  some 
of  the  prominent  diffusive  processes  at  work.)  Not  shown  in  the  figure  are 
the  associated  curve:  which  closely  resemble  those  for  the  gyre  center. 
This  is  consistent  with  the  numerical  result  that  the  level  of  tracer 
throughout  the  whole  gyre  rises  at  nearly  the  same  rate  (as  t  >•  »,  > 

6c  and  complete  homogenization  is  achieved). 

B)  With  Vertical  Flux 

We  now  consider  the  presence  of  vertical  mixing  and  explore  how  this 
modifies  the  various  results  that  have  been  obtained  thus  far.  This  corres¬ 
ponds  to  the  oxygen  case,  and  in  the  numerical  model  we  now  include  the  radi¬ 
oactive  decay  term  in  (2.2).  The  boundary  conditions  are  identical  to  those 
in  the  salt  case  as  we  take  eu  =  =  0  (i.e.  the  vertical  anomaly  is  just 

equal  to  the  concentration  in  the  center  of  the  layer).  Two  different  simu- 

2 

lations  were  done,  one  with  v  =  1  cm  /sec,  H  =  1000  m  and  the  other  with 
2 

v  =  2  cm  /sec,  H  =  750  m.  These  will  be  referred  to  as  v,  and  v?  respec- 

6  2  ^  ^ 
tively  (in  both  cases  *  =  10  cm  /sec). 

Figure  2.11  shows  a  snapshot  of  the  spin  up  of  .  The  first  thing  to 
notice  is  that  it  would  be  very  difficult  to  distinguish  between  this  distribu¬ 
tion  and  a  similar  snapshot  from  the  salt  case  (i.e.  from  the  spin  up  of  k^. 
The  plume  of  tracer  extends  from  the  boundary  and  spirals  into  the  gyre  in 
the  same  manner.  A  marked  difference  between  the  salt  and  oxygen  cases 
occurs  only  in  the  context  of  steady  state.  The  steady  state  distribution  of 
the  \>2  experiment  appears  in  Figure  2.12a,  this  is  to  be  compared  with 
Figure  1.17a  for  which  the  vertical  flux  is  zero.  The  most  striking  thing  in 
regards  to  such  a  comparison  is  that  the  level  of  tracer  in  the  gyre  in  the 
non-zero  vertical  flux  case  is  substantially  smaller.  This  difference  is 
highlighted  in  Figure  2.12b  showing  the  accumulation  of  tracer  in  the  gyre 
for  the  two  cases.  The  reason  for  this  discrepancy  is  straightforward  in 
that  for  oxygen  there  is  an  additional  mechanism  present  for  removing  tracer 
from  the  domain  (this  is  also  the  cause  for  the  decrease  in  time  required 
for  spin  up). 

Let  us  contrast  the  spin  up  of  v ^  to  that  of  The  concentration 

of  tracer  in  the  boundary  current  is  nearly  identical  to  that  when  there  is 
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Figure  2.11:  Instantaneous  distribution  of  tracer  during  spin-up,  for 
experiment  which  includes  vertical  mixing. 
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(b)  Time  history  of  tracer  accumulating  at  the  center  of  tne  g.y 
10°cm^/sec  comparing  the  case  with  no  vertical  mixing  to  tf>jt 
i  cnr/sec. 
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no  vertical  flux.  This  is  because  the  strong  flow  of  the  current  does  not 
allow  the  fluid  to  spend  enough  time  in  this  region  to  be  significantly  in¬ 
fluenced  by  vertical  mixing.  It  is  while  the  tracer  diffuses  laterally  from 
the  boundary  that  vertical  exchange  becomes  prominent,  acting  as  a  sink.  The 
small  amount  of  tracer  that  does  penetrate  into  the  gyre  does  so  in  the  same 
fashion  as  in  K^.  As  steady  state  is  approached  though,  whereas  in  K ^ 
the  plume  extending  from  the  boundary  becomes  much  less  pronounced,  here  it 
remains  prominent.  Along  with  this,  the  slight  bowl  shaped  tracer  distribu¬ 
tion  in  the  gyre  gives  way  to  homogenization  in  but  remains  a  permanent 
feature  in  v.,  (Figure  2.13). 

At  steady  state  there  can  be  no  net  flux  of  tracer  into  the  region 
bounded  by  a  streamline  of  the  gyre.  In  this  is  accomplished  in  the 
center  part  of  the  gyre  through  homogenization.  In  the  outer  part,  the  rem¬ 
nant  of  the  spiral  weakly  pumps  tracer  inward  in  the  southern  and  eastern 
regions,  and  outward  to  the  north  and  west.  In  the  oxygen  case  there  is  a 
persistent  vertical  flux  out  of  the  region  bounded  by  a  streamline  (true  for 
every  streamline)  which  must  be  balanced  by  an  inward  lateral  flux.  One  way 
to  view  why  the  plurne  of  salt  is  less  pronounced  than  that  for  oxygen  at 
steady  state  is  to  recall  that  each  successive  plume  leaving  the  boundary 
current  does  so  from  further  downstream.  This  tends  to  smear  the  intrusion 
as  time  progresses.  For  oxygen  though  spin  up  occurs  more  quickly  so  that 
this  effect  is  not  at  work  for  as  long. 

It  is  interesting  that  in  v2,  which  contains  a  rather  large  vertical 

diffusion  parameter  the  penetration  of  tracer  into  the  gyre  so  closely 

hr 

resembles  that  which  occurs  in  (which  has  no  vertical  flux).  This  is 
due  to  the  fact  that  so  little  tracer  enters  the  gyre  that  the  vertical  anom¬ 
aly  is  small,  hence  weak  flux.  Consider  the  steady  state  balance  of  terms  in 
(2.2)  that  applies  in  the  northward  flow  of  the  outer  part  of  the  gyre  in  v2, 
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where  it  is  assumed  that  alongstream  diffusion  is  negligible  with  respect  to 
cross-stream  diffusion.  In  this  region  the  spiral  is  significant  (recall 
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Figure  2.13:  Zonal  section  of  tracer  through  the  boundary  current  and  center 
of  the  gyre  for  the  steady  state  distributions  of  Figure  1.17a  (<  - 
IQ^citP/sec)  and  2.12a  {<  ~  lO^cn^/sec,  v  ~  2  cn^/sec)  Illustrating  the 
effects  of  vertical  mixing. 


that  In  the  example  analyzed  In  chapter  one  k  ~  10  cm  /sec  as  it  is  in  v^) 
so  that  the  cross-stream  length  scale  is  "L^.  The  size  of  the  terms  in 
(2.20)  are  then  estimated  as  follows, 
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where  V  =  representative  velocity 

L.  =  alongstream  length  scale. 

a 

Plugging  in  values  for  the  scales  from  \>g  it  is  found  that  the  first  two 
terms  of  (2.21)  balance  each  other  to  within  5  percent,  which  is  the  same 
balance  that  characterized  the  spiral  in  the  salt  case. 

This  result  has  an  Interesting  consequence.  It  says  that  a  point 
balance  in  this  region  is  one  of  advection  versus  lateral  diffusion.  From 
knowledge  of  the  various  scales  then  we  can  produce  an  estimate  of  the 
lateral  diffusivity. 


(2.22 


However,  we  know  that  in  terms  of  an  integrated  balance  within  a  streamline 
that  the  lateral  flux  of  tracer  inward  balances  the  vertical  flux  outward 
(the  contribution  from  advection  identically  vanishes).  In  particular, 
integrating  the  steady  form  of  (2.2)  within  a  streamline  and  applying  the 
divergence  theorem  gives. 
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where  A  »  area  enclosed  within  the  streamline  S, 
n  *  normal  to  the  streamline. 

Having  estimated  *  then  from  (2.22),  we  can  In  turn  use  this  value  in 
(2.23)  to  determine  v. 


93 


(2.24) 


This  procedure  was  applied  to  the  experiment  to  try  and  recover 

the  values  of  lateral  and  vertical  dlffuslvlty  from  knowledge  of  tne  steady 

state  tracer  field  and  the  velocity  field.  The  values  so  obtained  were  <  ~ 

6  2  2 

1.6  x  10  cm  /sec,  v  ~  2.3  cm  /sec  (actual  values  In  v9  were  1.0  x 

q  2  2  ^ 

10  cm  /sec  and  2.0  cm  /sec,  respectively). 

Transfer  Model  with  Vertical  Flux 

The  diffusive  transfer  model  can  also  be  modified  to  include  the 
effects  of  vertical  diffusion.  The  areas  directly  above  and  below  the  deep 
layer  (represented  by  ey  and  In  the  numerical  model)  can  be  thought  of 
as  together  forming  another  subregion  of  the  domain,  consistent  with  the 
premise  of  the  model  that  this  region  be  In  contact  diffusively  with  adjacent 
regions.  The  reasons  for  considering  this  additional  case  are  the  same  as 
before,  to  help  more  fully  understand  the  numerical  model  results  by  gaining 
insight  Into  how  the  various  diffusive  processes  Interact,  which  in  turn  will 
be  useful  when  considering  the  actual  oceanic  case.  For  Instance,  we  saw 
that  the  addition  of  vertical  flux  Influenced  spin  up  times,  as  well  as  the 
amount  and  distribution  of  tracer  in  the  gyre.  It  Is  desirable  to  quantify 
these  effects. 

Tne  specifics  of  the  expanded  transfer  model  are  as  follows.  The 
additional  region  Is  in  contact  with  only  two  of  the  previous  ones,  those 
comprising  the  gyre  (i.e.  the  edge  and  center).  The  reason  why  It  is  not  in 
contact  with  the  boundary  current  region  is  that  we  assume  the  tracer  flushes 
through  this  region  too  quickly  to  be  altered  by  vertical  diffusion.  It  Is 
not  in  contact  with  the  background  reservoir  for  the  trivial  reason  that  the 
reservoir  is  void  of  tracer.  Because  the  new  region  represents  a  vast  area 
as  well,  it  too  is  assumed  to  be  a  reservoir  whose  value  cannot  be  changed 
and  is  thus  set  equal  to  zero.  From  here  on  we  must  differentiate  between 
the  lateral  background  reservoir  and  the  vertical  background  reservoir. 
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As  such,  although  we  are  including  another  region  it  is  not  an  active 
one,  and  consequently  we  are  not  faced  with  a  third  governing  equation.  We 

H2 

do  however  need  to  introduce  another  time  scale,  t  =  -g-  ,  which  is  the 

V  ov 

time  it  takes  tracer  to  decay  into  the  vertical  reservoir.  The  appropriate 
set  of  equations  is, 
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(2.25) 


Note  that  we  have  set  xR  =  Tg  as  before.  (The  relation  (2.17)  must  still 
apply.)  The  initial  conditions  are  the  same  as  before,  dc(0)  *  dg(0)  *  0. 
When  ty  >  »  (no  vertical  reservol r)  the  set  (2.25)  collapses  to  (2.7). 

As  before.  It  is  convenient  to  discuss  results  in  terms  of  nondimen- 
slonal  parameters.  With  the  addition  of  xy  we  introduce  another  such 

T(i 

parameter,  e  =  —  ,  in  addition  to  a  (=  — )  defined  earlier.  Since  both 
Tv  tB 

6  2 

of  the  numerical  runs  with  vertical  diffusion  had  k  ~  10  cm  /sec  (the  value 
in  <j)  we  set  a  »  a^.  This  way  we  Isolate  the  effect  of  vertical  flux 
within  the  context  of  the  numerical  results  previously  analyzed.  Solutions 
of  (2.25)  are  thus  presented  In  terms  of  the  single  parameter  e,  which  mea¬ 
sures  the  relative  strength  of  lateral  diffusion  Into  the  gyre  versus  ver¬ 
tical  diffusion  out  of  the  gyre. 


1.  Spin  up  time. 


Because  the  expanded  transfer  model  does  not  require  a  third  governing 
equation,  the  solutions  for  and  i  are  still  of  the  form  (2.8),  with 

v  v 

two  transient  modes.  The  spin  up  time  is  defined  as  before  (the  time  scale 
of  the  slower  mode) , 
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where  tq(a^)  *  spin  up  time  of  the  system  when  there  Is  no  vertical  flux, 
evaluated  at  a  =  a^  We  saw  earlier  that  to(aj)  £  Ta,  so  the  function 
U  1)  measures  how  much  more  quickly  the  system  spins  up  relative  to  the 
boundary  diffusive  time.  The  function  S2  U  1)  measures  how  much  smaller 
the  spin  up  time  Is  relative  to  the  other  fundamental  time  scale,  the  vertic- 
al  diffusive  time.  In  the  limit  e  >  0,  which  characterizes  weak  vertical 
mixing,  t  >  t  in  the  opposite  limit  e  >  “  (strong  vertical  mixing), 

*  *  Ty.  The  graph  of  Sj  versus  e  appears  in  Figure  2.14. 

As  was  the  case  with  no  vertical  flux,  under  certain  conditions  the 
solutions  for  i  and  i  can  be  approximated  by  the  slower  transient  mode 

»  V 

(In  addition  to  the  steady  term).  Consider  the  function  which  measures  the 
relative  decay  of  the  two  transient  modes  according  to  the  strength  of  ver¬ 
tical  mixing  (the  analog  of  Rl), 


R2(  e ) 
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When  the  vertical  mixing  is  weak  U  ►  0)  the  faster  decaying  mode  can  be 
Ignored  (in  this  case  R2  collapses  to  Rl  evaluated  at  a  =  a^  <<  1,  and 
recall  that  for  small  a  the  single  mods  approximation  Is  a  good  one).  How¬ 
ever,  as  seen  in  Figure  2.15,  as  the  vertical  mixing  Increases  the  single 
mode  approximation  steadily  worsens,  and  in  the  limit  where  vertical  mixing 
dominates,  the  approximation  is  Invalid. 
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Figure  2.14:  Comparison  of  the  extent  to  which  the  transfer  model  spin-up 
time,  final  gyre  level,  and  homogenization  are  affected  by  the  presence  of 
vertical  mixing.  The  strength  of  the  mixing  Is  measured  by  the  size  of  the 
parameter  e.  The  value  of  e  associated  with  experiment  and  that 
with  experiment  are  Indicated  by  ej  and  z^  respectively. 


2. 


Level  in  gyre. 

The  other  major  difference  for  the  case  of  non-zero  vertical  flux,  in 
addition  to  the  shortened  spin  up  time,  is  the  reduced  level  of  tracer  in  the 
gyre  at  steady  state.  Recall  that  when  vertical  flux  is  absent  the  final 
level  is  independent  of  the  time  scales  (or  in  terms  of  the  numerical  model, 
the  lateral  diffusivity) .  It  happens,  however,  that  when  vertical  flux  is 
made  different  from  zero  (i.e.  Tv  becomes  finite)  not  only  is  this  level 
dependent  on  the  strength  of  that  flux,  it  becomes  dependent  on  the  strength 
of  the  lateral  fluxes  as  well.  In  the  steady  state,  the  expression  for  ic 
is. 
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The  function  f  U  1)  measures  the  fraction  of  tracer  in  the  gyre  relative 
to  how  much  would  oe  there  if  the  vertical  flux  were  identically  zero.  The 
function  is  plotted  versus  e  in  Figure  2.14. 


3.  Homogenization 

In  the  discussion  of  the  previous  transfer  model  it  was  mentioned  that 
the  evolution  of  was  nearly  identical  to  that  of  dc,  and  that  in  the 
large  time  limit  i  >  i  (i.e.  complete  homogenization  occurs).  In  light 

C  L 

of  the  numerical  model  results  with  vertical  mixing,  it  is  not  surprising 
that  this  no  longer  is  the  case  here.  In  the  final  state  of  v2  the  dis¬ 
tribution  of  tracer  within  the  gyre  is  characterized  by  a  slight  minimum  at 
the  gyre  center.  The  equivalent  to  this  in  the  transfer  model  is  the  fact 
that  *L(t  >  00 )  <  d  ( t  >  “) .  As  a  measure  of  how  significant  this  feature  is 

V  “ 

in  the  transfer  model  we  define  the  function  h,  which  compares  the  gradient 
of  tracer  off  of  the  boundary  to  the  gradient  across  the  gyre. 


The  plot  of  h  versus  e  appears  in  Figure  2.14. 

It  is  interesting  to  note  in  Figure  2.14  that  whereas  the  spin  up  time 
and  level  of  tracer  in  the  gyre  are  extremely  sensitive  to  the  strength  of 
vertical  mixing,  the  extent  of  homogenization,  as  measured  by  h,  remains 
basically  the  same,  i.e.  the  gyre  is  nearly  homogenized  regardless  of  the 
size  of  .  Closer  inspection  of  h  shows  that  the  point  at  which  the 
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gyre  is  least  homogenized  occurs  for  a  finite  value  of  e  (Figure  2.16).  The 
reason  for  this  is  straightforward.  At  very  small  values  of  the  diffusion 
parameter  ,  the  vertical  flux  is  so  weak  that  the  system  behaves  as  in 

H  8v 

the  salt  case  and  homogenization  occurs.  At  very  large  values  of  — X  ,  the 
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strong  flux  causes  the  vertical  reservoir  to  draw  away  nearly  all  the  tracer 
that  reaches  the  edge  of  the  gyre  (having  diffused  from  the  boundary).  Thus, 
very  little  tracer  ends  up  in  the  edge  of  the  gyre  and  even  less  in  the 
center,  so  homogenization  again  occurs  (in  the  sense  that  the  gyre  is  uni¬ 
formly  void  of  tracer).  The  point  at  which  the  gyre  is  least  homogenized 
then  occurs  between  these  limits. 

Having  reviewed  these  effects  of  vertical  flux  in  terms  of  the  param¬ 


eter  e,  let  us  see  where  in  this  domain  the  two  experiments  and  v2 
fall.  Using  the  previously  computed  estimate  of  tg  and  the  appropriate 
values  of  v  and  H,  we  find  that  corresponds  to  moderate  vertical 
mixing  and  v2  corresponds  to  strong  vertical  mixing  (the  values  of  e  so 
computed,  e1  and  e2,  are  marked  in  Figure  2.14).  As  was  done  earlier,  we 
can  use  this  knowledge  concerning  the  time  scales  to  predict  the  spin  up  of 
the  gyre  in  the  numerical  experiments.  Such  a  prediction  for  the  v2 
example  appears  in  Figure  2.17,  where  the  concentration  of  the  gyre  center  is 
compared  with  when  a  =  a^,  e  =  e2.  Also  shown  are  the  similar  compar¬ 
isons  for  < ■  and<g.  It  is  seen  that  for  this  case  as  well  the  transfer 
model  provides  a  reasonable  fit  to  the  numerical  data. 

There  are  other  interesting  differences  when  vertical  flux  is  added  to 
the  system  in  addition  to  those  discussed  above.  Some  of  these  are 


Figure  2.16:  Enlargement  of  the  graph  of  h  in  Figure  2.14  showing  the 
degree  of  homogenization  versus  the  strength  of  the  vertical  mixing.  Th< 
gyre  is  least  homogenized  for  e  ~  .23. 


Figure  2.17:  Comparison  of  the  accumulation  of  tracer  at  the  center  of  the 
gyre  for  the  different  values  of  lateral  and  vertical  diffusivity.  The  solid 
lines  are  the  corresponding  predictions  from  the  transfer  model. 
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illustrated  in  Figure  2.18,  within  the  context  of  the  transfer  model.  When 
there  is  no  vertical  reservoir  the  steady  state  path  of  tracer  is  from  the 
boundary  to  the  lateral  background  via  the  edge  of  the  gyre.  This  path  be¬ 
comes  more  complex  when  vertical  flux  is  introduced.  For  instance,  tracer  is 
then  fluxed  into  and  out  of  the  gyre  center,  although  the  strength  of  this 
flux  is  less  than  that  tnrougn  the  edge  of  the  gyre  (which  becomes  stronger 
with  increased  vertical  mixingj.  Also  illustrated  is  the  manner  in  which  the 
flux  away  from  the  edge  of  the  gyre  is  partitioned.  In  particular,  when  ver¬ 
tical  mixing  is  weak  the  flux  is  nearly  all  to  the  lateral  reservoir,  but 
when  vertical  mixing  is  strong  the  flux  is  equally  partitioned  between  the 
vertical  reservoir  and  the  center  of  the  gyre.  Mote  then  that  with  a  strong¬ 
er  vertical  flux  more  tracer  diffuses  laterally  from  the  boundary,  but  less 
makes  it  into  the  interior. 

C)  Conclusions 

Two  types  of  tracers  have  been  studied  in  some  detail  using  a  simple 
numerical  model  that  was  designed  to  represent  the  NRG/DWBC  system.  The 
first  type,  of  which  salinity  is  an  example,  is  governed  by  lateral  processes 
only  whereas  the  second  type,  such  as  oxygen,  is  influenced  by  vertical  mix¬ 
ing  as  well.  The  model  produced  various  interesting  results,  among  which  is 
the  fact  that  when  vertical  mixing  is  absent  the  final  level  of  tracer  in  the 
gyre  seems  to  be  independent  of  the  strength  of  the  lateral  mixing  (provided 
it  is  weak).  It  also  provided  a  method  for  recovering  the  values  of  the  eddy 
di ffusi vi ties.  Interpretation  of  the  results  though  became  clearer  when  an 
analogous  diffusive  transfer  model  was  applied  to  tiie  numerical  model. 

The  transfer  model  was  able  to  explain  quantitatively  what  factors 
determine  how  much  tracer  diffuses  into  the  gyre  from  the  boundary  current. 

In  particular,  it  led  to  an  expression  for  the  gyre  level  which  includes  the 
core  concentration  of  the  boundary  current  and  the  lateral  and  vertical  dif- 
fusivities.  When  the  vertical  mixing  is  identically  zero,  all  dependence  of 
the  level  on  the  strength  of  the  lateral  diffusion  drops  out  as  well,  pro¬ 
vided  the  nondimensional  parameter  ^  =  Tg/T^  remains  constant.  It  is  not 
immediately  obvious  why  this  happens  in  the  numerical  model,  and,  as  such, 
the  transfer  model  does  not  offer  an  explanation  as  to  why  the  level  is  inde¬ 
pendent  of  <  but  rather  explains  what  factors  are  involved. 


Figure  2.18:  The  effect  of  vertical  mixing  on  the  steady  state  flux  of 
tracer  between  the  regions  of  the  transfer  model. 

(a)  Comparison  of  the  flux  into  and  out  of  the  edge  of  the  gyre  versus 
the  center  of  the  gyre.  The  non-dimensional  functions  fp  and  fc  measure 
these  quantities  respectively,  relative  to  the  value  for  the  edge  wnen 
vertical  mixing  is  absent. 


The  transfer  model  does,  however,  nicely  explain  various  other  features 
of  the  numerical  results.  The  presence  of  homogenization  within  the  gyre  is 
a  rigorous  feature  of  the  system.  A  nearly  flat  plateau  develops  regardless 
of  the  strength  of  the  vertical  dlffuslvity  (provided  the  lateral  diffusivity 
is  small  enough).  The  spin  up  of  the  system  is  slower  than  might  be  suggest¬ 
ed  by  the  size  of  the  homogenized  region  in  the  different  examples.  This  is 
because  the  spin  up  is  dictated  by  the  slow  diffusion  of  tracer  from  the 
boundary  current  to  the  edge  of  the  gyre  (modified  by  diffusion  vertically 
out  of  the  domain).  In  addition  to  illuminating  some  of  the  results,  as  seen 
below  the  transfer  model  is  also  necessary  in  order  to  apply  them  to  the  data. 


Data  Comparison 

A  major  objective  of  the  numerical  model  study  (including  the  transfer 
model)  was  to  use  the  acquired  results  to  address  the  idea  that  this  type  of 
advective-diffusive  situation  exists  In  the  ocean.  Specifically,  does  tracer 
data  from  the  region  support  the  existence  of  the  northern  recirculation 
gyre?  In  this  section  we  examine  the  data  from  this  perspective.  From  the 
model  we  have  learned  what  things  to  be  looking  for  in  the  data,  how  to  in¬ 
terpret  various  features,  and  now  to  extract  specific  information  from  these 
features,  He  now  consider  in  some  detail  the  tracer  distributions  from  the 
data  set  described  in  Hogg  et  al .  (1986). 

The  first  thing  to  do  is  review  general  characteristics  and  see  if  they 
are  consistent  with  what  is  believed  to  be  the  flow  field.  This  was  done  in 
Hogg  et  al .  (1986).  The  distinguishing  feature  of  closed  circulation,  a 
region  of  homogenization,  was  present  in  all  of  the  various  tracer  fields. 

Here  we  comment  further  on  that  analysis.  The  main  thrust  of  this  section, 
however,  is  a  more  detailed  look  at  some  of  the  specific  features  of  the  data. 

In  the  model  study  we  needed  to  specify  boundary  conditions  and  certain 
parameters  (i.e.  the  diffusivities) .  In  order  to  apply  some  of  the  results 
so  obtained,  we  need  then  to  identify  the  analogous  boundary  conditions  in 
the  data.  It  is  results  of  this  type  which  are  examined  (for  instance, 
determining  how  much  of  a  given  tracer  has  diffused  from  the  boundary  current 
into  the  gyre).  As  a  consequence,  we  acquire  knowledge  regarding  the  oceanic 
values  of  the  diffusive  parameters. 


A  detailed  description  of  the  data  set  appears  In  Hogg  et  al .  (1986) . 
The  data  Includes  tracer  measurements  from  various  cruises  In  the  region 
north  of  35*N,  west  of  55  *W,  and  south  of  the  2000  m  Isobath  of  the  contin¬ 
ental  slope  (there  are  numerous  crossings  through  the  DWBC,  Figure  2.19). 

Most  of  the  data  comes  from  a  single  cruise,  QCEANUS  134,  which  was  part  of 
the  Abyssal  Circulation  Experiment  (1983-84)  designed  to  study  the  NRG.  The 
analysis  centers  on  the  distribution  of  oxygen,  which  is  assumed  to  be  in  a 
nearly  steady  state.  For  this  reason  the  results  from  the  model  regarding 
spin  up  are  not  directly  applicable.  We  will  discuss  in  turn  some  of  the 
other  Important  features. 

A)  Results 

1.  Homogenization 

We  consider  this  feature  first,  as  this  was  already  addressed  In  Hogg 
et  al.  (1986).  To  detect  homogenization  requires  only  Information  about  rel¬ 
ative  magnitudes  of  tracer  concentrations;  thus  we  need  not  concern  ourselves 
with  boundary  conditions  as  yet.  In  order  to  reduce  the  noise  level  In  the 
data,  Hogg  et  al .  employed  a  vertical  averaging  scheme.  For  each  station  In 
Figure  2.19  what  Is  desired  is  the  value  of  trace*-  at  the  density  level  cor¬ 
responding  to  the  core  of  the  DWBC  (consistent  with  the  notion  of  isopycnal 
spreading).  Rather  than  use  this,  Hogg  et  al .  used  the  average  concentration 
between  two  bounding  density  levels  of  the  core  (specifically,  an  anomaly 
averaged  over  the  deep  layer,  defined  here  as  *av) •  They  presented  lateral 

maps  of  e  for  various  tracers, 
av 

The  maps  for  oxygen  appears  In  Figure  2.20.  It  is  characterized  by  a 
strong  gradient  off  the  boundary  that  merges  into  a  broad  region  which  is 
nearly  uniform.  Further  to  the  south  the  level  once  again  drops  off.  Note 
the  plume  of  high  oxygen  that  extends  off  of  the  boundary  and  well  into  the 
interior.  From  this  map,  and  a  similar  one  for  salinity,  Hogg  et  al.  com¬ 
puted  a  corresponding  average  section  across  the  gyre  through  the  boundary 
current.  These  are  shown  in  Figure  2.21.  Hogg  et  al.  discussed  the  fact 
that  although  both  tracers  exhibit  homogenization,  the  oxygen  distribution 
has  a  slight  minimum  within  this  region,  suggestive  of  vertical  mixing. 
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Figure  2.19:  Station  positions  of  the  tracer  measurements,  from  six  differ¬ 
ent  cruises,  used  in  the  data  analysis  (from  Hogg  et  al.,  1986) . 
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Superimposed  on  the  lateral  map  of  oxygen  are  two  "streamlines."  The 
outer  one  was  constructed  using  the  closed  deep  layer  isopach  in  this  region 
(from  Hogg  and  Stommel ,  1985)  as  a  guide.  The  inner  one  was  constructed  to 
fit  entirely  within  the  region  covered  by  the  tracer  data.  Note  that  the 
oxygen  plume  emanates  from  the  boundary  beyond  the  region  where  the  stream¬ 
lines  turn  offshore.  This  Is  entirely  consistent  with  the  numerical  model 
results  and  further  supports  the  claim  for  a  tight  reel rculation.  The  ridge 
of  high  oxygen  then  proceeds  to  penetrate  into  the  gyre  in  a  manner  similar 
to  the  model.  Note  that  the  fact  that  this  ridge  is  very  slight  does  not 
imply  that  vertical  mixing  is  weak,  the  transfer  model  shows  that  any  ver¬ 
tical  mixing,  weak  or  strong,  exhibits  a  subtle  feature  such  as  this. 

2.  Flux  balances 

We  saw  earlier  that  the  point  balance  of  terms  in  (2.20)  in  the  region 
of  the  spiral  led  to  an  estimate  of  the  lateral  diffusivity.  In  the  oxygen 
distribution  of  Figure  2.20  we  see  evidence  of  such  a  spiral  in  the  data. 
However,  it  should  be  noted  that  the  distribution  is  a  synoptic  one,  and, 
this  being  an  area  characterized  by  strong  time- dependent  fluctuations,  such 
a  distribution  is  sure  to  differ  somewhat  from  the  average  pattern.  For 
this  reason,  applying  a  point  balance  calculation  involves  a  large  amount  of 
uncertainty,  and  therefore  we  do  not  make  use  of  (2.22). 

The  other  flux  balance  that  was  discussed  was  an  integrated  balance 
within  a  streamline  of  the  gyre.  The  lateral  flux  of  tracer  into  the  gyre 
matches  the  vertical  flux  Into  the  upper  and  lower  layers.  This  equality 
led  to  an  estimate  of  the  vertical  diffusivity. 


Since  (2.2 8)  Involves  spatial  integrals,  this  removes  some  of  the  noise  due 
to  transients  so  the  balance  can  be  applied  to  the  oxygen  distribution  with 
some  confidence.  In  (2.28)  recall  that  the  value  of  e  is  actually  the 


«u  *  ®L 

vertical  anomaly,  e  =  - 2 -  (ey  and  are  the  respective  values  at 

the  top  and  bottom  of  the  deep  layer,  ^  is  the  value  in  the  center  of  the 
layer).  In  the  model  study  we  set  ®u  =  ®l  =  °»  but  i n  the  ocean  this  is 
certainly  not  the  case. 

Consider  again  the  oxygen  profile  versus  depth  in  Figure  2.1  for  a  sta¬ 
tion  in  the  center  portion  of  the  gyre.  The  profile  is  asymmetric  in  that 
the  value  near  the  top  of  the  deep  signal  is  smaller  than  that  below  the  deep 
signal.  This  is  true  of  most  of  the  stations  in  the  domain.  As  a  result, 
the  flux  out  of  the  top  of  the  deep  layer  is  greater  than  the  flux  out  of  the 
bottom,  and  in  fact  this  downward  flux  can  be  ignored.  In  terms  of  (2.2), 
when  there  is  upward  flux  only  the  value  of  the  vertical  mixing  parameter 
decreases  by  a  factor  of  two,  and  the  anomaly  simplifies  to  ®  =  ®m  "  ®u 
(the  expression  (2.28)  is  altered  accordingly).  For  reasons  that  will  be 
explained  later,  the  layer  thickness  was  chosen  to  be  H  =  2500  m.  For  each 
station,  eM  was  identified  with  the  deep  relative  maximum  in  the  oxygen 
profile,  and  eu  the  value  1250  meters  above  this  (e^  so  defined  conforms 
approximately  to  an  isopycnal  surface).  With  this  definition  of  the  anomaly, 
the  denominator  in  (2.28)  was  estimated  for  the  region  within  the  inner 
streamline  of  Figure  2.20. 

To  compute  the  value  of  the  line  integral  in  (2.28)  we  broke  the 

streamline  into  two  parts,  a  northern  half  and  a  southern  half,  and  used  a 

single  representative  value  of  —  for  each  part,  estimated  from  the  average 

o  n 

section  in  Figure  2.21  (although  the  section  plots  values  of  eaw,  the  lat- 

a  V 

eral  gradients  of  closely  resemble  those  of  e).  Tnus  we  obtain  an 

estimate  for  the  term  in  brackets  in  (2.28)  (for  the  case  of  upward  flux 
only).  Since  there  exists  the  uncertainty  in  an  estimate  for  k ,  it  is 
best  to  keep  it  as  a  variable  and  think  in  terms  of  (<  ,  v)  pairs,  i.e.  for 
a  given  value  of  <  the  integrated  flux  balance  implies  an  associated  value 
of  v.  As  such,  (2.28)  predicts  the  following  pairs, 


2 

<  (cnr/sec) 

0(105) 

0(106) 

0(107) 


v  (cm^/sec) 

0(1) 

0(10)  (2.29) 

0(100). 


3.  Qy re  Level 

We  now  address  the  question  of  how  much  oxygen  is  contained  within  the 
NRG  with  regard  to  the  DWBC  source.  A  functional  relationship  between  the 
concentration  of  the  gyre  center  and  the  boundary  input  was  derived  from  the 
diffusive  transfer  model.  It  includes  the  values  of  the  various  time  scales 
but  also  relies  implicitly  on  the  fact  that  the  lateral  and  vertical  back¬ 
ground  reservoirs  are  void  of  tracer.  As  this  is  not  the  case  in  the  data, 
the  expression  (2.27)  must  be  revised  accordingly. 

We  will  assume  the  following  scenario  for  what  happens  in  the  ocean. 
Initially,  when  input  from  the  DWBC  begins  to  penetrate  the  interior,  there 
exists  a  non-zero,  uniform  background  oxygen  concentration  everywhere  below 
the  thermocline  resulting  from  remote  sources  and  sinks  (e.g.  diffusion  from 
the  surface  water).  Strong  lateral  diffusion  from  the  boundary  then  influ¬ 
ences  the  concentration  in  the  deep  layer,  but  vertical  diffusion  into  the 
layers  directly  above  and  below  is  not  substantial  enough  to  significantly 
raise  the  level  in  these  regions.  Thus  the  upper  and  lower  layers  remain  at 
the  concentration  of  the  initial  state.  This  scenario  is  extremely  oversim¬ 
plified,  but  such  a  representation  of  events  is  necessary  in  order  to  analyze 
the  data  within  the  context  of  the  transfer  model. 

To  incorporate  a  non-zero  background  into  the  transfer  model,  the  set 
(2.25)  must  be  altered  as  follows.  We  define  t  as  the  value  of  both  back¬ 
ground  reservoirs  as  well  as  the  initial  concentrations  of  and  de,  and 
write  the  equations  in  terms  of  the  rf-anomaly,  c ’  =  *ic  - 
$  '  =  <t>  -  i.  When  written  as  such  the  equations  for  the  primed  variables 
have  the  same  form  as  (2.25),  thus  we  use  those  solutions  derived  earlier 
with  the  appropriate  substitutions.  In  particular,  the  expression  for  the 
amount  of  tracer  at  the  gyre  center  (2.27)  becomes, 

*c(t  >  °°]  ~  f  7  +  (1  '  *  •  (2.30 

The  second  term  on  the  right  hand  side  represents  the  correction  due  to  non¬ 
zero  background. 

It  is  of  interest  to  contrast  (2.JUJ  with  the  analogous  expression  when 
vertical  mixing  is  absent.  In  that  case  the  amount  of  tracer  in  the  gyre  is 


Note  that  the  background  correction  for  the  two  cases  is  different.  Consider 
then  the  analog  to  the  function  f  when  there  is  non-zero  background  (recall 
that  f  measures  how  much  tracer  accumulates  in  the  gyre  relative  to  the 
amount  present  if  vertical  mixing  were  absent).  This  quantity  is  the  ratio 
of  the  expressions  (2.30)  and  (2.31).  When  written  in  terms  of  the  parameter 
*  =  -g  ,  winch  measures  the  strength  of  the  anomaly  supplied  by  the  DWBC,  it 
takes  the  form 


=  (»  -  l)f  ♦  2 
«  ♦  1 


(2.32) 


Figure  2.22,  which  plots  f  versus  <t>  ,  shows  that  for  $  >  1  (small  anom¬ 
aly)  f  >  1,  and  for  >  «  (large  anomaly)  f  >  f.  This  quantifies  the 
intuitive  result  that  when  there  is  non-zero  background  the  effect  of  ver¬ 
tical  mixing  is  lessened  (for  the  same  concentration  of  tracer  in  the  DWBC). 

In  the  data  we  can  readily  measure  the  equivalent  of  <5c(t  »  °°),  the 
concentration  at  the  center  of  the  gyre.  However,  in  order  to  apply  the 
result  (2.30)  we  have  to  identify  what  the  equivalent  boundary  conditions  are 
in  the  data,  i.e.  the  analogs  to  ♦  and  f i.  Consider  first  the  boundary  input 
value  *.  We  saw  from  applying  the  results  of  the  transfer  model  to  the 
numerical  model  that  in  terms  of  the  boundary  current,  the  appropriate  source 
concentration  that  dictates  the  level  in  the  gyre  is  the  value  of  tracer  at 
the  core  of  the  current.  It  is  a  matter  then  of  choosing  the  appropriate 
location  along  the  DWBC  for  which  this  relation  holds  in  the  data. 

If  this  location  is  chosen  so  far  upstream  of  the  gyre  that  some  of 
the  tracer  diffuses  directly  from  the  boundary  current  into  the  lateral  back¬ 
ground  (without  having  first  been  advected  around  the  gyre)  then  the  assump¬ 
tions  that  accompany  the  transfer  model  are  violated  and  the  relation  will  be 
inaccurate.  On  the  other  hand  if  the  location  is  not  far  enough  upstream 
then  some  of  the  tracer  that  enters  the  gyre,  having  diffused  from  further 
upstream,  will  not  be  accounted  for  and  again  the  relation  will  not  hold. 


Figure  2.22:  The  effect  of  a  non-zero  background  on  the  steady  state  gyre 
level  of  the  transfer  model.  The  level  of  the  gyre,  normalized  by  its  value 
when  vertical  mixing  is  absent,  is  shown  in  relation  to  the  strength  of  the 
input  anomaly. 
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The  correct  location  then  lies  just  in  between  these  limits,  and  in  light  of 
Figure  1.1  is  near  50°W.  Immediately  upstream  from  there  the  current  twists 
around  the  Grand  Banks  and  most  certainly  communicates  directly  with  the 
lateral  background.  (Note  that  in  the  numerical  model  the  northern  boundary 
could  be  anywhere  south  of  this  critical  location  and  the  relation  would 
still  hold.) 

In  regard  to  the  background  reservoir  value  t,  it  is  not  straightfor¬ 
ward  to  identify  such  a  quantity  in  the  data.  This  being  the  case  we  turn 
the  calculation  around  and  treat  i  as  the  unknown  (i.e.  solve  (2.30)  for 
d).  This  quantity  is  similar  to  a  quantity  such  as  the  diffusivity  in  that 
we  are  unable  to  directly  measure  it,  but  can  comment  on  whether  an  estimate 


is  reasonable  or  not.  As  a  guide  for  such  an  estimate  we  consider  the  fol¬ 
lowing  averaged  layer  value. 

Recall  that  in  the  description  of  the  numerical  model  governing  equa¬ 
tion  (2.2)  the  upper  and  lower  reservoir  values,  eu  and  eL,  were  indepen¬ 
dent  of  x  and  y.  In  order  that  the  model  behave  in  a  manner  similar  to  the 
ocean  (in  regard  to  vertical  mixing)  both  the  assumption  of  a  reservoir  and 
that  of  no  spatial  dependence  must  hold  to  some  extent  in  the  data.  In  view 
of  the  oxygen  profile  in  Figure  2.1  the  reservoir  assumption  for  the  upper 
layer  seems  valid  (keep  in  mind  that  the  shape  of  the  oxygen  profile  is  typ¬ 
ical  of  the  entire  data  set).  However,  the  lower  layer  value  never  levels 
off,  which  may  be  related  to  the  fact  that  this  value  is  consistently  higher 
than  its  upper  layer  counterpart.  The  idea  is  that  the  bottom  boundary  lim¬ 
its  the  size  of  the  lower  layer,  and  as  a  consequence  the  layer  begins  to 
fill  up  with  oxygen  that  has  diffused  from  above. 

In  order  to  address  the  question  of  spatial  dependence,  the  concentra¬ 
tion  of  the  02  maximum  was  compared  to  the  concentration  a  given  distance 
above  and  below  for  all  the  stations  within  the  inner  streamline.  It  was 
found  that  with  increased  separation  the  correlation  between  and  the  up¬ 

per  value  decreased  substantially,  whereas  the  correlation  between  eM  and 
the  lower  value  did  not.  This  is  consistent  with  the  reservoir  discrepancy 
mentioned  above,  and  suggests  that  when  modelling  the  effect  of  vertical 
mixing  we  should  in  fact  include  two  terms  in  the  governing  equation.  One 


should  be  a  decay  term  like  the  one  that  was  considered,  but  only  as  a  repre¬ 
sentation  of  mixing  with  water  above  the  deep  layer.  Mixing  with  the  bottom 
layer  is  characterized  more  appropriately  by  a  constant  flux  (independent  of 
the  concentration  of  the  deep  layer)  which  has  the  form  of  a  consumption 
term.  Recall,  however,  that  the  near  bottom  gradients  are  very  weak  and  the 
consumption  term  can  be  ignored  with  respect  to  the  radioactive  decay  term. 

As  a  value  then  with  which  to  compare  estimates  for  d,  the  average  of 
eu  at  a  distance  well  above  the  O2  maximum  was  computed  over  all  the 
stations  enclosed  by  the  inner  streamline.  The  distance  that  was  chosen  was 
H/2  =  1250  m,  as  a  distance  greater  than  this  does  not  cause  the  correlation 
to  decrease  significantly.  The  values  of  and  eu  for  the  group  of  sta¬ 
tions  appears  in  Figure  2.23  (also  shown  is  ej  .  The  two  quantities  exhibit 
similar  trends,  although  the  standard  deviation  of  eu  is  substantially 
smaller  than  that  for  eM.  It  should  be  noted  that  in  the  integrated  flux 

balance  calculation  above,  the  result  would  be  the  same  if  we  had  used  the 

value  of  (eM  -  ©u)  at  each  station  rather  than  (eM  -  eu). 

The  value  of  i  depends  on  the  value  of  f,  and  recall  that  the  func¬ 
tion  f  depends  on  the  parameters  a  and  e  (where  the  expression  for  e 

now  reflects  the  fact  that  there  is  flux  through  one  layer  only).  Estimating 
the  area  of  homogenization  from  the  data  set  and  using  a  value  of  H  =  2500  m, 
we  can  use  the  numerical  model  to  compute  an  associated  value  of  e  and  a 

for  each  pair  of  diffusivities  in  (2.29).  This  in  turn  gives  the  correspond¬ 
ing  magnitude  of  f.  It  is  the  case,  however,  that  the  values  of  d(f)  so 

predicted  from  (2.30)  are  significantly  larger  than  eu.  In  other  words  the 

transfer  model  says  that  there  should  be  less  oxygen  than  that  wnich  has 
diffused  into  the  NRG  according  to  the  data.  This  discrepancy  is  reconciled 
below  by  considering  variations  in  the  numerical  model  streamfunction. 

Western  Intensification 

Admittedly  there  is  some  question  as  to  the  validity  of  using  eu  as 
a  measure  of  the  "background"  oxygen  concentration  in  the  deep  layer.  This 
being  so,  an  inconsistency  concerning  the  level  of  oxygen  in  the  NRG  should 
not  in  itself  be  cause  to  completely  restructure  the  model  that  has  been  used 
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Figure  2.23:  The  value  of  the  oxygen  maximum  in  the  deep  layer,  and  the 
value  1250  m  above  this  for  each  of  the  stations  within  the  inner  streamline 
of  Figure  2.20. 
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In  addition,  there  is  uncertainty  in  the  integrated  flux  balance  calculation, 
and  it  is  possible  that  (2.29)  does  not  express  the  true  relationship  between 
the  two  diffusivities.  Nonetheless  it  is  desirable  to  try  and  fit  all  the 
pieces  of  the  calculation  together  in  a  consistent  fashion  and  offer  the 
result  as  one  possible  interpretation  of  the  data. 

The  velocity  field  of  the  numerical  model  was  chosen  as  a  simplified 
representation  of  the  DW6C/NRG  system,  and,  functioning  as  a  tool  to  reveal 
processes,  this  flow  field  in  conjunction  with  the  remainder  of  the  model 
proved  useful.  It  is,  however,  certainly  not  the  only  streamf unction  that 
could  have  been  used,  although  it  is  hoped  that  employing  a  similar  flow  pat¬ 
tern  would  not  alter  the  major  results  of  the  model  but  rather  have  only  a 
mild  effect  on  specific  features.  In  light  of  the  inconsistency  above,  one 
such  feature  that  we  concern  ourselves  with  now  is  the  strength  of  the  flux 
of  tracer  into  the  gyre. 

The  reason  why  we  choose  to  alter  the  velocity  field  is  twofold. 

First,  the  strength  of  this  flux  is  indeed  sensitive  to  the  specifics  of  the 
streamfunction.  Second,  the  current  meter  data  that  was  used  to  define  the 
NRG  is  somewhat  sparse,  thus  allowing  for  some  leeway  in  specifying  a  corres¬ 
ponding  simulated  flow  field,  i.e.  the  specific  details  are  really  an  unknown 
in  the  model  that  we  are  free  to  vary.  This  process  of  altering  the  stream- 
function  can  be  thought  of  as  fine  tuning. 

There  are  several  ways  that  the  velocity  structure  can  be  changed  to 
cause  tracer  to  fill  the  gyre  at  a  faster  rate.  Four  such  variations  appear 
in  Figure  2.24.  In  the  first  example  the  gyre  is  moved  closer  to  the  bound¬ 
ary  current  and  the  two  flows  partially  merge.  Tracer  now  diffuses  into 
stronger  flow  and  is  advected  more  quickly  into  the  interior.  In  the  second 
example  the  gyre  is  again  moved  towards  the  boundary  current  but  the  two  do 
not  merge,  rather  the  boundary  flow  converges.  This  strengthens  the  zonal 
gradient  of  tracer  there  causing  a  stronger  offshore  flux.  In  the  third 
example  the  gyre  is  turned  on  its  side  so  that  there  is  a  greater  region  of 
contact  between  the  flows,  which  also  serves  to  enhance  the  ventilation  of 
the  gyre.  In  the  final  example  the  gyre  is  no  longer  symmetric  but  western 
intensified.  This  causes  it  to  fill  up  more  rapidly  for  the  same  reason  as 
i n  the  fi rst  example. 
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Figure  2.24:  Schematic  showing  four  examples  in  which  the  streamline  pattern 
of  Figure  1.2  is  altered  in  a  way  that  causes  the  gyre  to  accumulate  tracer 
more  rapidly.  The  bounding  streamlines  of  the  boundary  current  and  gyre  are 
shown  alongside  a  zonal  velocity  section  through  the  current  and  center  of 
the  gyre.  The  unaltered  streamline  pattern  appears  at  the  top.  (a)  The 
gyre  is  shifted  towards  the  boundary  current  and  partially  merges  with  it. 

(b)  The  gyre  is  shifted  towards  the  boundary  current  and  causes  the  boundary 
current  streamlines  to  converge,  (c)  The  gyre  is  rotated  by  9U°.  (d)  The 

gyre  is  made  western  intensifed. 


It  is  possible  that  to  one  extent  or  another  all  of  these  effects  are 
present  in  the  ocean  and  have  not  been  accounted  for  in  the  model.  Certainly 
the  NRG  and  l)WBC  are  adjacent  over  a  greater  distance  than  is  portrayed  in 
the  model.  We  have  chosen  here  to  consider  the  effect  of  western  intensifi¬ 
cation,  as  the  data  implies  that  there  is  very  little  separation  between  the 
two  flows  wnen  they  are  side  by  side.  The  streamfunction  that  was  used  is 
shown  in  Figure  2.2b  compared  to  tne  previous  streamfunction.  The  boundary 
current  is  identical  in  both  cases,  and  the  only  major  difference  between 
the  gyres  is  that  now  the  maximum  southward  flow  of  the  gyre  is  comparable  in 
magnitude  to  the  core  speed  of  the  boundary  current  and  is  closer  in  proxim¬ 
ity  to  the  boundary  current. 

A  single  experiment  {abbreviated  W^ )  was  done  in  which  >  ~  1  x 
10bcm^/sec  (lateral  diffusion  only).  The  resulting  tracer  field  is  shown 
in  Figure  2.26  and  should  be  compared  to  that  of  the  experiment  (Fig¬ 
ure  l.I7a).  As  a  result  of  the  increased  speed  with  which  tracer  is  trans¬ 
ported  into  the  interior,  the  center  of  the  gyre  is  filled  at  a  faster  rate 
than  in  ■  ^.  However,  this  does  not  mean  that  more  tracer  ultimately  dif¬ 
fuses  into  the  gyre,  for  in  Figure  2.27  it  is  seen  that  the  level  asymptotes 
to  the  same  level  as  it  did  in  *  ^  (and  which  means  only  tnat  the  spin 

up  time  is  shorter.  Since  we  need  more  tracer  in  the  gyre  in  order  to  recon¬ 
cile  the  discrepancy  witn  the  data  it  seems  then  that  altering  the  stream- 
function  in  this  way  does  not  prove  useful.  Keep  in  mind,  however,  that  for 
oxygen  there  is  vertical  flux  present.  This  causes  the  system  to  spin  up 
rapidly  which  in  turn  preserves  any  differences  in  the  transient  states  of 
the  two  cases,  in  particular  the  fact  that  early  on  there  is  more  tracer  in 
the  western  intensified  gyre. 

In  order  to  apply  this  result  to  the  data  it  is  first  necessary  to 
couple  the  transfer  model  to  (note  that  we  do  not  have  to  alter  the 

set  up  of  the  transfer  model  to  reflect  the  change  in  the  numerical  model, 
because  the  specifics  of  the  advective  processes  are  not  included  within  the 
realm  of  the  transfer  model).  This  coupling  was  done  by  matching  (2.19a)  to 
the  spin  up  of  and  determining  a.  It  is  the  case  that  the  value  of 

is  again  small  as  it  was  in  the  •  ^  and  •  (j  experiments.  Tnis  tells  us 
immediately  then  that  is  in  the  same  limit,  namely  that  the  spin  up  is 
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Figure  2.2b:  Steady  state  distribution  of  tracer  for 
with  the  western  intensified  gyre  of  Figure  2.2S. 
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Figure  Z.Z7-  Time  history  of  tracer  accumulating  in  the  western  intensified 

gyre  (<  ~  I0°cm^/sec,  no  vertical  mixing)  compared  with  the  same  curve  for 
the  symmetric  gyre. 
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controlled  by  diffusion  from  the  boundary  and  behaves  as  a  single  exponen¬ 
tial.  Figure  2.28  shows  that  all  the  experiments  with  just  lateral  diffusion 
collapse  to  a  single  case. 

To  account  for  the  effect  of  vertical  flux  on  the  gyre  level  within  the 
context  of  western  intensification,  a  corresponding  value  of  e  must  accompany 
a.  Values  of  both  parameters  were  computed  for  the  diffusivities  and  length 

scales  relevant  to  the  data.  Recall  that  a  =  7—  and  e  =  —  .  We  need  the 

B  Tv 

numerical  model  to  compute  fg,  the  spin  up  time  of  the  system  (t  ~  i^) , 
as  this  quantity  can  not  be  deduced  from  the  data.  It  is  here  that  the  dis¬ 
tinction  appears  between  the  western  intensification  case  and  the  symmetric 
gyre  case.  The  time  scale  Tg  is  the  only  one  of  the  three  time  scales  that 
is  different  for  the  two  cases  (with  regard  to  the  data)  so  while  the  values 
of  e  that  accompany  (2.29)  will  be  the  same  here  as  before,  the  values  of 
a  will  be  somewhat  larger.  This  in  turn  alters  the  associated  values  of  the 
function  f  used  in  (2.30). 

As  seen  in  Figure  2.29  we  can  now  get  good  agreement  between  eu  and 

one  of  the  predicted  values  of  d(f).  The  pair  of  diffusivities  that  give 

6  2  2 

this  agreement  are  *■  ~  10  cm  /sec  and  v  ~  10  cm  /sec.  (In  order  to  get  a 

6  7 

prediction  for  <  ~  S  x  1U  cm  /sec  we  assumed  that  a  corresponding  western 

intensification  numerical  run  would  be  related  to  in  the  same  manner 

that  W^  was  related  to  <^).  Also  shown  in  the  figure  are  the  resulting 

predictions  if  we  disregard  the  flux  balance  constraint  and  arbitrarily  set 

v  *  1.  In  that  case  the  vertical  flux  is  extremely  weak  which  results  in  a 

comparatively  small  background.  When  the  gyre  of  Figure  2.27a  is  rotated 

90°  to  the  left  (which  combines  the  effect  of  increased  contact  with  western 

intensification)  a  spin  up  test  shows  that  tracer  fills  the  gyre  a  bit  more 

quickly  yet.  This  would  cause  even  closer  agreement  in  Figure  2.29. 

A  value  of  ~10  cm  / sec  for  v  is  substantially  larger  than  some 

earlier  estimates  of  the  vertical  diffusivity  at  thenmocline  depths  (for 

2 

instance  Rooth  and  Ostlund,  1972).  However,  values  as  large  as  3-4  cm  /sec 

have  been  calculated  for  the  Antarctic  bottom  water  flowing  northward  along 

the  western  boundary  of  the  South  Atlantic  (Hogg  et  al.,  ly82;  Whitehead  and 

6  7 

Worthington,  1982).  The  calculated  value  of  <  ~  10  cm  /sec  on  the  other 
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Figure  2.28:  The  accumulation  curves  for  each  of  the  three  experiments 
without  vertical  mixing  (Figures  2.5  and  2.27)  where  the  time  axis  has  been 
normalized  by  the  associated  spin-up  time. 
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Figure  2.29:  The  value  of  the  oxygen  background  concentration  as  predicted 
by  the  transfer  model.  The  particular  values  for  the  examples  discussed  are 
shown  in  comparison  to  the  measured  quantity  eu. 
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hand  is  smaller  than  most  of  the  estimates  found  in  the  literature.  For 

example,  direct  estimates  of  <  can  be  made  using  SOFAR  float  data,  and 

6  2 

Freeland  et  al .  ( 197 5)  calculated  a  value  of  7  x  10  cm  /sec  using  data  col¬ 
lected  for  the  MODE  area  (28°N,  69°W).  Using  the  same  technique  Price  (per- 

7  2 

sonal  communication)  computed  values  in  the  range  of  1.5-8. 2  x  10  cm  /sec 
for  the  LDE  region  (31°N,  70°W).  It  should  be  noted  though  that  these  esti¬ 
mates  are  for  the  thermocline,  and  Price  detects  a  decrease  in  the  size  of  < 
with  depth.  Still,  the  eddy  coefficients  calculated  here  are  for  the  area 
near  the  Gulf  Stream  which  is  an  area  of  increased  eddy  activity. 

B)  Conclusions 

The  review  of  the  data  in  Hogg  et  al .  (1986)  showed  that  in  the  deep 
layer  there  is  a  broad  region  characterized  by  very  weak  property  gradients 
to  the  south  of  the  DWBC  centered  at  roughly  60°W.  Such  homogenization  is 
suggestive  of  closed  circulation,  i.e.  the  Northern  Recirculation  Gyre.  In 
addition,  it  implies  that  the  Peclet  number  of  the  flow  is  large  in  the  sense 
discussed  in  chapter  one.  Closer  inspection  of  the  data,  within  the  context 
of  the  present  model  study,  yields  a  consistent  picture  in  regards  to  this 
feature  and  other  features  which  are  less  revealing. 

The  plume  of  oxygen  which  extends  from  the  DWBC  into  the  interior  de¬ 
limits  the  southern  extent  of  the  gyre  flow  that  has  turned  offshore.  This 
serves  as  additional  evidence  for  the  argument  of  tight  recirculation  of  the 
NRG.  In  contrast  to  the  salinity,  the  distribution  of  oxygen  has  a  slight 
relative  minimum  in  the  center  of  the  gyre  which  may  be  a  trademark  of  ver¬ 
tical  mixing.  The  Inward  flux  of  oxygen  associated  with  this,  along  with 
the  value  of  the  minimum  gives  rise  to  oceanic  mixing  coefficients  of  K  ~ 
10^cm2/sec,  v  ~  10  cm2/sec  for  this  region. 

The  value  of  the  vertical  diffusivity  seems  somewhat  large,  but  when 
viewed  in  terms  of  the  model  the  system  is  characterized  by  only  moderate 
vertical  mixing.  The  model  implies  that  the  amount  of  tracer  which  enters 
the  gyre  is  related  to  the  amplitude  of  an  appropriate  section  across  the 
DWBC,  most  likely  near  50  °W.  When  applied  to  the  data  this  says  that  the 
gyre  contains  approximately  40  percent  of  the  oxygen  that  is  available  to  it 


(that  is  to  say  40  percent  of  the  level  that  would  exist  if  there  were  no 
vertical  mixing  at  all).  In  light  of  how  sensitive  the  gyre  level  is  to 
vertical  mixing  (Figure  2.14)  it  is  the  case  then  that  this  effect  is  mild. 

This  should  in  turn  be  contrasted  to  the  case  of  the  symmetric  model 
gyre.  When  applying  that  case  to  the  data  it  turns  out  that  only  10  percent 
of  the  available  oxygen  makes  it  into  the  gyre  (recall  that  this  resulted  in 
an  inconsistency)  which  is  approaching  the  limit  of  strong  vertical  mixing. 

We  see  then  that  it  is  a  very  important  consequence  that  the  gyre  flow  come 
in  such  close  contact  with  the  boundary  current.  It  greatly  increases  the 
efficiency  with  which  tracer  is  transported  into  the  interior  and  enables  the 
lateral  flux  to  overcome  the  large  vertical  flux  in  supplying  the  NRG  with  a 
significant  amount  of  oxygen.  This  effect  has  no  bearing  on  the  salinity 
distribution  however.  The  results  of  the  numerical  model  imply  that  without 
vertical  flux  certain  characteristics  of  the  flow  field  have  no  effect  on  the 
level  of  the  gyre. 


CHAPTER  THREE:  TIME  DEPENDENT  INPUT 


Prel i mi naries 

In  the  previous  two  chapters  a  simple  model  was  presented  of  the  advec- 
tion  and  diffusion  of  a  passive  tracer  from  the  Deep  Western  Boundary  Current 
into  the  Northern  Recirculation  Gyre  of  the  deep  Gulf  Stream.  In  chapter  two 
when  the  results  were  applied  to  the  data,  only  the  distributions  of  salt  and 
oxygen  were  considered.  It  is  assumed  that  these  tracers  are  in  quasi-steady 
state,  that  is  to  say  the  source  of  tracer  provided  by  the  DWBC  is  not  chang¬ 
ing  substantially  in  time.  In  this  chapter  we  consider  what  happens  when 
this  source  grows  in  time,  as  is  the  case  with  freon. 

The  DWBC  fluxes  various  tracers  into  the  North  Atlantic,  among  which 
are  substances  called  chlorofluoromethanes.  The  two  which  are  most  readily 
measured  are  commonly  referred  to  as  freon-11  (F-ll)  and  freon-12  (F-12> . 
Freon  is  man-made  (for  instance  it  is  used  in  ref rigerants)  and  substantial 
amounts  began  to  accumulate  in  the  atmosphere  around  1950.  Since  that  time 
production  has  increased  rapidly.  As  is  the  case  with  oxygen,  the  surface 
waters  of  the  Norwegian-Greenl and  (N-Gj  Sea  that  sink  to  form  the  DWBC,  hav¬ 
ing  been  in  contact  with  the  atmosphere,  are  high  in  freon  content.  In  stark 
contrast  to  oxygen  though,  and  what  makes  freon  so  interesting  to  study,  is 
the  fact  that  they  have  been  present  for  only  a  short  time. 

By  keeping  track  of  the  production  of  freon  by  industry  versus  natural 
depletion,  the  atmospheric  source  function  from  1950  to  the  present  is  known 
with  some  certainty  (Smethie,  personal  communication).  This  presents  a  very 
good  opportunity  to  study  the  invasion  of  a  passive  tracer  into  the  world 
oceans  as  related  to  its  input.  This  is  not  the  only  appealing  aspect  of 
studying  freon  though.  Unlike  its  counterpart  tritium,  with  the  advent  of 
recent  technology  freon  is  relatively  easy  to  measure  and  can  be  done  so 
aboard  ship.  In  addition,  measurements  tend  to  be  characterized  by  a  high 
signal  to  noise  ratio  (Hogg  et  al . ,  198b). 

Within  the  realm  of  the  present  study  it  is  of  interest  to  consider 
freon  in  addition  to  the  steady  tracers  already  analyzed.  (For  instance  we 
know  how  long  freon  has  been  injected  at  high  northern  latitudes:  has  this 


been  long  enough  for  a  substantial  amount  to  accumulate  in  the  NRG?)  It  is 
however  not  trivial  to  incorporate  such  time  dependency  into  the  machinery 
that  has  been  developed.  This  is  because  there  is  much  involved  in  the  hi  s- 
tory  of  a  water  parcel  as  it  passes  from  the  surface  water  in  the  N-G  Sea 
via  the  DWBC  to  50°W,  at  which  point  it  enters  our  domain  of  study.  For  the 
cases  of  salinity  and  oxygen  we  could  ignore  this  previous  history,  it  was 
assumed  that  a  property  section  across  the  DWBC  remains  stationary  in  time. 
For  freon  this  is  obviously  not  true. 

In  the  numerical  runs  of  the  previous  chapters  the  input  of  tracer  into 
the  domain  was  specified  as  a  Gaussian  distribution  across  the  upstream  edge 
of  the  boundary  current  (a  step  function  in  time).  In  light  of  all  the  fac¬ 
tors  involved  in  the  overflow  process  as  well  as  the  shear  structure  of  the 
DWBC,  it  would  be  extremely  difficult  to  estimate  how  this  input  would  behave 
as  a  function  of  space  and  time  for  the  freon  case.  A:  such,  a  numerical 
model  with  a  regional  domain  is  not  useful  here.  However,  recall  that  the 
input  to  the  diffusive  transfer  model  corresponds  to  just  the  core  value  of 
the  distribution  across  the  boundary  current.  Therefore  to  make  use  of  the 
transfer  model  we  need  only  know  how  the  freon  core  value  evolves  in  time  at 
50°W  --  a  problem  which  is  much  more  tractable. 

It  is  the  case  then  that  before  we  make  use  of  the  transfer  model  for 
the  DWBC /NRG  system,  we  must  first  model  in  a  consistent  fashion  how  the 
atmospheric  input  of  freon  is  modified  during  the  processes  which  eventually 
lead  it  around  the  Grand  Banks  in  the  core  of  the  DWBC.  This  has  been  done 
in  two  stages:  the  deep  water  formation  stage  which  simulates  the  overflow 
process,  and  the  subsequent  boundary  current  stage  where  the  water  travels 
away  from  the  overflow  region  while  mixing  with  surrounding  water.  Both  of 
these  representations  are  greatly  simplified,  but  this  is  in  the  spirit  of 
the  transfer  model.  The  end  result  of  these  processes  is  finally  coupled  to 
the  transfer  model  and  contrasted  to  the  salt  and  oxygen  cases. 

In  the  literature  there  have  been  at  least  two  cases  in  which  the  age 
of  the  water  comprising  a  current  has  been  estimated  using  its  F-11:F-12 
ratio  values  {Smethie  and  Trumbore,  1984;  Weiss  et  al . ,  1985).  The  method 
assumes  that  the  flow  mixes  with  surrounding  water  which  is  mostly  void  of 
freon,  and  also  that  the  ratio  of  the  source  waters  is  not  altered  during  the 
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formation  process.  When  the  method  is  applied  to  the  DWBC  with  the  ABCE  data 
set,  the  age  of  the  water  so  calculated  suggests  an  unreasonably  slow  core 
speed.  Two  separate  boundary  current  models  are  presented  here,  each  based 
on  the  premise  that  as  time  progresses  the  mixing  that  occurs  is  no  longer 
with  freon-free  water.  When  considered  along  with  the  overflow  process,  this 
leads  to  substantially  larger  values  of  the  core  speed. 

Overflow  Basin  Model 

There  are  two  areas  at  high  latitude  where  convective  overturning 
occurs  supplying  deep  water  into  the  northern  North  Atlantic:  the  Norwegian- 
Greenland  Sea  and  the  Labrador  Sea.  The  water  which  is  formed  in  the  N-G  Sea 
is  more  dense  than  its  Labrador  counterpart,  and  it  is  this  water  which  upon 
entrainment  forms  the  DWBC  (Worthington,  1970).  There  are  three  sills  over 
which  the  newly  formed  water  in  the  N-G  Sea  flows  into  the  Atlantic.  These 
areas  are  the  Denmark  Straits,  the  Iceland-Faeroe  Ridge,  and  the  Faeroe  Bank 
Channel  (Figure  3.1).  The  latter  two  are  so  close  to  each  other  that  they 
are  usually  considered  as  a  single  overflow  region  called  the  Iceland- 
Scotland  Ridge. 

At  present  it  is  uncertain  as  to  what  extent  each  of  these  overflows 
contributes  to  the  water  which  comprises  the  DWBC  downstream  from  this  area. 
It  is  documented  that  the  Iceland-Scotland  overflow,  having  traveled  down 
the  slope  from  the  sill,  travels  westward  through  the  Gibbs  fracture  zone. 
According  to  Worthington  (1970)  this  flow  then  joins  the  Denmark  Straits 
overflow  and  the  combined  flow  eventually  makes  its  way  southward  as  the 
DWBC.  Swift  (1984)  argues  that  the  density  of  the  overflow  water  at  the  two 
locations  is  comparable,  but  the  Iceland-Scotland  flow  undergoes  more  intense 
mixing  as  it  progresses  from  the  sill  so  that  the  water  which  passes  through 
the  Gibbs  fracture  zone  is  less  dense  than  the  Denmark  Straits  contribution. 
As  a  result  the  two  flows  don't  really  combine  but  rather  influence  one 
another.  More  in  line  with  this,  Smethie  and  Trumbore  (1984)  present  a 
water  property  analysis  that  suggests  the  classical  DWBC  found  south  of  the 
Grand  Banks  is  comprised  mainly  of  Iceland-Scotland  overflow,  and  that  in 
this  area  the  flow  which  originated  from  the  Denmark  Straits  occurs  as 
weaker,  more  dense  filament-type  flows. 
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For  tni  s  work  it  has  been  assumed  that  tne  l)WBC  is  composed  solely  of  ; 

water  of  Iceland- Scotland  origin.  It  is  clear  though  that  there  are  other  j 

influences  at  work  with  regard  to  tne  freon  signal  of  the  current,  and  it 
would  oe  beneficial  to  study  this  point  in  greater  detail.  It  is  an  accepted 
idea  now  that  tne  fluid  which  overflows  the  Icel and-Scotl and  Kidge  (as  well 
as  that  which  overflows  the  Oenmark  Straits  sill)  originates  not  from  abyssal  1 

depths,  but  from  intermediate  depths  in  the  basin  (Bullister,  19B4) .  It  is 
consistent  then  that  the  water  at  these  depths  is  replenished  convectively 
[due  to  surface  cooling  (Worthington,  1970)  J  and  that  there  is  only  diffusive 
interaction  with  the  deepest  water. 

This  suggests  the  following  overflow  basin  model  to  compute  the  flux  of 
freon  over  the  sill.  An  influx  of  water  into  the  basin  occurs  in  the  surface 
layer,  the  water  acquiring  a  freon  concentration  according  to  the  atmospheric 
level  at  the  time.  As  a  result  of  heat  loss  to  the  atmosphere  the  fluid  then 
sinks,  ventilates  the  intermediate  layer,  and  resides  in  the  basin  for  a 
while  before  eventually  overflowing.  While  below  the  surface  the  fluid  can 
mix  diffusively  with  both  the  water  within  the  layer  and  with  the  abyssal 
water.  Based  on  equilibration  times  for  freon  it  is  reasonable  to  assume 
that  the  surface  layer  is  in  equilibrium  with  the  atmosphere  at  all  times 
(Smethie,  personal  communication).  Using  estimates  of  solubility  then,  the 
freon  concentrati on  of  the  surface  water  becomes  a  known  function  of  time. 

Tnis  information  is  the  starting  point  of  the  model,  it  is  from  thi s  function 
tnat  the  overflow  concentration  is  computed  and  ultimately  the  DWBC  core  con¬ 
centration  at  50 °W. 

The  details  of  the  model  are  outlined  in  Figure  3.2  which  shows  the 
simplified  scenario  of  ttie  overflow  process.  In  the  ocean  the  convective 
process  is  believed  to  happen  as  localized  events.  Because  such  events  are 
free  to  occur  anywhere  throughout  the  basin  we  envision  a  spatially  uniform 
flux  from  surface  layer  to  intermediate  layer,  which  is  represented  by  a 
specified  volume  transport  of  freon  concentration  e-(t).  Lateral  mixing 
also  occurs  uniformly  throughout  the  basin  stirring  up  tracer  in  the  inter¬ 
mediate  layer.  It  is  therefore  assumed  that  the  amount  of  freon  convected 
from  above  in  a  given  time  increment  is  "instantaneously"  diluted  to  a  uni¬ 
form  concentration  throughout  the  intermediate  layer.  This  representative 
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concentration  eQ(t)  is  then  also  the  concentration  of  the  overflow,  whose 
transport  matches  the  convective  transport.  To  keep  the  situation  as  simple 
as  possible  we  have  chosen  to  ignore  any  diffusion  between  the  intermediate 
layer  and  the  abyssal  layer,  arguing  that  the  effect  of  such  diffusion  is  of 
secondary  importance  relative  to  the  convection. 

The  freon  budget  for  the  intermediate  layer  that  results  from  this 
convective  source  and  overflow  sink  is  represented  as  follows. 


dt 


U(«i ( t)  -  eQ(  ti  i 


where  F  =  total  amount  of  freon  in  the  intermediate  layer, 

Q  =  transport  into/out  of  the  intermediate  layer, 
e^it)  =  surface  layer  concentration, 

and  eQ(t)  =  intermediate  layer  concentration  (overflow  concentration). 
The  quantity  of  freon  F  is  equal  to  VeQ ( t )  where  V  =  volume  of  the 
intermediate  layer,  and  constraining  V  to  remain  constant  the  above 
expression  can  be  written 


TR  \U>  *  “o11'  ■ 


(3.1) 


where  TR  =  residence  time  of  the  intermediate  layer. 

Assuming  that  the  intermediate  layer  is  initially  freon-free,  the  solution 
to  (3.1)  is 
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The  two  limits  of  (3.2!)  are  as  follows.  If  the  residence  time  of  the  layer 
is  inucti  smaller  than  the  time  scale  over  which  the  input  changes,  then  the 
overflow  concentration  mimics  the  surface  water  concentration,  or  in  terms  of 
(3.2), 


e0(t) 


(t)  . 


(3.3) 
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In  tne  opposite  extreme  when  the  input  time  scale  is  substantially  shorter 
than  the  residence  time,  (3.2)  gives 

1  r  t 

8  (t)  y  e.  (t')dt'  ,  (3.4) 

0  R  J0 

which  says  that  successive  input  concentrations  are  continuously  being  mixed 
in  the  intermediate  layer.  In  the  former  limit  note  then  that  there  is  no 
need  for  the  overflow  basin  model  in  so  far  as  determining  the  core  concen¬ 
tration  required  at  50°W. 

Figure  3.3  shows  the  surface  layer  F - 1 2  concentration  versus  time  from 
1950  to  1983  using  solubility  values  for  0“C,  3b°/°°  (1983  is  the  year  the 
freon  data  was  collected  in  the  region  of  the  NRu).  The  sill  depth  of  the 
Faeroe-Bank  Channel  is  ~850  m.  We  assume  that  the  surface  layer  is  approxi¬ 
mately  300  m  thick,  and  that  the  intermediate  layer  extends  from  tnere  to 
~1000  m.  Estimates  for  the  overflow  of  the  Icel and-Scotl and  System  put  the 
transport  at  roughly  1-2  Sverdrups  (for  example,  Worthi ngton ,  1970).  Approx¬ 
imating  tne  lateral  extent  of  the  region  which  supplies  water  to  this  over¬ 
flow  then  gives  an  estimate  of  the  residence  time  of  the  intermediate  layer 
Th  ~  10  years.  This  value  is  in  between  the  limits  (3.3)  and  (3.4),  and 
the  resulting  F-12  overflow  concentration  versus  time  is  shown  in  Figure  3.3 
compared  with  the  surface  water  concentration. 

There  has  been  only  one  cruise  which  has  sampled  freon  near  the  Faeroe- 
Bank  Channel.  This  was  the  HUDSON  82-001  expedition  conducted  in  1982. 

Bui  lister  (1984)  graphs  the  values  of  freon  versus  density  for  the  two  sta¬ 
tions  nearest  the  sill.  Using  the  estimate  of  the  overflow  density  there 
from  Swift  (1984)  this  gives  a  F - 12  overflow  concentration  of  ~1  p-mole/kg  in 
1982,  which  is  reasonably  close  to  the  basin  model  value  from  Figure  3.3. 
(Note  that  the  agreement  could  be  made  perfect  by  adding  a  small  amount  of 
mixing  with  the  abyssal  layer.) 

The  F-ll :F-12  Ratio 

One  aspect  of  freon  that  has  a  useful  application  is  the  fact  that  tne 
amount  of  F-ll  in  the  atmosphere  (and  surface  water)  has  been  increasing  at 


-  /  p-moles  \ 

Freon-12  — -  x  10 


Figure  3.3:  Concentrations  of  F - 1 2  in  the  surface  water  and  overflow  water, 
assuming  a  residence  time  of  10  years  for  the  intermediate  layer  of  the  over 
flow  basin.  The  first  measurement  (year  one)  corresponds  to  1950,  and  1983 
is  the  year  the  ABCE  freon  data  was  collected. 
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a  quicker  rate  than  that  of  F - 1 2  over  much  of  the  time  that  tne  two  have  been 
present.  Figure  3.4  plots  the  ratio  of  the  concentration  of  F-ll  to  that  of 
F - 1 2  in  the  surface  water  versus  time.  Suppose  the  water  which  forms  a  cur¬ 
rent  exhioits  this  trend  and  as  it  progresses  downstream  mixes  entirely  with 
fluid  that  is  freon-free.  This  means  that  while  the  F-ll  and  F - 1 2  content  of 
the  water  decreases,  the  relative  amount  remains  the  same  as  it  was  when  the 
water  left  its  source,  i.e.  the  value  of  the  ratio  remains  constant.  Since 
the  ratio  at  the  source  is  constantly  increasing,  this  means  the  water  down¬ 
stream  can  be  matched  to  the  time  when  it  was  isolated  from  the  source,  so 
indicating  the  age  of  the  water  and  the  speed  of  the  current. 

This  technique  has  been  used  by  Weiss  et  al .  (1985)  to  determine  the 
age  of  high  freon  intermediate  water  along  the  western  boundary  of  the  South 
Atlantic,  and  by  Smethie  and  Trumbore  (1984)  applied  to  the  UW8C.  In  both 
instances  input  functions  were  derived  from  the  atmospheric  concentration 
data  using  solubilities  appropriate  to  tne  source  waters  (the  Labrador  Sea 
and  Worwegian-breenland  Sea  respectively).  The  corresponding  ratio  curves 
were  then  used  to  determine  the  age.  Wote  that  this  corresponds  precisely  to 
using  ( t)  as  the  input  for  the  0W8C,  i.e.  ignoring  the  effect  of  the  for¬ 
mation  process.  In  Figure  3.4  it  is  seen  that  F-ll  :F-12  ratio  of  ( t ) 
flattens  out  as  of  1978.  The  above  ratio  curves  exhibit  this  feature  accord¬ 
ingly,  which  means  that  a  measured  value  of  the  ratio  near  that  level  can  at 
best  determine  the  age  within  a  range  of  possible  values.  This  difficulty 
arose  in  both  studies. 

In  the  overflow  basin  model  we  can  easily  see  how  the  F-11.-F-12  ratio 
is  affected  by  the  overflow  process.  In  the  small  residence  time  limit  the 
ratio  of  the  overflowing  water  is  identical  to  that  of  the  surface  water  (as 
are  the  concentrations  themselves).  In  the  large  residence  time  limit,  (3.4) 
gives 
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where  RQ(t)  =  ratio  of  overflow  water, 

e.  =  surface  water  concentration  of  F-ll, 

e,  =  surface  water  concentration  of  F-12. 

Mote  that  in  both  limits  the  ratio  is  independent  of  the  residence  time  T^. 
In  one  case  it  is  the  ratio  of  the  concentrations,  and  in  the  other  it  is  the 
ratio  of  their  integrals.  These  two  curves  are  plotted  in  Figure  3.5.  The 
difference  between  them  represents  the  maximum  extent  by  which  the  F-11:F-12 
ratio  can  be  altered  in  the  overflow  basin. 


In  between  these  limits  the  ratio  curve  does  depend  on  T^.  The  case 
for  TR  =  10  yr  is  also  plotted  in  Figure  3.5.  Note  that  because  of  the 
averaging  process  that  occurs  in  the  basin  the  10  year  overflow  curve  mono- 
tonically  increases  in  time  in  contrast  to  that  for  the  surface  water.  This 
means  then  not  only  are  the  ratios  lowered,  but  the  ambiguity  in  determining 
the  age  no  longer  exists.  These  results  suggest  that  it  is  important  to 
consider  the  overflow  process  when  applying  the  freon  dating  technique. 

Boundary  Current  Models 

In  the  ocean  as  newly  formed  water  flows  over  a  sill  it  proceeds  down- 
slope  and  entrains  water,  forming  a  current  of  increased  strength.  Through¬ 
out  tne  length  of  the  current  it  is  constantly  mixing  with  surrounding  fluid. 
The  fluid  with  which  it  mixes  contains  varying  amounts  of  the  properties 
which  are  being  advected  by  the  flow,  and  for  the  case  of  freon  it  is  usually 
assumed  that  the  level  in  the  neighboring  fluid  is  negligible  compared  with 
that  in  the  current.  We  are  presently  considering  only  the  core  of  the  flow 
and  within  this  context  can  devise  a  simulation  of  the  mixing  process  that 
enables  us  to  quantitatively  discuss  the  ratio  dating  process  and  the  freon- 
free  mixing  assumption  that  goes  with  it. 

We  assume  that  the  current,  which  has  no  "thickness"  and  is  represented 
by  a  speed  U,  flows  adjacent  to  a  reservoir  which  is  characterized  by  a 
given  (constant)  concentration  of  freon  eB.  The  mixing  that  occurs  is 
represented  by  a  diffusive  time  scale  (as  in  chapter  two).  The  governing 
equation  for  the  evolution  of  freon  in  the  current  is 
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Figure  3.5:  Value  of  the  F-11:F-12  ratio  in  the  overflow  water  compared  to 
that  in  the  surface  water  of  the  basin.  The  two  overflow  curves  are  those 
for  a  10  year  residence  time  and  an  infinitely  long  residence  time.  The 
solid  line  marks  the  ratio  value  of  1.97,  that  value  at  50°W  in  the  core  of 
the  DWBC  as  extrapolated  upstream  using  the  ABCE  data  from  1983.  It  is  seen 
that  the  TR  =  10  year  curve  reaches  this  value  ~8. 5  years  prior  to  1983, 
while  the  surface  water  curve  reaches  it  ”13  years  before. 
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where  e(x,t)  =  freon  concentration  in  the  current 
and  t  =  diffusive  decay  time. 

A  simple  scale  analysis  (using  the  diffusivities  estimated  in  chapter 
two  and  typical  lengths  scales  of  the  DWBC  freon  signal)  suggests  that  ver¬ 
tical  mixing  in  the  current  is  as  important  as  cross-stream  lateral  mixing, 
and  that  alongstream  lateral  mixing  is  negligible  with  respect  to  these.  In 
line  with  this,  in  (3.6)  we  have  ignored  the  latter  and  have  represented  the 
other  two  by  a  single  decay  term.  The  quantity  e0  should  be  thought  of 
as  a  basin-wide  average  concentration  of  the  surrounding  fluid.  The  upstream 
input  of  freon  to  (3.6)  is  the  overflow  concentration  versus  time  that  was 
computed  from  the  overflow  basin  model,  and  the  initial  condition  is  that 
the  level  of  tracer  in  the  current  (and  overflow)  is  the  same  as  that  in  the 
reservoi r.  The  solution  to  (3.6)  is 

x  >  Ut 

(3.7) 

x  _<  Ut  . 

Tne  data  set  that  was  described  in  chapter  two  contains  freon  measure¬ 
ments  at  the  station  locations  of  the  OCEANUS  134  cruise  in  1983  (Figure 
2.19).  This  provides  a  value  of  the  core  concentration  at  50°W  with  which  to 
compare  e(x  =  5500  km,  t  =  34  yr),  where  x  =  5500  km  is  the  estimated  dis¬ 
tance  along  the  DWBC  from  the  Iceland- Scotland  overflow  to  50°W,  and  t  = 

34  yr  corresponds  to  1983.  A  value  of  eB  to  use  is  the  freon  equivalent 
of  the  oxygen  background  that  was  used  in  chapter  two  (Figure  2.27).  With 
these  we  use  (3.7)  to  determine  a  value  of  t  =  1.06  yr  and  obtain  how  the 
core  concentration  of  freon  changes  in  time  at  50 °W,  which  is  the  information 
required  for  the  NRG  transfer  model. 

Consider  the  expression  for  the  F-ll :F-12  ratio  at  50°W  as  computed 
from  (3.7), 


where  e.  ,  en  =  overflow,  background  concentrations  of  F-ll 
°1  al 

eQ  ,  0B  =  overflow,  background  concentrations  of  F-12 
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L  =  distance  from  overflow  to  50°W. 

The  term  in  square  brackets  in  (3.8)  represents  the  deviation  of  the  ratio 

from  that  at  the  overflow  (Besides  from  the  delay  in  time  due  to  advection). 

a(p)eB 

In  order  for  the  deviation  to  be  small,  the  quantity  —  ^ ^  must  be  small 

o' 

for  both  F-ll  and  F-12.  The  parameter  p  is  a  Peclet  number  which  measures 
the  advective  strength  of  the  current  versus  the  mixing  that  occurs,  and  for 
large  values  of  p  the  function  a(p)  is  small.  So  there  are  three  condi¬ 
tions  which  will  cause  the  ratio  of  the  water  to  remain  unchanged  through  the 
course  of  progressing  downstream:  very  strong  input,  very  weak  background, 
or  very  weak  mixing  between  the  flow  and  surrounding  fluid. 

This  implies  that  the  freon  overflow  concentrations  have  indeed  become 
large  enough  by  1983  that  the  deviation  due  to  non-zero  background  is  neglig¬ 
ible,  i.e.  the  freon  dating  technique  is  applicable  to  the  present  data  set. 
However,  when  the  measured  ratio  is  matched  to  the  overflow  ratio  it  gives  an 
advective  time  of  8.5  yrs,  which  translates  to  a  core  speed  of  approximately 


2  cm/sec.  This  value  seems  too  small.  Below  we  discuss  two  separate  bound¬ 
ary  current  models  that  treat  the  mixing  differently  and  result  in  larger 
estimated  core  speeds  for  the  same  overflow  ratio  curve  of  Figure  3.5. 

Using  tne  value  of  t  and  a  value  of  U  =  8  cm/sec,  the  Peclet  number 
is  found  to  be  0(1).  This  means  that  the  freon  overflow  concentrations  have 
become  large  enough  by  1983  that  the  deviation  due  to  nonzero  background  is 
negligible,  i.e.  the  freon-free  mixing  assumption  appears  valid  in  the  pre¬ 
sent  case.  If  we  compute  the  age  of  the  water  as  was  done  in  the  studies 
described  above,  *his  results  in  a  DWBC  core  speed  of  ~1  cm/sec  (Figure  3.5). 
If  we  include  the  overflow  process  but  still  make  the  freon-free  mixing 
assumption,  this  increases  the  prediction  to  “2  cm/sec.  These  values  seem 
too  small.  For  example,  to  get  a  transport  of  10  Sverdrups  (the  value  esti¬ 
mated  by  Worthington  (1970)  for  the  DWBC)  for  a  current  200  km  wide  and 
1000  m  deep,  the  average  flow  speed  must  be  5  cm/sec.  Below  we  discuss  two 
separate  boundary  current  models  which  contain  different  mixing  schemes  that 
result  in  larger  estimated  core  speeds  for  the  same  overflow  ratio  curve  of 
Figure  3.5. 

A)  Back-Mixing  Model 

Expression  (3.8)  tells  us  just  how  large  the  background  concentration 
of  freon  would  have  to  be  in  order  to  affect  the  ratio  of  the  current  at 
50°W.  When  applied  as  above,  this  of  course  assumes  that  the  background 
value  computed  from  the  OCEANUS  134  data  set  is  valid  at  all  times  and  at 
every  location  along  the  current,  which  is  certainly  not  the  case.  Even  so, 
the  assumption  of  freon-free  mixing  results  in  an  unrealistic  core  speed  and 
suggests  that  we  explore  the  possibility  that  this  is  not  the  case.  Two 
scenarios  are  considered  here,  the  first  of  which  is  based  on  the  following 
simple  idea.  We  assume  that  the  first  bit  of  water  containing  freon  which 
overflows  the  sill  does  mix  entirely  with  freon-free  water.  However,  the 
amount  of  freon  that  it  imparts  along  its  path  to  the  surrounding  fluid  is 
enough  to  significantly  raise  the  freon  level  there.  Thus  the  next  bit  of 
water  that  flows  through  mixes  not  with  freon-free  water,  but  with  water  that 
has  a  concentration  reminiscent  of  the  previous  flow,  i.e.  some  of  boundary 
current  water  that  diffuses  into  the  neighboring  fluid  then  proceeds  to 


back-mix  with  newer  boundary  current  water.  It  is  evident  then  why  this 
effect  can  reduce  the  ratio  in  the  current,  for  its  water  is  constantly 
mixing  with  older,  lower  ratio  water.  The  difference  between  this  model  and 
the  process  described  by  equation  (3.6)  is  that  here  we  keep  track  of  the 
surrounding  fluid  as  well  as  the  fluid  in  the  current. 

Consider  a  short  section  of  boundary  current  (and  surrounding  fluid) 
as  depicted  in  Figure  3.6,  extending  from  the  overflow  region.  We  identify 
three  regions:  the  inner  core  which  moves  with  constant  speed  Uc  (i.e.  the 
current),  the  adjacent  shoulder  water  (which  is  still),  and  the  "infinite" 
amount  of  fluid  which  surrounds  this.  It  was  mentioned  earlier  that  vertical 
mixing  and  lateral  cross-stream  mixing  are  of  comparable  strength  in  the  cur¬ 
rent  (and  that  alongstream  mixing  can  be  ignored).  We  are  now  more  precise 
about  this.  Consider  the  advective-diffusive  equation 

e  +  ue  =  +  <e  +  ve„,  (3.9) 

t  x  xx  yy  zz 

where  x  =  alongstream  distance, 

y  =  cross-stream  distance, 
z  =  vertical  distance, 

k,  v  =  lateral,  vertical  diffusivity  (assumed  constant), 
u  =  flow  speed, 

and  e  =  freon  concentration. 

We  let  H,  L  ,  and  L  represent  the  vertical,  cross-stream,  and  alongstream 
y  * 

length  scales  of  the  tracer  distribution.  Inspection  of  the  ABCE  data  shows 

that  in  the  boundary  current  H  <<  L  <<  L  .  If  (3.9)  is  scaled  accordingly 

y  * 

it  is  seen  that  the  ratio  of  the  alongstream  diffusive  term  to  the  cross- 

stream  diffusive  term  is  0 ( L  /L  )  <<  1,  and  so  to  lowest  order  the  balance 

y  * 

is 

9t  +  U9x  =  K0yy  +  v9zz  •  (3-10) 

We  now  introduce  a  stretched  coordinate  in  the  vertical  direction.  In 
particular,  we  define 


Figure  3.6:  Schematic  depicting  the  DWBC  as  being  composed  of  many  short 
sections  of  length  La.  An  enlargement  is  shown  of  the  first  such  section 
originating  from  the  overflow  basin.  The  flow  is  confined  to  the  innermo; 
region. 


The  flow  is  confined  to  the  innermost 


1  ’*  •  K.  *  -  ■  •. 


■  .''.-■■.Vv’vV>V/v‘vV'/vv/v\>Vv' 


in  terms  of  which  (3.10)  becomes 


9t  *  u8x  =  *8yy  * 


l/v 

tf"* 


Vz* 


(3.11) 


In  chapter  two  we  constrained  the  values  of  k  and  v  by  the  relation  (2.28), 

which  says  that  </v  ~  10^  (see  (2.29)).  Using  H  ~  1000  m,  L  ~  300  km 

2  2  5  “ 

(Figure  3.7)  this  results  in  /H  ~  10  as  well,  which  implies  that 
both  types  of  mixing  are  equally  as  important.  We  therefore  set  the  ratio 
< H^/vL  ^  =  1,  and  rewrite  (3.11)  as 

8t  *  u8x  =  ‘<8yy  *  8zV>  l3’12) 

Ly2 

where  <  =  (-*— )v  . 

H 


The  stretched  coordinate  then  increases  the  vertical  diffusive  length  to  tne 
size  of  the  cross-stream  diffusive  length,  so  that  diffusion  acts  isotropic¬ 
ally  in  the  y-z'  plane. 

It  is  this  coordinate  frame  in  which  the  freon  signal  associated  with 
the  boundary  current  appears  circular  (Figure  3.7)  and  for  which  the  schem¬ 
atic  in  Figure  3.b  applies.  Accordingly,  the  section  of  boundary  current  is 
described  using  cylindrical  coordinates  (x,  r,  x) .  The  x  direction  is  taken 
alongstream,  and  r  and  x,  the  radial  and  azimuthal  coordinates,  are  defined 
by 

y  =  r  cosx, 
z'  =  r  sinx. 

Equation  (3.12)  is  then  written 


8t  *  uex  -  7  ,r*r)l-  '  u-13’ 

It  is  assumed  that  the  concentration  of  tracer  is  independent  of  x.  (In  the 
outer  reservoir  the  concentration  is  taken  to  be  zero.)  The  alongstream  and 


radial  dependency  is  treated  using  a  finite-difference  representation  so  that 
the  concentration  in  each  region  is  characterized  by  a  single  value. 

The  length  of  the  section  La  is  taken  small  enough  that  fluid 

a 

flushes  through  the  core  region  very  quickly.  Correspondingly,  it  is  assumed 
that  the  region  adjusts  instantly  to  changes  in  input,  and  so  is  governed  by 
the  steady  state  balance 


Ucex  =  7  ^rerV  at  a11 


(3.14a) 


The  shoulder  water  acquires  tracer  by  way  of  diffusion  from  the  core,  and 
loses  tracer  through  mixing  with  the  fluid  in  the  outer  reservoir.  The 
shoulder  water  is  governed  by  the  equation 


9t  -7  <r»r>r 


(3.14b) 


Using  upstream  differencing  in  the  alongstream  direction,  (3.14a)  and 
(3.14b)  become 


ejr.t)  -  (j  i(r,t)  .  _ 

Uc  (-1 - [-—Li - ,  (r.tl) 

a 


ft9j(r-*’  -7T7  (r  jr  9j(r't,)  ■ 


(3.15a) 


(3.15b) 


where  the  subscript  j  denotes  the  x-grid  point  index  (L,  is  the  grid 

a 

spacing).  Integrating  (3.15a)  and  (3.15b)  over  the  cross-sectional  area  of 
the  core  and  shoulder  regions  respectively  gives 


drdxr(e  .(r,t) 

J 


-  9j_1(r,t)) 


<irr  —  e.(r  ,t) 
c  ar  j  c’ 


(3.16a) 


rs  rX2 


rc  A1 


drdxr Sj(r,t)  =  <1,(rs  7?  ej^rs,t:)  “  rc  7r  ej  (rc,t} }  *  (3.1bb) 


where  the  limits  of  integration  are  depicted  in  Figure  3.6. 
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Using  centered  differencing  in  the  radial  direction,  (3.16a)  and 
(3.16b)  become 


A  J 
c  c 


‘V1’ 


'j-1 


(t))  = 


Kirr 

C  J  J 


d 

■at 


(t) 


Kirr  icur 

-rjr1  (0  -  e  ( t) )  -  (e  (t)  -  e.  (t))  , 

°  5  j  Q  f  C  j  ^  j 


(3.17a) 


(3.17b) 


where  ar,  and  ar$  are  the  grid  spacings,  A  denotes  cross-sectional  area, 

and  the  subscripts  c  and  s  are  the  radial  grid  point  indices  referring  to 

the  core  and  shoulder  regions  respectively  (Figure  3.6). 

The  input  into  the  core  region,  er  ( t)  ,  is  set  equal  to  the  over- 

Cj-1 

flow  ©0(  t ) .  Also,  we  do  not  want  the  shoulder  region  to  be  too  vast  (to 
allow  it  to  readily  accumulate  tracer),  so  A$  is  set  equal  to  Ac.  (This 
serves  as  a  simplification  as  well.)  Then  with  some  rearrangement  (3.17a) 
and  (3.17b)  become 


e  (t)  =  0  (t)  -  (er  (t)  -  e,.  ( t) ) 
cj  0  c  cj  sj 


(3.18a) 


ft*.  •*>  ■ 

J 


(e  (t)  -  y9  (t)) 

J  SJ 


(3.18b) 


where  Tr  =  =  residence  time  of  the  fluid  in  the  core  region, 


Vrc  i 

t .  =  ( — n - )  -  =  diffusive  time  between  the  core  and  shoulder, 

c  pc  K 


Vrs  1 

t  =  (-p - )  —  =  diffusive  decay  time  of  the  shoulder  into  the  reservoir 

s  ■  5 


=  1  +  t  * 


and  p  ,  p  =  outer  perimeter  of  core,  shoulder  region. 


The  parameter  y  is  called  the  reservoir  parameter.  The  value  of  1/y 
represents  the  fraction  of  the  core  concentration  that  the  shoulder  water 
would  equilibrate  to  if  given  sufficient  time.  This  parameter  is  tne  key  to 
the  model.  When  >  0  (  y  >  00 )  the  outer  reservoir  prohiDits  any  tracer 
from  accumulating  in  the  shoulder  region,  so  the  core  mixes  with  freon-free 
water.  We  are  interested  then  in  finite  values  of  y.  When  x>  ~(y=  1) 
this  corresponds  to  the  case  when  the  shoulder  water  is  most  readily  filled 
with  tracer  (the  outer  reservoir  is  "turned  off").  Thus  the  larger  the  value 
of  y,  the  more  the  shoulder  water  acts  as  a  reservoir  sink  for  the  core. 
When  y  =  1  this  is  referred  to  as  the  finite  reservoir  case,  and  when  y  > 

00  this  is  called  the  infinite  reservoir  case.  In  between  these  limits,  the 
shoulder  water  can  be  thought  of  as  a  semi-infinite  reservoir. 

In  the  limit  when  the  mixing  between  the  core  and  shoulder  goes  to  zero 
(t  »  00 )  the  concentration  of  freon  leaving  the  core  is  the  same  as  that 
which  enters.  In  the  limit  of  very  strong  mixing  (t  »  0)  the  coupled  set 
(3.18)  gives 

e  ( t)  =  e  ( t) 

j  j 

e  (t)  =  y9  (t) 

j  j 

(remember  that  ec  (ti  and  es  (t)  vary  on  the  time  scale  of  the  input).  For 
j  j 

y  >  1  this  implies  that  e  (t)  =  e  (t)  =  0,  which  is  consistent  in  that 

j  j 

no  freon  can  accumulate  in  either  the  core  or  shoulder.  When  y  =  1  the 
shoulder  will  instantly  equilibrate  to  the  level  in  the  core,  and  the  govern¬ 
ing  equations  tell  us  only  that  e  (t)  =  e  (t).  It  is  obvious  that  we  do 

ct  sj 

not  wish  to  consider  either  of  these  cases,  thus  tc  is  constrained  to  be 

>  T  ( and  f i ni te) . 

—  r 

Because  the  fluid  spends  only  a  short  time  inside  the  core  and  because 
the  ratio  of  the  shoulder  continually  increases  as  well,  the  amount  by  which 
the  ratio  of  the  boundary  current  is  lowered  is  minimal.  However  if  this 
process  occurs  along  the  entire  length  of  the  current  the  effect  can  be  sub¬ 
stantial.  So  we  think  of  the  boundary  current  (and  shoulder  water)  then  as 


being  made  up  of  many  of  these  sections,  each  one  represented  by  a  different 
value  of  the  index  j.  The  output  from  the  core  of  a  given  section  is  in 
turn  the  input  to  that  of  the  adjacent  one.  However  we  do  not  allow  interac¬ 
tion  between  adjacent  shoulder  regions.  (Note  that  this  Is  consistent  in 
that  we  are  ignoring  alongstream  diffusion.)  It  is  evident  that  the  amount 
by  which  the  core  ratio  decreases  depends  on  the  length  of  time  the  fluid 
spends  in  each  section,  which  will  enable  us  to  track  the  ratio  versus  advec- 
tive  time  and  match  the  observed  ratio  in  the  data. 

The  solution  to  (3.18)  is 


9  At)  ft  r+. 

0(11  -ttt* - 5 - 7  ert  e  ( t'  )e_rt  dt' 

Cj  Tc(l  ♦  «)Z  0 


(3.19a) 


..  ft  e,  ( t* ) 

V(tl  *  e  t  |1  »  i)  e  dt  • 

J  Jc 


(3.19b) 


where  r  =  -  7-  (y  -  t  ~-.-) 
c  1  4 


and  6  =  —  . 

c 


These  represent  the  concentrations  for  the  first  section  (j=l),  and  (3.19a) 
should  in  turn  be  used  as  the  input  to  the  next  core.  Considering  successive 
sections  as  such,  however,  results  in  an  insurmountable  amount  of  algebra, 
and  a  simplifying  assumption  must  be  made.  What  we  do  is  approximate  the 
overflow  concentration  curve  of  Figure  3.3  (for  both  F-ll  and  F - 12 )  by  an 
exponential . 


When  the  expression  eQ(t)  =  Aet/T  is  substituted  into  (3.19),  this 


gives 


,  90(t) 

ec  (t)  =  + - 

Cj  1  *  TC<1  *  6 


le  It)  -  Aert) 


(3.20a) 


e,  (t)  = 


t(1  ♦  6 )  ( -  r  ♦▼) 


(e0(t)  -  V  )  * 


(3.20b) 
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The  important  thing  here  is  that  the  terms  proportional  to  e  decay  very 
quickly  relative  to  those  proportional  to  eQ(t)  (which  grow  in  time).  This 
means  that  to  a  good  approximation. 


•c.(t)  ~  Gc(tc,  Tr,  T,  Y)e0U) 
J 


e$  (t)  '  Gs(tc,  Tr,  T,  Y)eQ(t)  , 
3 


(3.21a) 

(3.21b) 


where 


GC  -T^T 


r  * 

T  (1  ♦  6 j  (-r  ♦  j) 


and 


G_  = 


tc(1  ♦  6)(-r  ♦  f) 


With  this  being  the  case  it  is  straightforward  to  combine  successive  sec¬ 
tions.  After  n  sections  (3.21)  gives  that  the  concentrations  are 


e_  (t)  ~  (G  )ne  (t) 
c j=n  c  0 


(3.22a) 


.  n-1 


e  (t)  ~  ( G  ) ( G  )  e  ( t) 
s j=n  s  c  ° 


(3.22b) 


The  functions  Gc  and  Gs 


are  called  decay  factors  (whose  magnitudes 

are  _<  1).  For  F-il  the  overflow  concentration  growth  scale  T  is  ~  5  yr, 

and  for  F-12,  T  ~  7  yr.  It  is  this  slight  difference  which  is  the  reason 

why  the  ratio  proceeds  to  decrease  as  the  fluid  progresses  downstream.  The 
fG(T=5)ln 

quantity  ^  (j  -  77  measures  this  decrease,  and  in  line  with  what  was 

mentioned  above,  only  for  a  large  enough  n  is  the  difference  substantial. 
Figure  3.8  plots  the  core  concentration  decay  factor  and  ratio  decay  factor 
versus  the  strength  of  the  mixing  between  core  and  shoulder  for  the  three 
different  types  of  reservoirs  (there  is  no  ratio  decay  for  the  infinite 
reservoir).  As  expected  the  concentration  decays  more  with  smaller  Tc> 
Note  the  difference  between  the  ratio  decay  of  the  finite  reservoir  case 


v.V 


finite  reservoir 
semi-infinite  reservoir 
infinite  reservoir 


Tc( yr) 


Figure  3.8:  The  dependence  of  the  core  decay  factor  on  the  strength  of  the 
mixing  between  the  core  and  shoulder  for  three  different  values  of  the  reser 
voir  parameter  y:  y  =  1  (finite  reservoir),  y  =  1.25  (semi-infinite 
reservoir),  T  >  »  (infinite  reservoir).  The  value  of  the  residence  time 
of  the  sections  is  taken  to  be  .2  years,  (a)  The  freon  concentration  decay 


factor 


u  (T  =  5  yr) 


♦  GC(T  =  7  yr) 


versus  the  semi-infinite  case.  In  the  latter,  not  only  is  the  decay  less 
pronounced,  but  for  small  values  of  tc  the  decay  decreases  with  increased 
mixing.  This  is  because  for  very  strong  mixing  no  tracer  is  able  to  accumu¬ 
late  within  the  shoulder  water  (in  contrast  to  a  finite  reservoir  where  the 
shoulder  level  approaches  that  of  the  current). 

As  seen  in  Figure  3.9,  for  a  given  value  of  mixing,  as  the  size  of  the 
reservoir  parameter  gets  smaller  the  amount  of  ratio  decay  increases,  whereas 
the  amount  of  concentration  decay  decreases.  So  depending  on  the  value  of  y 
which  is  chosen,  the  model  will  predict  a  different  core  speed.  Note  also 
that  within  the  context  of  this  model  the  freon  dating  technique  involves  not 
only  matching  the  ratio  but  the  concentration  as  well.  The  technique  is 
appl ied  as  follows. 

First  we  choose  the  magnitude  of  the  residence  time  that  characterizes 
the  component  sections.  The  value  that  was  used  is  Tp  =  .2  yr  (this  is 
arbitrary,  although  the  smaller  the  value  of  T  ,  the  higher  the  resolu¬ 
tion).  For  a  given  value  of  y  then  there  are  two  unknowns,  the  strength  of 
the  mixing  x  and  the  value  of  n.  It  is  then  just  a  matter  of  determining 
these  unknowns  by  forcing  the  core  concentration  of  F-12  and  the  value  of  the 
ratio  (or  equivalently  the  two  concentrations  themselves)  to  match  the  data 
at  50°W .  (The  ABCE  region  extends  only  to  55°W;  an  extrapolation  was  used  to 
obtain  a  core  value  of  3.34  p-moles/kg  x  10  for  F-12  and  1.97  for  the 
F - 11  :F12  ratio  at  50°W.)  It  should  be  noted  that  even  though  we  make  use  of 
(3.22),  the  actual  eQ( t )  curves  are  used  in  the  calculation  (rather  than 
the  exponential  approximations). 

For  the  case  of  a  finite  reservoir  (y  =  1)  it  is  not  possible  to  get 
a  consistent  solution.  This  is  because  when  the  ratio  is  matched  to  the 
value  in  the  data  there  is  too  much  tracer  in  the  shoulder  water  to  drive 
down  the  F-12  concentration.  In  the  infinite  reservoir  case,  Gc  =  y — 

1  n 

so  the  concentration  decay  amplitude  after  n  sections  is  (  r  ■+■  . )  .  From 

L 

'  UT 

(3.7)  we  see  that  in  the  continuous  case  the  equivalent  quantity  is  e 
(there  is  no  background  so  9g  =0).  In  the  limit  as  n  >  °°  the  finite- 
difference  solution  approaches  the  continuous  result.  This  can  be  seen  by 


y 


Figure  3.9:  The  dependence  of  the  core  concentration  decay  factor  and  core 
ratio  decay  factor  on  the  size  of  the  reservoir  parameter  y.  The  strength 
of  the  mixing  between  the  core  and  shoulder  is  held  constant  (tc  =  1.5  yr) 
as  is  the  residence  time  of  the  sections  (Tr  =  .1  yr). 
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noting  that  &  = 


and  by  definition  lim  - - r 

n»=°  l  ♦  (-L-)  - 
uxc  n 


(t  and  x  are  equivalent  in  this  case).  Thus  the  scenario  which 


accompanies  the  freon-free  mixing  assumption  is  contained  as  a  special  case 
within  the  present  boundary  current  model. 

In  Figure  3.10  the  resulting  core  speeds  are  plotted  versus  various 
values  of  y.  For  y  nearly  equal  to  one,  the  model  gives  a  core  speed  as 
large  as  6  cm/sec  (remember  this  is  an  average  over  the  length  of  current, 
from  the  overflow  to  50 °W) .  As  a  possible  guide  in  determining  which  of 
these  examples  might  be  the  most  realistic  in  terms  of  the  ocean,  we  can  con¬ 
sider  the  corresponding  values  of  t  and  e  (t)  for  each  example.  Smaller 

j 

values  of  the  reservoir  parameter  require  stronger  mixing  (in  each  example 


tc  >  T  )  which  accompany  higher  shoulder  water  concentrations  (Figure  3.11). 


An  observation  that  has  led  to  the  freon-free  mixing  assumption  is  that 
successive  locations  along  a  current  show  little  variation  in  the  core  value 
of  the  F-11:F-12  ratio  (Smethie  and  Trumbore,  1984).  Because  „he  ratio  is 
not  a  conserved  quantity  in  this  model  (for  finite  y)  it  might  seem  then 
that  the  resulting  downstream  variation  contradicts  this  observation.  This 
is  not  the  case  however.  In  fact,  even  though  the  amount  of  decay  varies 
with  y,  the  alongstream  change  in  ratio  is  comparable  in  each  case. 

The  OCEANUS  134  data  set  contains  several  crossings  of  the  UWBC  with 
which  to  compare  the  model  predictions  concerning  this  point.  In  the  limit 
of  freon-free  mixing  (y  >  °°)  the  core  ratio  varies  slightly  alongstream 
because  the  ratio  at  the  source  (overflow)  is  changing  in  time  and  the  core 
speed  is  only  2  cm/sec.  With  such  a  small  core  speed  the  input  at  a  given 
time  cannot  propagate  very  far  downstream  before  another  change  in  input 
occurs.  The  resulting  downstream  trend  in  ratio  agrees  well  with  what  is 
observed  in  the  data.  In  the  opposite  limit  of  y  >  1,  the  ratio  decays  as 
the  fluid  proceeds  downstream,  but  because  it  flows  so  quickly  a  signal  will 
propagate  a  long  way  before  the  input  changes  appreciably.  Thus  the  effect 
that  caused  the  variation  for  the  y  >  °°  case  has  little  consequence  here, 
and  it  turns  out  that  the  decay  process  results  in  an  alongstream  gradient  of 
nearly  identical  magnitude.  In  between  these  limits  the  importance  of  the 
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Figure  3.11:  Graph  showing  the  mixing  that  is  required  for  the  different  core 
speeds  and  the  resulting  cross-stream  variation  in  concentration  for  the  back- 
mixing  model.  The  solid  line  measures  the  strength  of  the  mixing  between  the 
core  and  shoulder  relative  to  what  it  would  be  if  the  core  were  mixing  with  a 
freon-free  shoulder.  The  dashed  line  plots  the  concentration  of  the  shoulder 
relative  to  that  of  the  core  (x  10). 


F-12  Differential  X10 


two  effects  varies  inversely  making  it  hard  to  distinguish  between  any  of  the 
cases.  The  observed  feature  in  the  data  can  therefore  be  explained  by  any  of 
the  examples  of  the  model  (as  can  the  variation  in  concentration,  Figure  3.12) 


B)  Shear  Model 

In  the  above  representation  it  was  assumed  that  the  boundary  current  had 
no  structure  to  it,  i.e.  the  flow  was  represented  by  a  single  core  vel¬ 
ocity.  It  was  seen  that  the  mixing  between  the  core  and  shoulder  water  re¬ 
sulted  in  decay  of  the  core  ratio  when  the  shoulder  did  not  act  as  an  infin¬ 
ite  reservoir  sink.  It  is  easy  to  envision  a  different  type  of  scenario  which 
could  also  decrease  the  freon  ratio  of  the  core.  We  consider  a  boundary  cur¬ 
rent  that  has  cross-stream  shear  (so  that  the  core  refers  to  only  the  fastest 
part  of  the  flow)  and  assume  that  the  core  mixes  entirely  with  the  outlying 
weaker  flow.  Because  this  slower  travelling  water  is  older  and  has  a  lower 
freon  ratio,  this  will  tend  to  drive  down  the  core  ratio.  We  structure  the 
model  the  same  way  as  the  previous  back-mixing  model  in  that  the  boundary 
current  is  represented  as  a  number  of  sections. 

Consider  the  first  such  section,  which  receives  its  input  from  the 
overflow  as  depicted  in  Figure  3.13.  As  before  it  is  comprised  of  the  core, 
the  shoulder,  and  the  vast  amount  of  surrounding  fluid  (which  again  is  taken 
to  be  an  infinite  sink).  The  difference  here  though  is  that  the  shoulder 
water  is  moving  as  well  --  the  core  is  characterized  by  speed  Uc  and  the 
shoulder  by  U$  (<  Uc).  So  whereas  previously  the  current  corresponded  to 
just  the  core  region,  now  it  is  made  up  of  the  core  and  shoulder  together. 

We  assume  again  that  the  concentrations  are  uniform  in  the  azimuthal  direc¬ 
tion  and  apply  the  same  finite-difference  approximations  radially  and  along- 
stream. 

In  the  core  region  the  governing  equation  is  as  before. 


U  e  =  —  ( re  ) 
c  x  r*  '  r  r 


(3.23) 


Because  of  the  advection  in  the  shoulder,  (3.2 3j  is  the  governing  equation 
there  as  well,  with  \i^  replacing  Uc>  After  integrating  over  the 
cross-section  and  finite-differencing,  the  two  governing  equations  become 
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Figure  3.12:  The  alongstream  variation  in  core  concentration  and  ratio  as  pre¬ 
dicted  from  the  back-mixing  model,  compared  with  that  observed  in  A BCE .  Three 
different  examples  are  shown.  The  curves  are  normalized  by  the  values  at  i>0  W. 


Figure  3.13:  Schematic  depicting  the  l)W BC  as  being  composed  of  many  short 
sections,  as  In  Figure  3.b.  In  the  shear  model  flow  occurs  in  both  the  core 
and  shoulder  regions. 
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U  ic*r  Kir. 

Ac  T5-  (ec  <*>  -  ec  (*)  =  T^°  -  ®c  (t))  -  — =■  (ec  (t)  -  e  (t))  (3.24b) 

S  La  sj  sj-l  Ars  sj  Arc  sj  cj 


which  should  be  compared  to  (3.17a,b).  [All  labels  here  are  the  same  as  in 
(3.17a, b) . J 


Both  e  (t)  and  ec  (t)  are  set  equal  to  the  overflow  e  ( t). 
cj-l  sj-l  0 


Also,  we  take  r  <<  r  so  that  the  last  term  on  the  right  hand  side  of 
c  s 


(3.24b)  can  be  ignored.  (Note  that  as  rc  >  0,  ar£  remains  finite.)  This 


says  that  the  area  of  contact  between  the  core  and  shoulder  is  so  small  that 
the  diffusive  flux  into  the  shoulder  across  this  surface  is  negligible 
compared  with  the  flux  across  the  outer  surface  of  the  shoulder.  The  set 
(3.24a,b)  can  then  be  written  as 


e  (t)  =  e  (t)  -  —  (e  (t)  -  e  ( t) ) 
J  c  j 


(3.25a) 


0  ( t)  =  0  ( t)  -  - —  0  ( t)  , 

j  s  \j 


(3.25b) 


where  Tr  =  ^j—  =  residence  time  of  fluid  in  core  region. 


T^.  =  j—  =  residence  time  of  fluid  in  shoulder  region. 


and  and  x^  are  defined  as  before.  Note  that  whereas  in  the  previous 
back-mixing  model  the  amount  of  shoulder  water  is  comparable  to  the  amount 
of  core  water,  in  the  present  shear  model  the  shoulder  region  is  much  larger 
than  the  core  region.  Also  note  that  in  the  previous  model  the  shoulder 
receives  all  of  its  tracer  by  way  of  diffusion  from  the  core,  whereas  in  this 
model  the  shoulder  receives  all  of  its  tracer  through  advection  from  the 
overf 1 ow. 
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The  solutions  to  (3.25)  are  determined  straightforwardly , 


®o(t) 

V(t)  =  T^T 


*j 


(3.26a) 


e  (t)  = 
CJ 


«Jt) 


4r  e  (t) 

U  +  c  /  o  \ 
1'  +  6  1  ♦  6  (i  ♦  6j 

CCS 


(3.26b) 


where 


C  * 


and  5$  .  ^ 


However,  we  must  account  for  the  fact  that  once  freon  has  filled  the  core, 

the  fluid  there  mixes  with  freon-free  water  until  that  time  later  when  freon 

also  fills  the  shoulder.  The  time  over  which  this  freon-free  mixing  occurs 

is  Tr  -  T  =  aTr.  The  term  in  parenthesis  in  (3.2t>b)  represents  the 
s  c 

shoulder  concentration  that  mixes  with  the  core.  It  should  therefore  be 
delayed  in  time  by  the  increment  ATr, 


e  ( t) 

•c.(t'  -T-r^ 


,90(t  *  lTr\ 

1  W  6 - ] 


(3.27) 


The  degree  to  which  the  core  speed  and  shoulder  speed  vary  can  be  taken 

as  a  parameter  in  the  model.  Specifically  let  T  =  bT  (b  _>  1).  Also,  we 

s  s  c 
define  the  ratio  —  '  a,  which  measures  the  discrepancy  in  decay  times 

C  a 

between  the  core  and  shoulder.  It  follows  that  =  (—  )«„.  Dropping  the 

S  a  C  J 

subscripts  on  the  core  variables,  the  expression  (3.27)  then  becomes 
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(3.28) 


This  represents  the  core  output  after  the  first  section  (j=l)  and  must 
In  turn  be  used  as  the  input  into  the  next  core.  A  similar  coupling  occurs 
with  the  shoulder,  but  note  that  the  shoulder  evolves  independently  of  the 
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core.  After  n  sections  the  shoulder  concentration  with  which  the  core 
mixes  is 


e.i  t  -  n(  8  -  l)Tr) 

V  - n - ~ 

J‘"  (1  *  (-)s) 

a 


wnere  n(  b  -  1 ) Tp  is  the  discrepancy  in  advective  time  between  the  core  and 
shoulder.  It  follows  that  the  concentration  of  the  core  after  n  sections  is 


v  (t)  *  'rH1  eoUi 

j=n 


n 

v 

u 

i=l 


n-l i-1) 


e0(t-i(8-l)Trl 

(1  ♦  (-16)1 
a 


(3.29) 


Mote  that  expression  (3.^9)  involves  no  simplifying  assumption  concerning  the 
overflow  concentration,  as  was  necessary  in  the  back-mixing  model  (the  reason 
for  this  is  the  governing  equations  (3.2ba,b)  are  algebraic,  and  the  task  of 
coupling  sections  together  remains  manageable). 

As  was  the  case  in  the  back-mixing  model,  the  core  ratio  in  the  shear 
model  is  only  altered  significantly  when  n  is  large  ( 1 . e .  the  process  must 
occur  throughout  the  length  of  the  current).  As  the  magnitude  of  b  is  made 
larger  the  discrepancy  in  ratio  between  the  core  and  shoulder  increases,  and 
the  mixing  then  causes  a  greater  ratio  decay  in  the  core.  However,  b  even¬ 
tually  becomes  so  large  that  for  the  majority  of  the  time  the  core  mixes  with 
freon-free  water.  For  this  range  of  B  then  the  ratio  decay  decreases  with 
increasing  b.  The  effect  of  making  a  larger  is  to  cause  the  shoulder  con¬ 
centration  to  decay  more  slowly  relative  to  the  core.  This  means  simply  that 
tnere  is  more  freon  with  which  to  influence  the  core,  so  the  ratio  decay  of 
the  core  will  be  more  pronounced. 

Figure  3.14  illustrates  these  effects  on  the  ratio  decay  of  the  core. 
The  figure  plots  the  freon  ratio  as  predicted  from  ( 3.^9)  at  t  =  yr 
(i.e.  the  year  the  data  was  collected)  versus  the  magnitude  of  b.  Note  that 
when  i  is  made  greater  than  1  not  only  is  the  ratio  lowered,  but  the 
value  of  b  correspondi ng  to  the  maximum  decay  is  shifted.  As  s  ►  >■  all 
the  curves  will  converge  to  the  input  ratio.  For  the  case  when  b  =  1  the 
curves  also  maintain  tne  value  of  the  input.  In  this  case  the  core  and 
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Figure  3.14:  Three  different  examples  showing  the  value  of  the  core  ratio 
after  a  number  of  sections,  as  a  function  of  the  discrepancy  in  speeds 
between  the  core  and  shoulder. 


shoulder  progress  as  one  unit,  surrounded  by  a  freon-free  reservoir  (hence  no 
ratio  decay).  Note  that  this  is  the  same  set  up  as  the  infinite  reservoir 
case  of  the  back-mixing  model. 

The  manner  in  which  the  shear  model  is  applied  to  the  data  to  determine 
age  is  exactly  the  same  as  was  done  previously  with  the  back-mixing  model, 
only  in  this  case  there  is  an  additional  parameter.  If  a  is  maintained  at 
a  fixed  value,  then  for  each  e  the  matching  process  described  earlier  gives 
an  estimate  of  the  mixing  coefficient  tc  and  the  core  speed.  As  expected, 
the  core  speed  estimate  will  at  first  increase  with  increasing  b  but  even¬ 
tually  fall  back  to  the  freon-free  value  of  2  cm/sec.  Similar  curves  can  be 
generated  for  different  values  of  a.  Three  such  curves  are  shown  in  Figure 
3.15  (T  is  set  equal  to  .2  yr  as  before)  and  it  is  evident  how  successive 
curves  with  larger  a  take  longer  to  reach  their  peak.  As  seen,  it  is  pos¬ 
sible  to  obtain  core  speeds  in  the  range  of  6-8  cm/sec. 

Figure  3.16  shows  two  examples  of  the  cross-stream  variation  in  F-12 
(at  50  °W)  and  the  strength  of  the  mixing  as  a  function  of  core  speed.  It  is 
the  analog  to  Figure  3.11.  In  the  present  case  the  cross-stream  variation 
increases  with  larger  core  speed:  this  is  in  contrast  to  the  previous  back- 
mixing  case  in  which  it  decreases.  Also  recall  that  stronger  mixing  is 
required  at  larger  core  speeds  in  the  previous  case.  Here  the  mixing  remains 
nearly  constant.  (In  each  example  the  diffusive  times  are  larger  than  the 
residence  times  as  required.) 

Although  the  core  speed  is  sensitive  to  the  value  of  b,  the  corre¬ 
sponding  shoulder  speed  stays  in  the  range  of  1-2  cm/sec.  This  can  be  under¬ 
stood  as  follows.  The  shoulder  water  mixes  with  a  freon-free  reservoir,  thus 
if  it  flows  at  2  cm/sec  then  at  50 °W  its  ratio  will  match  that  which  is 
observed  in  the  data  for  the  core.  However,  the  core  has  to  decay  by  way  of 
mixing  with  water  of  a  lower  ratio.  This  puts  an  upper  bound  on  the  shoulder 
of  2  cm/sec.  Recall  that  in  the  back-mixing  model  the  alongstream  variation 
in  core  ratio  was  nearly  independent  of  core  speed.  For  the  same  reasons 
tnis  is  also  true  in  the  shear  model. 


Core  speed  (cm/sec) 


Figure  3.16:  Graph  showing  the  mixing  and  cross-stream  variation  in  concen 
tration  corresponding  to  the  different  core  speeds,  as  in  Figure  3.11,  for 
the  shear  model.  Two  examples  with  different  values  of  a  are  compared. 


A  Comparison  of  the  Back-mixing  and  Shear  Models 

For  Doth  of  the  above  boundary  current  models  there  is  a  whole  range  of 
possible  solutions,  each  one  associated  with  a  different  value  of  the  core 
speed  and  mixing  coefficient  x  It  is  informative  to  compare  a  solution 
from  each  model  with  regard  to  what  is  known  about  the  DWBC  and  its  associ¬ 
ated  freon  signal.  This  can  help  us  better  understand  the  distinguishing 
aspects  of  the  two  models. 

It  is  obvious  that  we  do  not  wish  to  consider  solutions  with  a  small 
core  speed,  as  such  a  value  does  not  seem  to  be  representative  of  the  DWBC. 
(This  was  the  point  in  developing  the  two  models.)  However  we  also  cannot 
use  the  large  core  speed  extremes  in  either  of  the  cases,  for  a  different 
reason.  In  the  back-mixing  model,  in  the  limit  of  small  y,  the  concentra¬ 
tion  of  the  shoulder  water  approaches  that  of  the  core  so  that  there  is  no 
cross-stream  variation  in  freon,  which  is  certainly  not  the  case  in  the  data 
(Figure  3.7).  In  the  large  8  limit  of  the  shear  model  the  opposite  hap¬ 
pens,  and  tne  concentration  of  the  shoulder  water  approaches  zero.  This  is 
unacceptable  as  well.  In  fact,  there  is  really  only  a  very  small  range  of 
core  speed  for  which  the  solutions  are  reasonable  physically. 

The  two  examples  which  are  compared  are  the  5  cm/sec  solution  from  the 
back-mixing  model  and  the  6  cm/sec  (a  =  36)  solution  from  the  shear  model. 
Referring  to  Figure  3.6  of  the  back-mixing  model,  we  are  free  to  specify  the 
value  of  r  ,  i.e.  the  width  of  the  freon  signal  (r  is  constrained  by 
the  fact  that  the  core  and  shoulder  cross-sectional  areas  are  the  same).  The 
width  is  also  freely  specified  in  the  shear  model.  (There  r£  must  satisfy 
the  relation  r  <<  r  ,  Figure  3.13).  We  consider  a  variety  of  widths  and 

V  J 

snow  how  the  two  solutions  vary  accordingly. 

Recall  that  each  solution  was  determined  by  matching  the  value  of  the 
core  concentration  and  core  ratio  to  the  data.  This  was  accomplished  by 
evaluating  the  core  speed  (the  value  of  n)  and  the  mixing  time  scale  x  . 
^hen  written  entirely  in  terms  of  the  radial  distances,  the  mixing  time  scale 

rcrs  1 

taxes  the  form  tc  =  (— ^ — )—  .  It  is  therefore  evident  tha*  n  specifying 
tne  *idth  (and  r. )  ,  this  determines  both  the  transport  of  the  cur- 

and  the  value  of  the  diffusivity  <  .  Figure  3.17  plots  the  transport 
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Figure  3.17:  Transport  of  the  DWBC  associated  with  the  value  of  diffusivity 
as  predicted  by  the  two  boundary  current  models. 
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versus  <  for  different  choices  of  the  width.  (Recall  that  the  vertical  dif- 
fusivity  v  is  related  to  <  by  (2.28).)  Table  3.1  shows  the  corresponding 
values  of  the  radial  distances  r  ,  r  and  grid  spacings  ar  ,  tr  .  It  is 

v  j  L  o 

seen  that  both  models  produce  comparable  transports,  however  we  can  use 
another  consistency  check  with  the  data  to  pin  down  which  values  are  most 
realistic  for  each  model. 

Consider  again  the  trends  in  the  cross-stream  variation  of  freon  for 
the  two  models  (Figures  3.11  and  3.16).  For  the  two  examples  here,  the  back- 
mixing  model  predicts  the  shoulder  concentration  to  be  ~7l)  percent  that  of 
the  core;  the  shear  model  predicts  20  percent.  For  the  freon  section  of  Fig¬ 
ure  3.7,  as  with  any  such  section  across  the  DWBC,  a  Gaussian- type  distribu¬ 
tion  is  associated  with  a  radial  cut  extending  from  the  core.  For  purposes 
of  comparison,  a  Gaussian  was  fit  to  the  horizontal  cut  for  the  section  in 
the  figure.  We  take  the  average  value  of  the  cut  from  r  =  0  to  r  =  rQ  as 
a  measure  of  the  core  concentration,  and  the  average  from  r  =  r  to  r  =  r 
as  the  shoulder  concentration.  Figure  3.18  shows  the  comparison  of  these 
concentrations  to  the  corresponding  model  predictions  for  a  range  of  trans¬ 
ports  (i.e.  range  of  widths).  It  shows  that  the  back-mixing  model  gives  bet¬ 
ter  agreement  for  smaller  transports  (~5  Sv) ,  whereas  the  shear  model  does  so 
for  larger  transports  (~15  Sv) .  (A  smaller  width  clearly  means  a  larger  per¬ 
cent  concentration  for  the  shoulder,  as  this  means  the  shoulder  region  is 
closer  in  proximity  to  the  peak  of  the  Gaussian.) 

We  can  now  more  fully  appreciate  how  the  scenarios  of  the  two  boundary 
current  models  differ.  The  current  in  the  back-mixing  model  is  characterized 
by  thin,  uniform  flow  with  a  small  transport.  High  concentrations  of  tracer 
are  found  in  the  water  directly  surrounding  the  flow.  In  the  shear  model  the 
boundary  current  is  much  more  broad  with  more  intense  flow  at  the  core,  and 
has  a  large  transport.  Only  a  small  amount  of  tracer  accumulates  in  the 
water  just  outside  of  the  current.  These  differences  are  depicted  in  Figure 
3.19.  In  the  former,  the  freon  ratio  of  the  core  decays  because  of  a  diffus¬ 
ive  mechanism;  in  the  latter  it  decays  because  of  an  advective  mechanism 
(although  diffusion  is  still  important). 

At  the  present  time  there  is  some  uncertainty  as  to  what  constitutes  a 
representative  value  of  the  core  speed  of  the  DWBC,  as  well  as  a 
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TABLE  3.1:  Magnitudes  of  the  radial  distances  and  grid  spacings  (Figures  3.6 
and  3.13)  associated  with  the  different  values  of  transport  in  Figure  3.17, 
for  the  two  boundary  current  models. 
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Figure  3.18:  Predicted  concentrations  relative  to  those  measured  from  Fig¬ 
ure  3.7  for  the  core  and  shoulder  regions,  for  different  values  of  the  DWBC 
transport,  (a)  Back-mixing  model . 
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Figure  3.19:  The  core  and  shoulder  regions  (denoted  by  thick  lines)  for  an 
example  from  each  of  the  two  boundary  current  models  superimposed  on  the  F-12 
section  of  Figure  3.7.  (a)  Back-mixing  model  case  In  which  the  DWBC  trans¬ 

port  Is  5  Sverdrups. 
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(b)  Shear  model  case  in  which  the  DWBC  transport  is  15  Sverdrups. 
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representati ve  value  of  the  transport.  The  DWBC  is  known  to  pulsate  in  time 
(for  example  Richardson,  1977)  as  well  as  meander  up  and  down  the  slope 
(Luyten,  1977)  which  complicates  the  determination  of  these  quantities. 
Estimates  of  the  transport  vary  from  values  as  small  as  4  Sverdrups  (Pierce, 
1986)  to  those  as  large  as  24  Sverdrups  (Richardson,  1977),  although  the 
former  is  a  synoptic  estimate  and  the  latter  includes  slope  water  as  shallow 
as  200  m.  More  of  the  estimates  are  in  the  range  of  8-12  Sverdrups  (for 
example  Joyce  et  al . ,  1986;  Hogg,  1983;  Richardson  and  Knauss,  1971). 

In  regard  to  the  core  speed,  numerous  direct  measurements  have  been 
made.  Instantaneous  values  have  been  recorded  as  large  as  50  cm/sec 
(Richardson,  1977),  although  mean  speeds  appear  to  be  more  in  the  range  of 
5-8  cm/sec  (Luyten,  1977;  Richardson,  1977).  Jenkins  and  Rhines  (1980) 
found  a  mean  DWBC  core  speed  of  21  cm/sec  near  the  Blake-Bahama  outer  ridge, 
but  tnis  region  is  characteri zed  by  strong  convergence  of  the  isobaths  which 
should  accelerate  the  flow. 

Both  of  the  boundary  current  solutions  mentioned  here  fall  within  this 
range  of  transports  and  core  speeds.  Although  neither  of  the  models  is  en¬ 
tirely  satisfactory  by  itself,  they  do  suggest  that  to  one  extent  or  another 
both  of  the  freon  ratio  decay  mechanisms  may  be  present  in  the  DWBC.  For  a 
transport  of  5  Sverdrups  in  the  back-mixing  model  and  15  Sverdrups  in  the 
shear  model,  the  diffusivities  so  predicted  are  in  the  range  of  2-4  x 

c  p  p 

10  cm  /sec  for  <  (20-40  cm  /sec  for  v) ,  which  is  in  reasonable  agreement 
with  the  estimate  from  the  flux  balance  calculation  of  chapter  two. 


Time  Dependent  Transfer  Model 

In  the  previous  sections  we  have  investigated  various  advection/diffu- 
sion  schemes  which  influence  the  signal  of  freon  being  carried  by  the  core  of 
the  DWBC.  The  boundary  current  represents  the  primary  source  of  freon  into 
the  deep  layer  in  this  part  of  the  ocean,  and  it  is  of  interest  to  study  the 
process  by  wnich  the  freon  eventually  fills  the  entire  layer  (throughout  the 
various  basins).  Uf  the  components  of  mean  circulation  in  the  interior  abys¬ 
sal  western  North  Atlantic  are  regions  of  closed  circulation,  and  here  we 
comment  on  the  penetration  of  freon  into  the  cyclonic  recirculation  gyre  of 
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the  Gulf  Stream,  or  the  Northern  Recirculation  Gyre  (NRG)  as  it  has  come  to 
be  called. 

OCEANUS  cruise  134  (part  of  the  Abyssal  Circulation  Experiment  in  1983- 
84)  collected  water  sample  data  in  the  region  of  the  NRG,  and  it  included 
several  crossings  of  the  DWBC  which  nicely  show  the  signal  of  F- 12  associated 
with  the  current.  The  entire  cruise  track  is  shown  in  Figure  3.20.  The  ver¬ 
tical  sections  of  F - 12  and  potential  temperature,  for  lines  1-4,  are  present¬ 
ed  in  Figure  3.21.  Note  in  each  section  that  there  are  two  distinct  regions 
in  which  the  freon  is  higher  near  the  inshore  edge  than  on  the  same  tempera¬ 
ture  surfaces  further  in  the  interior.  A  strong  upper  signal  appears  in  the 
potential  temperature  range  4°-6°C,  which  is  interesting  in  its  own  right.  A 
water  property  analysis  suggests  that  this  water  is  of  Librador  Sea  origin, 
yet  it  is  not  in  the  temperature  range  of  classical  Labrador  Sea  Water  (which 
is  ~3.5*C,  Talley  and  McCartney,  1982).  This  feature  warrants  further  inves¬ 
tigation,  although  it  is  not  addressed  in  this  stucty. 

Deeper  in  the  water  column,  against  the  slope,  is  a  somewhat  weaker 
signal  associated  with  the  DWBC.  As  explained  in  chapter  two  it  is  believed 
that  the  eastward  flow  of  the  NRG  pulls  a  plume  of  tracer  from  this  DWBC 
signal.  If  this  is  indeed  the  case  then  the  offshore  signal  from  the  gyre 
should  become  more  distinct  from  the  DWBC  signature  as  one  progresses  from 
section  1  to  4.  Evidence  for  such  a  separation  is  clear  in  Figure  3.21. 

The  difference  between  tracers  such  as  salt  and  oxygen  (whose  charac¬ 
teristics  within  the  NRG  were  discussed  in  chapter  two)  and  freon  is  that  the 
source  of  freon  has  only  just  been  "turned  on,"  so  we  are  presently  in  the 
midst  of  the  spin  up  process.  This  presents  an  ideal  situation  in  the  ocean 
to  monitor  the  accumulation  of  a  passive  tracer  into  a  closed  circulation, 
knowing  that  initially  there  was  no  freon  present.  In  this  section  we  make 
use  of  the  machinery  that  was  developed  in  chapter  two  regarding  the  flux  of 
tracer  into  the  NRG,  and  incorporate  a  time-dependent  source.  The  data  to 
which  these  ideas  are  applied  is  a  single  synoptic  view  of  the  region.  It 
would  be  beneficial  to  add  to  both  its  spatial  and  temporal  resolution,  thus 
allowing  for  a  more  extensive  investigation  into  the  accumulation  process. 
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Figure  3.20:  Station  positions  of  0C134,  June/July  1983.  Sections  1  through 


4  are  marked. 
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Figure  3.21:  Vertical  sections  of  F - 1 2  (p-moles/kg  x  10)  and  potential  tem 
perature  (  Cj  corresonding  to  sections  1-4  of  Figure  3.20. 
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A)  Boundary  Input  of  Freon 

As  Is  the  case  with  oxygen,  the  deep  layer  in  this  region  is  character¬ 
ized  by  higher  concentrations  of  freon  than  those  found  In  the  layers  above 
and  below.  This  means  that  we  must  consider  the  effects  of  vertical  mixing. 
In  particular,  we  will  apply  the  diffusive  transfer  model  developed  earlier 
to  the  distribution  of  F-12.  Recall  that  the  input  to  the  transfer  model 
corresponds  to  the  value  at  tne  core  of  the  DWBC  near  50 °W.  The  time  history 
of  the  F-12  core  is  provided  by  the  boundary  current  models  discussed  above 
(although  it  is  seen  that  there  is  more  than  one  possible  input  curve). 

Before  making  use  of  this  information  we  first  examine  the  case  of  an  arbi¬ 
trary  source  in  the  context  of  the  transfer  model. 

The  only  change  that  must  be  made  in  the  model  is  to  allow  the  input 
value  to  vary  in  time.  Thus  the  appropriate  equations  are  (2.25)  where  $(t) 
replaces  The  solutions  can  still  be  obtained  easily,  and  are  of  the  form 


rf(t)  =  ^(tle  1  ♦  d2(t)e  2 


(3.30) 


Note  the  similarity  of  (3.30)  to  (2.8).  Recall  that  in  the  constant  input 

t_ 

Ti 

case  the  transient  mode  proportional  to  e  was  negligible.  In  the 
present  case  the  two  amplitudes  4j(t)  and  grow  In  time,  yet  It  is 

still  true  (over  the  time  period  being  considered)  that  the  term  proportional 
t_ 

T1 

to  e  can  be  Ignored.  Accounting  for  the  fact  that  there  Is  a  slight 
non-zero  background,  the  expression  for  the  level  of  freon  in  the  gyre  center 
is  then  given  approximately  by 
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A  ♦  1  ♦  e  -  /  tT  *  1 
4  •  (constant)  background  concentration. 

It  Is  of  Interest  to  consider  how  the  F-11:F-12  ratio  varies  In  the 
center  of  the  gyre.  Using  (3.31),  the  expression  for  the  ratio  Is 
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♦  .(t),  t)  »  Input  concentrations  of  F-ll,  F-12 

dp  d2  -  background  concentrations  of  F-ll,  F12. 

Mute  the  similarity  between  (3.32)  and  (3.8),  as  the  bracketed  term  In  both 
expressions  measures  the  effect  of  non-zero  background  on  the  ratio.  In  the 
present  case  the  quotient 
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must  be  small  for  both  freons  In  order  to  Ignore  the  background. 
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The  quantity  t  measures  how  quickly  the  system  would  spin  up  If  the 
Input  were  to  become  steady  (l.e.  It  Is  just  the  spin  up  time  for  the  oxygen 
case).  When  t  is  large  (3.33)  Is  approximated  by 


i 


(3.34) 


♦  ( t* ) dt* 


This  limit  applies  when  either  or  is  large  (t  must  be  large  In 


either  case).  It  happens  that  regardless  of  which  one  Is  specified  as  such 
the  term  in  brackets  >  °°.  This  means  that  even  for  a  large  Input  and  small 
background  the  ratio  Is  altered  by  the  background  concentration  (because  It 
takes  so  long  for  tracer  to  accumulate  In  the  gyre). 

In  the  opposite  case  when  t  Is  small  (3.33)  reduces  to 


(3.35) 


In  order  for  this  to  be  true  either  r  must  be  small  or  both 


and  tc 
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have  to  be  small.  When  the  former  Is  true  the  term  In  brackets  >  • ,  but  when 

the  latter  applies  the  quotient  approaches  a  finite  value.  Thus  In  the  first 

Instance  the  background  contribution  Is  substantial,  but  In  the  second 

Instance  It  will  become  negligible  If  In  time  the  Input  concentration  becomes 

significantly  larger  than  that  of  the  background.  (This  makes  sense  In  that 

when  >  0  all  the  tracer  diffusing  off  the  boundary  lonedlately  gets 

fluxed  vertically  out  of  the  deep  layer,  but  when  tq  and  v  »  0  tracer 

D  <1 

rapidly  accumulates  In  the  gyre.) 

If  It  Is  such  that  the  background  concentration  can  be  ignored  then 
Rc(t)  "  Kt(t).  Rote  that  R^t)  has  the  exact  same  form  as  the  expression 
for  the  freon  ratio  In  the  overflow  basin  (l.e.  the  ratio  formed  from  (J.<J)J. 
Even  though  diffusion  Into  the  gyre  Is  a  different  process  than  that  which 
occurs  in  the  overflow  basin,  because  In  both  instances  the  mixing  Is  assumed 
to  occur  Instantaneously  throughout  the  region  It  turns  out  that  the  two  pro¬ 
cesses  are  described  by  the  same  type  of  equation.  This  means  that  the  same 
limits  concerning  the  freon  ratio  that  were  discussed  In  the  context  of  the 
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overflow  apply  here  (only  with  the  diffusive  time  t  replacing  the  residence 
time).  In  particular.  If  the  Input  time  scale  varies  more  slowly  than  t  , 
then 


4»i(t) 

V*1 


(3.36) 


whereas  In  the  opposite  case  when  the  Input  time  scale  varies  more  quickly 
than  t , 
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B)  F - 1 2  Level  In  the  MRii 

In  chapter  two  when  the  numerical  model  and  transfer  model  were  applied 
to  the  distribution  of  oxygen  It  was  determined  that  roughly  40  percent  of 
the  available  oxygen  from  the  OWBC  had  accumulated  In  the  NRG.  We  are  now  In 
a  position  to  obtain  the  analogous  Information  concerning  F- 12  (or  F-ll).  It 
Is  just  a  matter  of  applying  (3.31)  to  determine  the  level  In  1983  using  the 
values  for  the  unknowns  suggested  by  the  data.  In  regard  to  the  OMBC  Input 
value  we  consider  only  the  back-mixing  case,  l.e.  (3.22a). 

As  was  done  with  oxygen,  a  western  Intensified  gyre  Is  considered,  and 
the  value  of  the  background  Is  taken  to  be  the  average  concentration  at 
1250  m  above  the  deep  freon  maximum.  In  light  of  the  discussion  of  the 
previous  section  It  seems  that  we  should  use  for  the  Input  a  boundary  current 
solution  that  has  a  core  speed  near  5  cm/sec.  As  It  turns  out  however  we  can 
use  any  of  the  solutions,  for  they  all  result  In  nearly  Identical  predictions 
for  the  F-12  level  In  the  gyre  versus  time.  The  reason  for  this  Is  that  the 
F-12  level  (at  50*W)  In  the  boundary  current  Is  Itself  nearly  the  same  for 
each  case.  When  the  core  speed  Is  small  the  level  In  the  current  Is  zero 
over  a  large  extent  of  time,  but  when  the  core  speed  Is  large  the  Increased 
mixing  that  Is  required  keeps  the  level  near  zero  for  almost  the  same  amount 
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of  time  (and  keep  In  mind  that  all  the  curves  are  constrained  to  reach  the 
same  value  In  1983).  Although  this  result  may  seem  surprising  remember  that 
for  small  levels  a  small  change  In  concentration  can  mean  a  significant 
change  In  ratio.  It  Is  the  case  that  after  1983  the  curves  will  begin  to 
separate,  so  that  subsequent  freon  measurements  In  this  region  of  the  DWBC 
would  help  Indicate  If  the  5  cra/sec  solution  Is  In  fact  the  correct  one. 

The  level  of  F-12  that  Is  predicted  for  the  gyre  agrees  closely  with 
that  measured  In  the  data.  However,  this  level  Is  hardly  distinguishable 
from  the  small  amount  of  background  present.  The  model  predicts  that  the 
amount  of  F-12  in  the  gyre  Is  only  ~7  percent  of  that  which  Is  available  from 
the  boundary  as  of  1983  (the  available  amount  Is  one-half  the  concentration 
of  the  core).  The  data  says  that  the  level  Is  “9  percent.  Figure  3.22  com¬ 
pares  the  gyre  level  for  oxygen  and  freon.  That  the  F-12  level  Is  only 
around  10  percent  (versus  40  percent  for  02)  Is  consistent  with  the  fact 
that  the  spin  up  time  f  Is  In  fact  longer  than  the  freon  Input  time  scale 
(l.e.  the  gyre  cannot  respond  quickly  enough  to  changes  In  the  Input).  Note 
that  this  also  means  the  limit  (3.37)  would  apply  for  the  gyre  ratio. 

However  the  level  Is  so  close  to  the  background  that  the  quotient  (3.33)  is 
large  and  thus  causes  Rc( t )  to  deviate  significantly  from  t) .  In 
time,  although  the  level  may  remain  near  10  percent  of  the  available  amount. 
It  will  nonetheless  increase  substantially  above  the  background 
concentration. 

Conclusions 

The  process  by  which  freon  accumulates  In  the  Northern  Recirculation 
Gyre,  the  freon  having  originally  come  from  the  surface  water  In  the 
Norweglan-Greenland  Sea,  Is  quite  complicated.  It  has  been  described  here 
very  crudely  as  being  composed  of  three  stages.  In  the  first  stage  the  sur¬ 
face  water  In  the  N-G  Sea  sinks  to  mid-depths  where  it  undergoes  continual 
mixing  for  roughly  10  years  before  overflowing  a  sill  and  forming  a  boundary 
current.  It  has  been  assumed  that  the  major  contributor  to  the  DWBC  Is  the 
Iceland- Scotland  overflow,  but  this  point  remains  unclear.  The  second  stage 
corresponds  to  the  time  during  which  the  freon  travels  In  the  boundary  cur¬ 
rent,  making  its  way  around  the  Grand  Banks  to  the  region  near  the  NRG. 
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Figure  3.22:  Average  sections  of  freon  and  oxygen  through  the  DWBC  and  NRG 
as  In  Flaure  2.21.  The  concentration  of  the  deep  anomaly  Is  plotted,  normal 
lied  by  its  value  at  the  boundary. 


Finally,  in  the  last  stage,  the  freon  diffuses  from  the  DWBC  and  becomes 
entrained  into  the  gyre. 

For  the  case  of  oxygen  which  was  discussed  in  chapter  two  it  was  only 
necessary  to  consider  the  last  stage,  but  unlike  oxygen  the  level  of  freon 
in  the  ocean  is  increasing  at  a  rapid  rate.  It  is  assumed  that  the  surface 
water  concentration  of  freon  versus  time  in  the  N-G  Sea  is  known  with  some 
certainty.  This  information  represents  the  boundary  condition  for  the  three 
coupled  components.  Because  there  is  a  discrepancy  in  the  Increase  of  F-ll 
versus  F-12  this  means  that  the  F-ll :F- 12  ratio  In  the  surface  water,  and 
consequently  the  water  which  overflows  the  sill,  changes  (increases)  continu¬ 
ally  in  time. 

In  the  description  of  the  second  stage  two  different  schemes  were  pre¬ 
sented  in  which  the  core  of  the  current  mixes  with  surrounding  water  of  a 
lower  F-ll :F-12  ratio.  Extreme  cases  in  both  models  predict  that  the  ratio 
of  the  core  remains  unchanged,  and  this  is  in  line  with  previous  ideas  in 
which  the  boundary  current  core  speed  is  estimated  using  the  value  of  the 
ratio.  However  there  are  many  other  solutions  in  the  two  models  which  alter 
this  prediction  by  varying  amounts.  At  present  each  of  these  solutions  pro¬ 
duces  an  accurate  estimate  of  the  F-12  level  in  the  NRG,  as  the  level  is 
barely  above  that  of  the  background. 
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Summary 

The  diffusive  transfer  of  properties  from  the  Deep  Western  Boundary 
Current  (DWBCJ  into  the  Northern  Recirculation  Gyre  of  the  Gulf  Stream  was 
Investigated  in  some  detail.  An  advective-diffusive  numerical  model  was  im¬ 
plemented  to  aid  in  the  study.  The  model  consisted  of  a  specified  gyre  cir¬ 
culation  located  alongside  a  boundary  current,  both  flows  being  steady. 
Tracer  was  injected  continuously  into  the  boundary  current  upstream  of  the 
gyre,  and  the  manner  in  which  the  tracer  eventually  fills  the  domain  was 
monitored. 

The  entrainment  of  tracer  into  the  elliptical  model  gyre  from  the 
boundary  current  is  characterized  by  a  plume  encircling  the  gyre,  while  spi¬ 
ralling  in  towards  its  center.  The  extent  of  the  spiral  can  be  related  to 
the  characteristics  of  the  gyre  velocity  field.  In  particular,  where  the 
flow  is  strong  and  the  cross-stream  shear  is  large  the  spiral  is  hardly 
noticeable,  whereas  in  the  weak  flow  where  the  cross-stream  shear  is  small 
the  spiral  is  more  pronounced.  The  inward  progression  of  the  spiral  eventu¬ 
ally  halts  as  the  plume  penetrates  the  interior,  faster  flowing  portion  of 
the  gyre.  It  is  this  region  which  becomes  homogenized  in  the  steady  state, 
provided  the  lateral  diffusivity  is  small  enough.  As  the  diffusivity  is 
increased,  the  manner  of  entrainment  is  altered  and  the  pool  of  homogenized 
tracer  shrinks  in  size. 
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A  simple  box-model  representation  of  the  processes  occurring  in  the  nu¬ 
merical  model  was  solved  analytically,  which  led  to  a  more  quantitative  un¬ 
derstanding  of  some  of  the  numerical  results.  Two  different  types  of  tracers 
were  studied,  one  of  which  was  subject  to  lateral  diffusion  only,  the  other 
subject  to  lateral  and  vertical  diffusion.  It  was  found  that  the  presence  of 
vertical  mixing  has  little  effect  on  the  homogenization  that  occurs  in  the 
gyre,  but  has  a  profound  effect  on  the  amount  of  tracer  which  gets  entrained 
and  how  quickly  the  steady  state  is  reached.  With  no  vertical  mixing,  the 
level  of  the  homogenized  pool  in  the  gyre  is  independent  of  the  strength  of 
the  lateral  mixing  and  is  equal  to  one-half  the  core  concentration  of  the 
boundary  current  at  the  upstream  boundary.  When  vertical  mixing  is  added  the 
gyre  level  decreases  and  also  becomes  dependent  on  the  size  of  the  lateral 
diffusivity.  The  equilibration  time  of  the  gyre  decreases  as  well. 
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Various  results  from  the  numerical  and  box-model  study  were  applied  to 
the  A8CE  tracer  data  set.  The  presence  of  vertical  mixing  appears  to  be  the 
reason  for  a  discrepancy  between  the  deep  lateral  distributions  of  salt  and 
oxygen,  the  latter  being  characterized  by  a  slight  minimum  in  the  region  of 
the  NRG.  The  values  of  diffusivity  which  give  the  best  agreement  between 
data  and  model  concerning  the  level  of  oxygen  in  the  NRG  and  the  associated 
diffusive  fluxes  are  <  ~  106cm2/sec  (lateral)  and  v  ~  10  cm2/sec  (ver¬ 
tical).  Although  seemingly  large,  such  a  value  for  v  results  in  only  a 
moderate  effect  on  the  o*ygen  level  for  an  NRG  which  is  western-intensified. 

The  evolution  of  freon  within  the  NRG  was  studied  as  well,  but  only  in 
the  context  of  the  box-model .  Because  of  the  time  dependency  of  freon  a 
broader  perspective  had  to  be  taken  which  addresses  what  happens  to  the  freon 
prior  to  when  it  diffuses  from  the  DWBC  into  the  NRG.  A  simple  representa¬ 
tion  of  the  overflow  process  in  the  Norwegian-Greenland  Sea  predicts  that 
both  the  freon  concentration  and  freon  ratio  are  altered  by  this  process. 

Two  different  boundary  current  models  show  how  the  ratio  may  be  further 
altered  by  way  of  mixing  in  the  DWBC.  These  mixing  schemes  suggest  a  DWBC 
core  speed  of  5-6  cm/sec.  When  the  prior  history  of  the  freon  is  viewed 
together  with  the  regional  box-model,  it  implies  that  freon  is  only  now 
beginning  to  accumulate  in  the  NRG. 


7\(Y, ,  i,  iifcfc- j ^ .vvVva. 'fc.-Tv  v- «•- » •r-  ■ 


*«ya 

m 

s::-! 


i 


lh 


I  would  like  to  thank,  my  committee  members  Dale  Haidvogel,  Bill  Young, 
Harry  Bryden,  and  Mike  McCartney  for  their  guidance  with  this  work.  Bill 
and  Dale,  by  way  of  numerous  discussions,  gave  me  very  useful  ideas  and 
suggestions.  Bill  Smethie  kindly  allowed  me  to  use  unpublished  freon  data 
and  provided  much  insight  concerning  its  application.  I  would  also  like  to 
thank  Fred  Grosse,  whose  enthusiasm  for  teaching  increased  my  enthusiasm 
for  learning.  My  parents,  Stan  and  Louise,  have  always  given  me  support 
and  confidence.  My  wife  Anne  has  been  patient  and  giving.  I  thank  Mary 
Ann  Lucas  very  much  for  preparing  the  manuscript. 

Lastly,  I  would  like  to  thank  my  advisor.  Nelson  Hogg,  who  constantly 
shared  his  time,  his  ideas,  and  his  experience  with  me.  I  have  learned 
very  much  from  him. 

This  work  was  supported  by  the  office  of  Naval  Research  through 
contracts  N00014-76-C-0197  and  N00014-84-C-0134,  NR  083-400;  and 
N00014-82-C-0019  and  N00014-85-C-0001,  NR  083-004,  and  the  National  Science 
Foundation  through  grant  OCE82-14925. 


/.  v-.-v v\'  .•■•/•.■  v'v  vv'v  .  /jo  v  v  -.-.v- ■.  v  v- . .  .  - 

•\  \  \  v  '.v  • v  v  *.  v  w  *.  *,*  v  *,•  .*  v  ..  %  v  v  vV 

VJVjIijW  rjiAj*  V  ,v* W ^  ^  -  r  *  *  »  r  **  VV  -  * •  *  ■  » V  *  r  m  Jr  » V 


References 


Bullister,  J.  L.,  1984.  Atmospheric  chlorofluoromethanes  as  tracers  of  ocean 
circulation  and  mixing:  measurement  and  calibration  techniques  and 
studies  in  the  Greenland  and  Norwegian  Seas.  Ph.D.  Thesis,  University 
of  California,  San  Diego. 

Freeland,  H.  J. ,  P.  B.  Rhines,  and  T.  Rossby,  1975.  Statistical  observations 
of  the  trajectories  of  neutrally  buoyant  floats  in  the  North  Atlantic. 
Journal  of  Marine  Research,  33,  383-404. 

Harrison,  D.  E.,  1982.  On  deep  mean  flow  generation  mechanisms  and  the  abys¬ 
sal  circulation  of  numerical  model  gyres.  Dynamics  of  Atmospheres  and 
Oceans,  2*  289-325. 

Hogg,  N.  G.,  1983.  A  note  on  the  deep  circulation  of  the  western  North 
Atlantic:  its  nature  and  causes.  Deep-Sea  Research,  30,  945-961. 

Hogg,  N.  G. ,  P.  Biscaye,  W.  Gardner,  and  W.  J.  Schmitz,  Jr.,  1982.  On  the 
transport  and  modification  of  Antarctic  Bottom  Water  in  the  Verna 
Channel.  Journal  of  Marine  Research,  Supplement  to  40 ,  231-263. 

Hogg,  N.  G.,  and  H.  Stommel ,  1985.  On  the  relationship  between  the  deep 
circulation  and  the  Gulf  Stream.  Deep-Sea  Research,  32,  1181-1193. 

Hogg,  N.  G.,  R.  S.  Pickart,  R.  M.  Hendry,  and  W.  J.  Smethie,  1986.  The 
northern  recirculation  gyre  of  the  Gulf  Stream.  Deep-Sea  Research, 

33,  1139-1165. 

Holland,  W.  R.,  1978.  The  role  of  mesoscale  eddies  in  the  general  circula¬ 
tion  of  the  ocean  --  numerical  experiments  using  a  wind-driven  quasi- 
geostrophic  model.  Journal  of  Physical  Oceanography,  8,  363-392. 

Holland,  W.  R.,  and  L.  B.  Lin,  1975.  On  the  generation  of  mesoscale  eddies 
and  their  contribution  to  the  oceanic  general  circulation.  Parts  I 
and  II.  Journal  of  Physical  Oceanography,  S,  642-669. 

Holland,  W.  R.,  and  P.  B.  Rhines,  1980.  An  example  of  eddy-induced  ocean 
circulation.  Journal  of  Physical  Oceanography,  10,  1010-1031. 

Jenkins,  W.  J. ,  and  P.  B.  Rhines,  1980.  Tritium  in  the  deep  North  Atlantic 
Ocean.  Nature,  286,  877-880. 

Joyce,  T.  M.,  C.  Wunsch,  and  S.  D.  Pierce,  1986,  Synoptic  Gulf  Stream  veloc¬ 
ity  profiles  through  simultaneous  inversion  of  hydrographic  and  acous¬ 
tic  Doppler  data.  Journal  of  Geophysical  Research,  91,  7573-7585. 


Luyten,  J.  R.,  1977.  Scales  of  motion  in  the  Deep  Gulf  Stream  and  across  the 
continental  rise.  Journal  of  Marine  Research,  35,  49-74. 

McCartney,  M.  S.,  L.  V.  Worthington,  and  M.  £.  Raymer,  1980.  Anomalous  water 
mass  distributions  at  55°W  in  the  North  Atlantic  in  1977.  Journal  of 
Marine  Research,  38,  147-172. 

McDowell,  S.,  P.  Rhines, and  T.  Keffer,  1982.  North  Atlantic  potential  vor- 
ticity  and  its  relation  to  the  general  circulation.  Journal  of 
Physical  Oceanography,  12,  1417-1436. 

Musgrave,  D.  L.,  1985.  A  numerical  study  of  the  roles  of  subgyre-scale  mix¬ 
ing  and  the  western  boundary  current  on  homogenization  of  a  passive 
tracer.  Journal  of  Geophysical  Research,  90,  7037-7043. 

Niiler,  P.  P. ,  1982.  The  general  circulation  of  the  oceans.  Proceedings  of 
the  50th  Anniversary  Symposium,  Woods  hole  Oceanographic  Institution. 

Pierce,  S.  D.,  1986.  Gulf  Stream  velocity  structure  through  combined  inver¬ 


sion  of  hydrographic  and  acoustic  doppler  data.  M.  S.  Thesis,  Massa¬ 
chusetts  Institute  of  Technology,  65  pp. 

Rhines,  P.  8.,  1983.  Lectures  in  geophysical  fluid  dynamics.  In:  Fluid 
Dynamics  in  Astrophysics  and  Geophysics,  Norman  R.  Lebovitz,  editor; 
Lectures  in  Applied  Mathematics,  20,  American  Mathematical  Society, 
Providence,  Rhode  Island;  pp.  3-68. 

Rhines,  P.  8.,  and  W.  R.  Young,  1982a.  Homogenization  of  potential  vorticity 
in  planetary  gyres.  Journal  of  Fluid  Mechanics,  122,  347-367. 

Rhines,  P.  B.,  and  W.  R.  Young,  1982b.  A  theory  of  the  wind-driven  circula¬ 
tion.  I.  Mid-ocean  gyres.  Journal  of  Marine  Research,  Supplement  to 
40  ,  559-596. 

Rhines  P.  B.,and  W.  R.  Young,  1983.  How  rapidly  is  a  passive  scalar  mixed 
within  closed  streamlines?  Journal  of  Fluid  Mechanics,  133,  133-145. 

Richardson,  P.  L.,  1977.  On  the  crossover  between  the  Gulf  Stream  and  the 
western  boundary  undercurrent.  Deep-Sea  Research,  24,  139-159. 

Richardson,  P.  L.,  1985.  Average  velocity  and  transport  of  the  Gulf  Stream 
near  55°W.  Journal  of  Marine  Research,  43,  83-111. 

Richardson,  P.  L.,and  J.  A.  Knauss,  1971.  Gulf  Stream  and  western  boundary 
undercurrent  observations  at  Cape  Hatteras.  Deep-Sea  Research,  18, 
1089-1109. 


194 


Roache,  P.  J.,  1972.  Computational  Fluid  Dynamics.  Hermose,  434  pp. 

Rooth,  C.  G. ,  and  H.  G.  Ostlund,  1972.  Penetration  of  tritium  into  the 
Atlantic  thermocline.  Deep-Sea  Research,  19,  481-492. 

Smethie,  W.  M.,  and  S.  Trumbore,  1984.  Chlorofluoromethanes  (F-ll  and  F- 12) 
in  the  western  North  Atlantic  Ocean  and  the  deep  western  boundary 
undercurrent.  Eos,  Transactions,  American  Geophysical  Union,  64,  p. 
1089. 

Smith,  R.,  1982.  Dispersion  of  tracers  in  the  deep  ocean.  Journal  of  Fluid 
Mechanics,  123,  131-142. 

Smolarkiewicz,  P.  K.,  1983.  A  simple  positive  definite  advection  scheme  witn 
small  implicit  diffusion.  Monthly  Weather  Review,  111,  47y-486. 

Swift,  J.  H. ,  1984.  The  circulation  of  the  Denmark  Strait  and  Iceland- 
Scotland  overflow  waters  in  the  North  Atlantic.  Deep-Sea  Research, 

31,  1339-1355. 

Talley,  L.  U.,  and  M.  S.  McCartney,  1982.  Distribution  and  circulation  of 
Labrador  Sea  Water.  Journal  of  Physical  Oceanography,  12,  1189-1205. 

Warren,  B.  A.,  1981.  Deep  circulation  of  the  world  ocean.  In:  Evolution 

of  Physical  Oceanography,  Scientific  Surveys  in  Honor  of  Henry  Stommel , 
B.  A.  Warren  and  C.  Wunsch,  editors;  The  MIT  Press,  Cambridge,  Massa¬ 
chusetts;  pp.  6-41. 

Weiss,  R.  F.,  J.  L.  Bullister,  R.  H.  Gammon,  and  M.  J.  Warner,  1985.  Atmos¬ 
pheric  chlorofluoromethanes  in  the  deep  equatorial  Atlantic.  Nature, 
314,  608-610. 

Whitehead,  J.  A.,  Jr.,  and  L.  V.  Worthington,  1982.  The  flux  and  mixing 

rates  of  Antarctic  Bottom  Water  within  the  North  Atlantic.  Journal  of 
Geophysical  Research,  87,  7903-7924. 

Worthington,  L.  V. ,  1970.  The  Norwegian  Sea  as  a  Mediterranean  basin.  Deep- 
Sea  Research,  17,  77-84. 

Worthington,  L.  V. ,  1976.  On  the  North  Atlantic  circulation.  The  Johns 
Hopkins  Oceanographic  Studies,  6,  110  pp. 

Wunsch,  C. ,  and  B.  Grant,  1982,  Towards  the  general  circulation  of  the  North 
Atlantic  ocean.  Progress  in  Oceanography,  1 1,  1-59. 

Zill,  D.  G. ,  1979.  Differential  Equations.  Prindle,  Webster,  and  Schmidt, 
486  pp. 


195 


APPENDIX  A 


Abbreviations  and  symbols: 

aX,aY  zonal,  meridional  grid  spacing 
at  temporal  grid  spacing 
ADE  approximate  difference  equation 
FT  forward  In  time 

US  upstream  In  space 

CT  centered  In  time 

CS  centered  in  space 

Interior  Scheme 

The  governing  equation  Is, 

jl  e(x,y,t)  ♦  u(x,y)  *  v  «|x,y,t)  =  v  *  <v  e(x,y,t)  (A.l) 

Several  finite-difference  schemes  were  tested. 

1.  CTCS  advection/FTCS  diffusion 

The  diffusive  test  produced  accurate  results.  However,  It  is  a  proper¬ 
ty  of  the  CTCS  scheme  that  advectlve  phase  error  occurs.  Briefly,  this  means 
that  the  phase  speed  of  a  wave  depends  on  its  wave  number,  so  rather  than  a 
disturbance  propagating  as  a  unit  at  the  advectlve  speed  (the  analytical 
result)  it  disperses  into  its  component  as  shorter  waves  lag  longer  ones. 
Figure  A. la  shows  the  result  of  a  zonal  advectlve  test.  There  are  two  ways 
to  reduce  phase  error.  One  is  to  choose  the  parameters  ax  and  At  such  that 
the  courant  number,  Cx  =  uAt/Ax  ,  is  close  to  one.  (Cx  =  1  eliminates 
phase  error  altogether.)  The  other  is  to  increase  the  spatial  resolution. 

The  effect  of  these  adjustments  is  Illustrated  in  Figures  A. lb  and  A.lc. 

Many  numerical  schemes  include  a  stipulation  on  the  parameters  of  the 
ADE  In  order  to  avoid  numerical  Instability  (an  artifact  of  the  ADE  that 
gives  rise  to  an  unstable  solution).  Some  schemes  are  unconditionally 


mm 


Figure  A.l:  (a)  The  effect  of  phase  error  on  a  Gaussian  distribution  of 
tracer  being  advected  by  a  uniform  flow.  The  amplitude  of  the  strongest 
trailing  wave  is  ~2  percent  that  of  the  Gaussian.  Note  that  in  addition  to 
the  dispersion,  phase  error  manifests  itself  through  distortion  of  the  (ini 
tially  symmetric)  distribution. 
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(b)  The  same  distribution  and  flow  field  as  In  (a)  only  the  courant 
number  uat/ax  has  been  made  close  to  one. 
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unstable  as  such.  (An  example  Is  CTCS  for  diffusion,  which  Is  why  FTCS  was 
chosen  even  though  It  Is  a  less  accurate  approximation.)  For  CTCS  advectlon, 
stability  places  an  upper  bound  on  the  size  of  the  courant  number.  For  this 
application  of  the  scheme,  ax  and  ay  cannot  be  made  small  enough  to  keep 
the  phase  error  at  an  acceptable  level  while  still  satisfying  the  stability 
criterion.  As  for  the  other  adjustment  to  minimize  phase  error,  the  velocity 
field  Is  too  complex  (among  other  reasons)  to  keep  the  courant  number  close 
to  one  everywhere  In  the  domain. 

2.  FTUS  advectlon/FTCS  diffusion 

Although  phase  error  is  absent  In  the  lower  order  scheme,  another  phen¬ 
omenon,  implicit  diffusion.  Is  present.  This  error  source  causes  a  distribu¬ 
tion  of  tracer  to  spread  as  If  It  were  being  acted  upon  by  explicit  diffusion 
(Figure  A. 2).  As  with  the  higher  order  scheme.  Increased  spatial  resolution 
reduces  the  error.  However,  It  Is  Impossible  to  keep  the  implicit  diffusion 
less  than  the  explicit  diffusion  without  causing  numerical  Instability, 
within  the  realm  of  desired  parameters. 


It  is  a  property  of  higher  order  advectlve  schemes  to  exhibit  phase 
error  and  those  of  lower  order  to  exhibit  implicit  diffusion  (Smolarkiewicz, 
1983).  The  approximation  that  was  applied  to  (A.l)  makes  use  of  a  lower 
order  advectlve  scheme  developed  by  Smolarkiewicz  (1983)  that  includes  a  sep¬ 
arate  step  to  counteract  the  Implicit  diffusion.  The  procedure  is  based  on 
tne  fact  that  the  form  of  the  Implicit  dlffusivlty  can  be  determined.  (It 
depends  on  the  parameters  In  the  problem.)  Therefore  a  process  that  acts  as 
the  reverse  of  diffusion  can  be  adjusted  to  offset  the  implicit  diffusion. 

The  anti-diffusion  procedure  Is  cast  into  the  form  of  advection  (with  an 
effective  velocity). 

The  composite  scheme  that  was  used  then  consists  of  three  steps: 

1)  FTUS  advection 

2)  FTUS  advection  with  effective  velocity  (i.e.  anti-diffusion) 

3)  FTCS  diffusion 

The  accompanying  stability  criterion  is  (Roache,  1972) 


‘/j.WAV .  .  • :  v-V-V- '.v'riA-.VvVx.'jvVo,- 


Although  the  problem  of  Implicit  diffusion  Is  addressed,  this  scheme 
Is  not  without  errors.  For  example,  artificial  diffusion  occurs  during  the 
second  advectlve  step  (Smolarklewlcz,  19113).  Another  drawback  of  FTUS  for 
advectlon  Is  the  presence  of  a  directionally  dependent  distortion.  This 
error  reveals  Itself  only  when  both  u  and  v  are  non-zero  and  Is  most 
prevalent  when  |u|  and  |v|  are  equal  {see  Figure  A. 3  for  a  discussion). 

Boundary  Conditions 

The  edges  of  the  domain  do  not  correspond  to  solid  boundaries  but 
rather  open  ocean.  There  are  three  flow  regimes  that  exist  along  different 
portions  of  the  boundary. 

A)  Inflow 

Where  the  boundary  current  flows  Into  the  domain  the  concentration  of 
tracer  Is  specified. 

B )  Outflow 

Where  the  boundary  current  flows  out  of  the  domain  diffusion  normal  to 
the  boundary  Is  omitted,  as  Is  anti-diffusion.  Advectlon  and  cross-stream 
diffusion  are  carried  out  as  In  the  Interior.  (Figure  A. 4  shows  results 
from  an  experiment  testing  the  outflow  boundary  condition.) 

C)  Mo  Flow 

The  Inflow  and  outflow  boundary  conditions  are  applied  where  the  bound¬ 
ary  current  velocity  Is  significantly  greater  than  zero  (which  was  chosen  to 
be  greater  than  .1  cm/sec).  Along  the  remaining  part  of  the  boundary  veloci¬ 
ties  are  smaller  than  this  and  are  set  equal  to  zero,  so  that  only  the  diffu¬ 
sion  step  requires  a  boundary  condition.  However,  CS  diffusion  requires  that 
the  concentration  be  known  at  neighboring  grid  points,  which  is  not  possible 
at  the  boundary.  An  extrapolation  boundary  condition  was  developed  which  at¬ 
tempts  to  simulate  open  boundary  diffusion.  It  is  based  on  the  fact  that 
the  normal  gradient  of  tracer  along  the  no-flow  segments  will  be  everywhere 
Inward,  l.e.  there  will  be  an  outward  diffusive  flux. 


Figure  A. 3:  (a)  Gaussian  In  a  uniform  flow  In  which  |  u|  *  |v|  demonstrat¬ 

ing  another  type  of  error  associated  with  the  FTUS  scheme.  (There  Is  neglig¬ 
ible  Implicit  diffusion  because  the  corrective  step  has  been  employed.) 


(b)  Same  conditions  as  in  (a)  only  the  advectlon  is  performed  In  two 
steps  --  each  dimension  separately.  This  eliminates  the  distortion.  (Im¬ 
plicit  diffusion  Is  present  because  the  corrective  step  was  omitted.)  This 
technique  Is  not  feasible,  however,  because  it  must  be  accompanied  by  two 
anti-diffusive  steps,  and  the  resulting  computational  expenses  would  be  too 
high. 


The  simplest  approximation  to  an  open  boundary  consists  of  specifying 
a  constant  flux  at  the  boundary  (a  passive  flux).  This  means  that  the  con¬ 
centration  there  would  evolve  in  time  (an  active  concentration).  Allowing 
the  flux  to  change  in  time  according  to  a  specified  criterion  (an  active 
flux)  forms  an  improved  estimate.  A  natural  way  to  come  up  with  such  a  cri¬ 
terion  would  be  to  make  use  of  what  is  known  about  the  flux  just  inside  of 
the  boundary;  in  particular,  compute  the  trend  of  outward  flux  approaching 
the  boundary  along  lines  perpendicular  to  it,  and  extrapolate  to  compute  the 
flux  at  the  boundary  --  from  which  the  concentration  is  then  determined.  In 
line  with  this,  consider  the  following  procedure  to  be  performed  at  each 
time  step  after  the  interior  solution  has  been  determined. 

i)  Tne  value  of  the  flux  at  a  point  is  defined  in  an  upwind  sense 


where  e1+1  is  the  concentration  adjacent  to  in  the  direc¬ 
tion  of  the  inward  normal,  and  as  is  the  grid  spacing  (equal  to 
ax  or  Ay).  The  value  of  F.  near  the  boundary  must  always  be 
<_  0.  F.j  is  computed  at  the  three  points  prior  to  the  boundary 
along  the  normal  (i  =  1  corresponds  to  the  point  closest  to  the 
boundary) . 

i  1 )  Let  Fq  denote  the  value  of  F..  at  the  boundary.  If  Fj  =  J, 

Fq  Is  set  =  0. 

iii)  If  F1  is  non-zero,  Fq  is  predicted  using  a  three-point 
extrapolation:  Fq  =  3F^  -  3F2  +  Fg. 

iv)  If  the  predicted  Fq  Is  >0,  Fq  is  reset  =  0. 

v)  The  concentration  of  tracer  at  the  boundary  is  then  determined 

from  the  value  of  F  .  (If  the  calculated  concentration  is 

o 

<  0,  then  it  is  reset  =0.) 

Tnis  procedure  was  applied  in  a  test-run  which  consisted  of  having  a 
Gaussian  spot  of  tracer  diffuse  near  a  boundary.  Results  of  the  test  are 
shown  in  Figure  A. 5,  Note  that  as  time  progresses  the  contours  right  next 
to  the  boundary  become  artificially  squished,  indicating  that  the  gradients 


there  are  too  large.  This  problem  can  be  corrected  by  extrapolating  a  dif¬ 
ferent  quantity  other  than  the  flux.  . 

Consider  the  parameter  R^  =  (— - ,  which  is  the  ratio  of  adja¬ 

cent  concentrations  in  the  direction  of  the  Inward  normal  (e  is  a  small  num¬ 
ber  to  prevent  division  by  zero).  In  the  region  near  the  boundary  will 
always  fall  between  zero  and  one  (the  former  means  tracer  is  just  beginning 
to  penetrate  the  area,  the  latter  corresponds  to  no  flux).  If  R^  is  ex¬ 
trapolated  in  the  manner  described  above  in  order  to  calculate  the  boundary 
concentration,  this  results  in  weaker  gradients  than  the  with  flux  extrapola¬ 
tion.  The  reason  for  this  is  that  this  procedure  involves  percents  rather 
than  differences.  For  example,  if  there  were  a  linear  decrease  in  concentra¬ 
tion  towards  the  boundary,  the  flux  prediction  for  e  at  the  boundary  would 
continue  this  trend,  whereas  the  ratio  prediction  would  cause  the  trend  to 
flatten  out,  which  is  more  realistic  of  diffusion. 

Tne  active  flux  condition  that  was  used  in  the  model  then  is  outlined 
as  follows. 

i)  The  value  of  R.  is  computed  at  the  three  points  prior  to  the 

boundary  along  the  normal  (i  =  1  corresponds  to  the  point  closest 
to  the  boundary). 


d 


> 

•A 


ii)  Let  Rq  denote  the  value  of  R^  at  the  boundary.  If  R^  =  0, 
Rq  is  set  =  0. 

iii)  If  Rj^  is  non-zero,  RQ  is  predicted  using  a  three-point 
extrapolation:  RQ  =  3R^  -  3R£  +  Rg. 

i v)  If  the  predicted  R  >  1,  R  is  reset  =  1. 

o  o 

v)  If  the  predicted  RQ  <  0,  RQ  is  reset  =  0. 

vi )  The  concentration  of  tracer  at  the  boundary  is  then  determined 

from  the  value  of  RQ. 

Results  from  the  same  test-run  applied  above,  using  this  scheme,  are  shown 
in  Figure  A. 6.  It  is  seen  that  the  isolines  near  the  boundary  no  longer 
appear  to  get  distorted. 


Corners 

The  four  corners  of  the  domain  lie  in  the  no-flow  region;  the  active 
flux  condition  is  applied  at  each,  with  the  extrapolation  performed  along 
the  diagonal. 


Table  A.l  contains  a  summary  of  the  boundary  conditions. 
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TABLE  A. 1 

BOUNDARY  CONDITIONS 


I. 

First  Step: 

Advectlon 

(1) 

Inflow, 

specify  concentration 

(11) 

outflow, 

— 

(111) 

no-flow. 

II. 

Second  Step: 

Anti-diffusion 

(1) 

Inflow, 

specify  concentration 

(11) 

outflow. 

omit  normal  flux 

(111) 

no-flow. 

III. 

Third  Step: 

Diffusion 

(1) 

Inflow, 

specify  concentration 

(11) 

outflow. 

omit  normal  flux 

(ill) 

no-flow, 

active  flux 

I* 
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I  hr-  distributions  of  tracer  associated  with  the  Northern  Recirculation  Gyre  of  the  Gulf  Stream 
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flow  parameters.  Homogenization  occurs  at  steadv  state,  consistent  with  recently  collected  tracer 
data.  The  presence  of  vertical  mixing  is  considered  in  an  attempt  to  explain  a  difference  between 
salimt \  and  oxvgen  observed  in  the  data.  (  ompanson  of  the  model  results  to  the  oxygen  data 
is  favorable,  and  leads  to  an  estimate  of  the  lateral  and  vertical  diffusivity.  The  time  dependent 
nature  t  ireor.  is  addressed  using  a  coupled  model  of  the  deep  water  overflow  process, 
act c  e.  non  mixing  in  the  1)WB(  ,  and  subsequent  entrainment  into  the  NRG.  Comparison  with  the 
data  shows  that  very  little  freon  has  accumulated  in  the  NRG,  and  that  these  processes  feet 
th<-  Ire  or  IDIreon  1 ratio  as  well. 


..  •..imogenizat  ion  d  tru-  er> 

.  <  •  ui  t  >’•  c  ;j ir .  ul a  t  ion 

c.  Dec  Acstcrr.  Boundary  1  irrent 


t>  i  Opi  ■  ( 


COUTl 


14  »*•«*••»<•»»  PpMNl1 


Approval  for  ptif'U-  ation;  list riLut ion  inhmiuM. 


14  ct**»  rr**  mpudi 

UNCLASSIFIED 

>4  ClM4  fTMp 


